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- XXIX. BAND DRITTES HEFT 1960 


Drehzahlabhangige Selbsterregung des unsymmetrischen Kreisels 


Von R. Grammel 


Kurt von Sanden zum fiinfundsiebenzigsten Geburtstag gewidmet 


1. Einleitung. Ein (unsymmetrischer) Kreisel, also ein beliebig gestalteter Kérper, der sich 
um seinen Schwerpunkt beliebig und (in unserem Falle) unabhangig von einer etwaigen Bewegung 
des Schwerpunktes drehen kann, heift selbsterregt, wenn seine Drehung durch ein inneres Moment 
erzeugt oder geandert wird, d. h. durch ein Moment, dessen Vektor im Kérper fest ist oder in ihm 
in vorgeschriebener Weise wandert. Ein solches Moment kann als Reaktion von ausgeschleu- 
derten Massen entstehen, deren Betrag je Zeiteinheit hier vernachlassigbar klein sein soll gegen- oa 
tiber der Gesamtmasse, damit wir die zeitliche Massenanderung des Kreisels auBer Betracht lassen 
kénnen, wogegen die Ausschleudergeschwindigkeit hinreichend grof sein mu8, um ein merk- 
liches Reaktionsmoment zu erzeugen. Es soll hier nur angedeutet werden, daB ein inneres Mo- 
ment auch durch einen eingebauten Motor hervorgerufen werden kann, der einen kleinen und 
(wie das heute méglich ist) freischwebenden, dynamisch kugelsymmetrischen Rotor im Inneren 
des Kérpers antreibt. 

Fiir den Fall, daB das Moment einer solchen Selbsterregung konstant ist oder in zeitlich vor- 
geschriebener Weise gesteuert wird, habe ich die Theorie in mehreren Arbeiten! entwickelt. Im 
Folgenden will ich mich mit dem praktisch oft noch wichtigeren Fall befassen, das der Moment- 
vektor eine Funktion des Drehvektors ist, etwa ein Moment, das den Kreisel um eine Hauptachse 
von der Ruhe oder von einer schon vorhandenen Drehgeschwindigkeit aus bis zu einer vorge- 
schriebenen Drehgeschwindigkeit oder auch iiber alle Grenzen schnell antreiben soll. 

Zumeist ist aber gar nicht dieses Antriebsproblem (das ja véllig elementar ist) das Wichtigste, 
sondern die Frage, wie eine solche, zunachst stérungsfrei gedachte Antriebsbewegung durch zu- 
satzliche Stérungen (wie etwa kleinere oder gréBere seitliche DrehstéBe) geandert wird, Wir werden 
namlich sehen, daB (im Gegensatz zum symmetrischen Kreisel) derartige Stérungen die eigentliche 
Antriebsbewegung ganz entscheidend beeinflussen kénnen, so daf die erstrebte Drehgeschwindig- 
keit nicht und nicht einmal annahernd erreicht wird, ja sogar daB die Bewegung sich umkehrt und 
in der Drehgeschwindigkeit periodisch wird, oder da die Drehgeschwindigkeit nach einer Bewe- 
gungsumkehr iiber alle Grenzen wachst. Die Mannigfaltigkeit der zum Teil paradoxen Erschei- 
nungen ist gréBer, als man wohl erwarten kénnte. 3S 

DaB unsere Fragestellung fiir manche Probleme der Raketentechnik bedeutsam ist, brauche 

_ ich kaum zu betonen. 

Die hier zustandigen Eulerschen Gleichungen in dem kérperfesten, rechtshindigen Haupt- 
achsenkreuz (x, y, z) durch den Schwerpunkt haben mit den Hauptdrehmassen A, B, C, den Kom- 
ponenten w,, w,, @, und M,, M,, M, des Drehvektors und des Momentvektors der Selbsterregung 


y2 
und also mit den Tragheitszahlen 
B—C A—C A— B 
Gs Vir aR gm Epa 2 ee CG (1) 
sowie den Momentfunktionen oe 
: _M. 2 
egy Meme ee) 
die Form ; 
mH Wy Oy = M, » (3) 
dw 
a a Os De = My > (4) 
do} 
: Fe Hg Wy Wy = M,- (5) 


ll 
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Schreiben wir noch die Rangordnung 
ASB G (6) 


vor, schlieBen also den symmetrischen Kreisel aus — er wird sich nahher durch eine kurze Bemer- 
kung erledigen lassen —, so sind die drei Tragheitszahlen wesentlich positive GréBen: 


“; 0 (i=1,2, 3). 


Im allgemeinsten Falle waren nun also m,,m, und m, vorgeschriebene Funktionen von @,, Wy, 0;- 
Es hatte aber wohl wenig Reiz, diese Funktionen beliebig zu wahlen (und dann die Lésung des 
Problems etwa elektronisch durchzufiihren); wir wollen hier lieber fiir verniinftig vorgeschriebene 
Funktionen, wie sie sich auch technisch leicht verwirklichen lassen, auBer der zahlenmaBigen Lo- 
sung des Problems allgemeingiiltige Eigenschaften der Bewegung finden. 

Unsere erste Einschrankung fiir die Funktionen m,, m,, m, besteht darin, da wir nur Selbst- 
erregungs-Momente betrachten, deren Vektor in einer der drei Hauptachsen liegt; unsere zweite 
Einschrankung besagt, daB ein solches Moment nur eine Funktion der Drehkomponente jener 
Hauptachse sein soll. 

Wir beginnen mit der x-Achse, setzen also 


m, =™m,(@,), my=0, m=0. (7) 


Obwohl es nicht schwer ware, die Lésung dann fiir ziemlich allgemeine Funktionen m,(w,) zu 
erledigen, so wollen wir — dies ist die dritte Einschrankung — die Momentfunktion m, mit einer 
positiven Konstanten 4 in der Form 


Mm, = p(w — wz) (8) 


ansetzen, weil sie in dieser Gestalt nicht nur eine sehr durchsichtige explizite Lésung des Problems 
erlaubt, sondern auch technisch die haufigste Momentart scin mag, die man verwirklicht, um eine 
gewiinschte positive Drehgeschwindigkeit w, um die x-Achse zu erzielen. Denn die zu wz propor- 
tionale Fliehkraft kann in bekannter Weise zur Regulierung eines solchen Momentes dienen; und 
man erkennt leicht — wir werden dies natiirlich alsbald bestatigen —, daB ein solches Moment 
mindestens im stérungsfreien Falle das Bestreben hat, jede Anfangsgeschwindigkeit w,, > — Wo 
gegen den positiven Wert j, hin zu steuern, eine Anfangsgeschwindigkeit w, , << — Wy aber negativ 
iiber alle Grenzen zu treiben. 

Die zusammengehorenden Verabredungen u > 0 und w, > 0 sind dabei nicht irgendwie ein- 
schrankend, weil man die positive Richtung der x-Achse stets geeignet wahlen kann. 


2. Die Grundbewegung. Wir erledigen zuerst den zwar ganz elementaren, aber als Grundlage 
fiir unsere spatere eigentliche Untersuchung wichtigen Fall, daB auBer der Vorschrift (8) fir m, 
auch noch w, = 0 und w, = 0 sein soll, was nach (4) und (5) wegen (7) mdglich ist, also den stérungs- 
freien Antrieb. Er gehorcht zufolge (3) und (8) der Differentialgleichung 


dax 
ap = (Op — 5) - (9) 
Die Lésung lautet je nach dem Anfangswert w,, zur Zeit t = 0, wie man leicht nachpriift, 
I) far, , > @,: Wz = Wy Cig f Wy t + Ar Ctg ay | 
0 
2) Pur 2 = We: O, = Wy» 
3) fir — 0) <.0;,< Wy: O, =o, XG f Wot + Wr Ty), (10) 
0 

4) fir w,, = — Wp: Wx =— Wy» 


5) fiir @,)< — @p: On = Wy Stg (uct + We tg 2) 


und wird dargestellt durch Kurven, die in Abb. 1 die entsprechenden Nummern tragen. Die Zeit t, 
der gré8ten Winkelbeschleunigung in der Kurve 3” und die Grenzzeit t,, In der Kurve 5 sind 


———-O) 


aes (11) 


pe OAs eee 1 
PO We 8 Wy boo = Tg, Ut Ctg 


Wo 


Wir wollen diese stérungsfreie Antriebsbewegung die Grundbewegung nennen. 


yas 
1h 


REO Fi 
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Ks ist ur di dj 
Sarees pope ae Patethonns der durch zusatzliche seitliche DrehstoBe gestérten 
eutungsvoll, die Lésung fiir die Grundb in ei i 
4 s fir die | ewegung in einer zweiten, anderen 
estalt zu wiederholen, die zudem alle fiinf Falle in einer Gleichungsform ae cients it indem 
man statt der Zeit t den von t = 0 aus gemessenen Drehwinkel 


t 
Ps = fo, dt (12) 


als unabhangige Variable einfiihrt (die nachher den Schliissel zur Lésung des Stérungsproblems 


bilden wird). Multiplizi : = é 
Be ibeae tipliziert man (9) mit @, = dy,/dt, so kommt als neue Differentialgleichung der 


d(x) 
dpx 


toe We == 2 we: (13) 


Fr0 
Abb. 1, Abb, 2. 


und diese hat die schon an die Anfangsbedingung w, = w,, fiir t = 0, also fir y, = 0 angepabte 
Lésung 


wo, = V U9.) (14) 


Up(p.) = 02 + (w2y — w2) 7", (15) 


wobei iiber das jeweilige Vorzeichen von /U,(%.) in (14) spater noch Festsetzungen zu treffen 


sein werden. 
Wir merken nur noch an, daB aus (14) durch eine Quadratur 


mit 


Px 
dx 
i= —-— = t(9,) 16 
J itor oy 
und dann durch Funktionsumkehr 


folgt, womit w, gema® (14) wieder als Funktion von t dargestellt ist. Wiirde man die Integration 
(16) und die Funktionsumkehr (17) ausfiihren, so kame man natiirlich gerade wieder auf die Lé- 
sungen (10) der Grundbewegung. 

Die Deutung der Lésung (14), (15) ist in Abb. 2 angegeben, wo die Bewegungen 1, 2, 3, 3’, 3", 4 
und 5 von Abb. 1 mit den gleichen Nummern als Kurven dargestellt sind. Dabei ist 2u beachten, 
daB dp, gem&B (12) stets das gleiche Vorzeichen wie w, hat, ferner dab y, = 0 zu t= 0 gehért, 
weiter daB gemaB (14) nur positive Werte der Ordinate Uj zu einer reellen Bewegung gehéren, und 
endlich daB das Vorzeichen von VU, wegen (14) stets das gleiche ist wie das von w,. 

: 11* 
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So sind die Kurven 1, 2, 3 und 3’ ohne weiteres verstandlich. Bei der Kurve 3” ist anfanglich 
W, <0, also auch //U, negativ bis zu dem Wert @,,, fiir welchen Uy = 0, also die Drehung w, 0 
wird, namlich gemaB (15) 
ee ath 2b (18) 
eas fr haat we 
Dann kehrt die Drehrichtung von w, um, da nunmehr /U, und dy, positiv werden; der rechts- 
gerichtete Pfeil deutet hier wie schon bei den Kurven 1, 2, 3 und 3’ an, daB der Drehwert w asympt- 
otisch angestrebt wird. Bei den Kurven 4 und 5 ist dauernd /U, < 0 und dy, < 0, was negativ 
konstante bzw. negativ iiber alle Grenzen wachsende Drehgeschwindigkeit w, bedeutet. Alles 
dies ist natiirlich in voller Ubereinstimmung mit den in Abb. 1 dargestellten Bewegungen. 


3. Die gestérte Bewegung. Wahrend fiir einen symmetrischen Kreisel (B = C, also a, = 0) 
diese Grundbewegung w, gar nicht von etwa hinzukommenden seitlichen Drehungen w, und @, 
beeinfluB&t wird, weil mit «, = 0 die Gleichung (3) nicht mehr mit den Gleichungen (4) und (5) ver- 
koppelt ist, so kann die Grundbewegung beim unsymmetrischen Kreisel durch seitliche DrehstéBe 
ganz wesentlich gestért werden, wie wir jetzt zeigen wollen. 

Der bisherige Nullpunkt t = 0 der Zeit soll also der Beginn einer Storbewegung w,, w, sein, 
die den Gleichungen (4) und (5) gehorcht und mit den Anfangswerten w,, und w,, einsetzt, von denen 
mindestens einer von Null verschieden sein soll. 

Mit der unabhangigen Variabeln ¢, (12) anstelle der Zeit t gelingt es nun, das in w,, w,, @, nicht- 
lineare System (3) bis (5) in ein System iiberzufiihren, das in wz sowie w, und q, linear ist. Man 
erkennt das sofort an den beiden Gleichungen (4) und (5), die, wie man fiir m,= 0 und m,= 0 
nachprift, ausgedriickt im Drehwinkel ¢, die Lésungen haben 


Oy =, VX» cos (og %3 Pz ot 0) ; O, = WM, xg sin (/og 3 Px a a) (19) 


mit den Konstanten @, und g, die durch 


2 2 / 
2_ Fo Orzo a W 
Ce oh OW (20) 
Xe as a3 Hyg 


definiert sind, wobei wir also jetzt w, == 0 voraussetzen und die Quadratwurzeln in (19) und (20) 
sowie @, positiv wahlen diirfen. 


Mit (19) und (8) geht die Gleichung (3) iiber in 


d(w3 — , = 
i +2uoi =2uoarz+ wf? oy Vo. O& sin 2 (|/%» 3 Dy + 0)» (21) 


also in der Tat in eine nun in w; lineare Differentialgleichung. Ihre Lésung lautet 


o, = VUQ,) ; | (22) 
dabei ist 


U(px) = Ulex) — Ui(pz) (23) 
mit der friiheren Grundfunktion U,(p,) (15) und mit der Stérungsfunktion 


I 1 Xs Xs ee = 
Us(p,) = 4 of 2122 [eos 2 (Vay 05¢2 +e +5) —cos2(gto)e2"%=], (24) 


und hierin sind (neben den fritheren Konstanten w, und @) die beiden weiteren Konstanten x und o 
definiert durch 


Vos as 


mit positiven Vorzeichen der Quadratwurzeln. 


Die Stérung driickt sich also darin aus, daB von der Grundfunktion U, dic Sta 4 
tion U,(p,) abzuziehen ist. o(Pz) die Stérungsfunk 


Wir bemerken noch, daB fiir Y, = 9, alsot = 0 
U(0) =o, U,(0) =0 (26) 
wird, so daB die Lésung (22) die Anfangsbedingung erfillt. 
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i a ooo der Bewegungen @,, die durch (22) bis (24) dargestellt sind, insbesondere 
ie Mannig altigkeit der Abweichungen von den zugehérigen Grundbewegungen (14) und (15) 
wie a infolge der durch (20) ausgedriickten Stérungen auftreten kénnen, ist erstaunlich prob, 
oo ooo unserer fritheren Abb. 2 der Grundbewegungen mit den Abb. 3 bis 6 der oS: 
storten Dewegungen zeigt, die wir jetzt erklaren wollen (und di ine A i a 
risch berechneten Kurvenblattern wiedergeben). ioe 


# 
dy 


Abb. 3. & 


PY 
% 
a ~ 
G 
/ * 
5 


Abb. 4. 


Abb. 3 und 4 gelten fiir cos 2 (9 + 0) > 0, Abb. 5 und 6 fiir cos 2 (9 + 0) < 0, und zwar Abb. 3 
und 5 fiir eine kleine, Abb. 4 und 6 fiir eine groBe Stérung (je verglichen mit der angestrebten Dreh- 
_ geschwindigkeit «,). Die Nummern der Kurven in Abb. 3 bis 6 entsprechen den Nummern in 
Abb. 2. Die Ordinatendifferenz der Kurven U,(p,) und U,(p,) ist nach (22) und (23) ein Maf fiir 
w2, und die Kurven U,(y,) sind jeweils nur soweit eingezeichnet, als U > 0 bleibt; und auch hier 
ist wieder zu beachten, daB dy, und VU stets das gleiche Vorzeichen wie w, haben. 

Im Einzelnen bemerken wir zu Abb. 3: Bei den Kurven I, 2, 3, 4 und 5 unterscheiden sich die 
gestérten Bewegungen von der Grundbewegung nur dadurch, daf noch eine (von q, abhingige) 
Schwankung in w, hinzukommt, so daB also der angestrebte Wert «, auch nicht mehr asymptotisch 
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genau erreichbar ist, sondern schlieBlich mehr und mehr genau umtanzt wird — man kénnte das 
eine periodisch-asymptotische Annaherung nennen. 
Bei hinreichend kleinen Anfangswerten |w,,| tritt das ein, was fiir w,, = 0 in Kurve 3’ darge- 
stellt ist, namlich ein periodisches Hin- und Herschwanken von o, zwischen einem positiven und 
einem ebenso groBen negativen Wert, aber keinerlei Annaherung mehr an den angestrebten Wert Wp. 


44 


Abb. 6. 


Die Bewegung 3” zeigt das gleiche Verhalten wie die Grundbewegung 3” in Abb. 2, nur daB 
sie sich nach einer Richtungsumkehr dem Wert w, bloB periodisch-asymptotisch nahert 
Ganz neuartig und aus dem Rahmen der Grundbewegung vollig herausfallend sind die Bewe- 
gungen 3* und 3** in Abb. 3. Fangt die Bewegung 3* mit einem positiven Wert w.. an, so nahert 
sich die Kurve dem Punkt A, und zwar, wie man zeigen kann, erst mit t > 0d. Dort ist die Dreh- 
komponente w, = 0, was aber keineswegs Ruhe bedeutet, da die pulsierenden Drehungen w, und 
@, nach (19) vorhanden sind, die einen in der (y, z)-Ebene umlaufenden Drehvektor hilden De 
Zustand im Punkte A ist, wie man ebenfalls zeigen kann, labil in dem Sinne, daB w jetzt von Null 
an in positiver oder negativer Richtung zu wachsen vermag, allerdings theoretisch auch erst nach 
unendlich langer Zeit, praktisch aber bei einem zufalligen kleinen AnstoB schon nach endlicher 
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Zeit. Das bedeutet Fortschreiten auf der Kurve 3* entweder nach rechts mit periodisch-asym- 
ptotischer Annaherung an ,, oder nach links mit Umkehr im Punkte B und Wiederholung des 
ganzen Spiels. Ahnliches gilt fiir die Kurve 3**. ‘ 

Nach dieser ausfiihrlichen Erklarung von Abb. 3 werden Abb. 4 bis 6 von selbst verstindlich 
sein. AuBer den schon bei Abb. 3 besprochenen Effekten kann infolge einer starken Stérung noch 
der Fall eintreten, daB bei keiner Anfangsgeschwindigkeit @,9 eine wenigstens periodisch-asym- 
ptotische Annaherung an den erstrebten Wert «, méglich ist, sondern daB (auBer periodischen Be- 
wegungen) stets nach einer Bewegungsumkehr die Drehgeschwindigkeit «, negativ tiber alle Grenzen 
wachst; dies zeigen in Abb. 4 die Kurven 1, 1*, 2 und 3’’, und in Abb. 6 die Kurven I] und 1*. 

Ubrigens sind die Doppelstriche, wie wir sie in Abb. 2 bei Kurve 3” angebracht hatten, weiter- 
hin tiberall (auBer bei Kurve 3” in Abb. 3) weggelassen, um die Ubersichtlichkeit nicht zu gefahrden, 
und auch weil sie an manchen Stellen offensichtlich sogar mehrmals oder beliebig oft hatten wieder- 
holt werden miissen. In Abb. 5 kénnte die Kurve 5* natiirlich nach rechts fortgesetzt werden 
(wie eine Kurve 1), und in Abb. 6 kénnte man die Kurve I’ auch 5* bennnnen. Man sieht 
leicht, daB es noch viele weitere Kurventypen gibt, als die eingezeichneten, z. B. in Abb..6 
eine Kurve 5**, deren linker Teil vom Typ der Kurve 5* in Abb. 5 ware. 

Wir erwahnen noch, da8 man w,, Dy: Oz ,explizit als Funktionen von t ganz nach dem Schema 
von (16) und (17) darzustellen hatte, nur natiirlich mit dem Unterschied gegen Ziff. 2, daB die bei 
der Quadratur (16) und bei der Funktionsumkehr (17) auftretenden Funktionen bisher nicht tabu- 
fiert sind und wegen der groBen Anzahl der Parameter auch kaum allgemein tabuliert werden 
kénnten. 


4, SchluBbemerkungen. Der Vollstandigkeit halber merken wir noch an, daB der Fall 
m,=0, m=0, m, =p(o5— 2) (27) 


eines Antriebsmomentes in der z-Achse sich durch die gleiche Rechnung erledigen la8t und zu gleich- 
artigen Ergebnissen fiihrt; man braucht nur die entsprechende Vertauschung der Indices x und z 
sowie 1 und 3 vorzunehmen. 


Aber auch der Fall des Antriebs um die y-Achse 


m,=0, m,="Uoj—o}), m,=0 (28) 


2 


gehorcht einem ebenso gebauten Formelsystem, jedoch mit einem wesentlichen Unterschied. Man 
mu jetzt x mit y sowie x, mit — o, vertauschen, und dann gehen die Kreisfunktionen in Hyperbel- 
funktionen iiber, und es sind dabei die zwei Unterfialle w,,/x, < w,,/x; in der gleichen Weise zu 
unterscheiden, wie ich das an einem anderen Problem friiher getan habe1, wobei der Grenzfall 
W,9/%, = ,9/%3 durch einen einfachen Grenziibergang folgt. 

Die weitere Diskussion ergibt, daB das Moment m, die erstrebte Drehgeschwindigkeit «, bei der 
kleinsten Stérung stabil erzielen kann. Dies steht ganz im Einklang mit der bekannten Tatsache, 
daB die Drehung um die mittlere Hauptachse y selbst ohne Moment stets labil ist, und daher 
wird dem Antrieb (28) kaum je eine praktische Bedeutung zukommen (weswegen wir darauf ver- 
zichten, sein Formelsystem explizit anzugeben). 

Es soll noch darauf hingewiesen werden, daB die hier nicht beachteten Reibungseinfliisse unsere 
Ergebnisse teilweise einzuschranken vermégen, weil sie die Schwankung (19) des Drehvektors der 
Stérung mehr oder weniger rasch abklingen lassen, so dafs dann die Bewegungen der starken Sté- 
rung (Abb. 4 und 6) in die der schwachen Stérung (Abb. 3 und 5) iibergehen oder sich sogar genahert 
oder praktisch genau in die Grundbweegung zuriickverwandeln kénnen. 

In vielen Fallen ist mit der Kenntnis der GréBe und Lage des Drehvektors (w,, @,, @,) im 
Kérper das technische Problem schon gelést. Will man auch die Bewegung des Kérpers im Raum 
als Funktion der Zeit kennen, so hat man aus den Komponenten @,, w,, , des Drehvektors etwa 
die drei Eulerschen Stellungswinkel des Kérpers zu bestimmen. Die Methode zur Lésung dieser 
einfachen kinematischen Aufgabe ist aber wohlbekannt? und mu8 nur noch auf die jetzigen Funk- 
tionen ©,(t), ,(t) und w,(t) angewendet werden, was allerdings recht umstandliche Rechnungen 
mit sich bringt. 

(Eingegangen am 17. Juli 1959) 
Anschrift des Verfassers: Professor Dr. R. Grammel, Stuttgart N, Robert-Bosch-StraBe 101. 


1 Vel. S. 85 der in FuBnote ! von S. 153 zitierten zweiten Untersuchung. ; 
2 Sithe etwa: R. Grammel, Der Kreisel, Bd. 1, § 9, Ziff. 6, 2. Aufl. Berlin/Géttingen/Heidelberg 1950. 
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Zur Theorie der Unterwassertragfliigel bei wellenformiger Anstromung 


Von W. H. Isay 


1. Einleitung. In meiner vorangehenden Arbeit! habe ich das stationare Problem von nahe 
der Wasseroberflache fahrenden Tragfliigeln (homogene Anstrémung) behandelt. Vielfach be- 
finden sich Unterwassertragfliigel jedoch in einer wellenférmigen Anstrémung, so daf} sich der am 
Profil wirksame Anstrémwinkel dauernd andert, und somit ein instationares Problem vorliegt. 

Die ersten Untersuchungen iiber instationére Unterwassertragfliigel wurden von Kaplan? 
durchgefiihrt. Jedoch wird in diesen Arbeiten (wie ja verschiedentlich auch in der stationaren 
Theorie, vgl. Einleitung Arbeit I) die Strémungsrandbedingung am Fliigelprofil nicht erfiillt, und 
daher bleiben auch die Wirbelstarke (Zirkulation) und die GréBe der Fligelkrafte unbestimmt 
und werden durch Ausdriicke aus der Hydrodynamik der Fliigel im unbegrenzten Medium er- 

anzt. 
: In einer kiirzlich veréffentlichten Arbeit? wird unter Vernachlassigung der Schwerkraft [g — 0, 
d.h. in dem Sonderfall sehr hoher Froudescher Zahlen] die Randbedingungsintegralgleichung 
aufgestellt und gelést. 

In der vorliegenden Arbeit wird eine Theorie des instationaren Unterwassertragfliigels angege- 
ben, die auch die Auflésung der aus der Strémungsrandbedingung am Fliigelprofil entstehenden 
Integralgleichung fiir die Fligelzirkulation einschlieBt. Ferner werden die Fliigelkrafte und die 
Form der freien Wasseroberflache in der Nahe des Tragfliigels berechnet. Wie in Arbeit I wird 
auch hier die Strémung als eben vorausgesetzt (unendliche Spannweite). Die numerische Auswertung 
dieser Theorie durch Zahlenbeispiele soll einer weiteren Arbeit vorbehalten bleiben. 


2. Das Geschwindigkeitspotential eines Wirbels unter der Wasseroberflache bei wellenformiger 
Anstrémung. Wir betrachten eine Tragflache in der mittleren Tiefe h unter der bei y ~ 0 liegenden 
Wasseroberflache. Die Tragflache liege in einer Anstrémung, die sich aus einem homogenen 
Hauptanteil uy und den noch genauer zu definierenden wellenférmigen Anteilen u,(x, y,t) und 
v,(x, y, t) zusammensetzt, wo u und v die x- und y-Komponente der Geschwindigkeit bedeuten. 
Dabei sei u? << u2, vi << u?. Die Tragflache werde wie iiblich durch eine gebundene Wirbeldichte 
y(é, t) (—a S& <a) ersetzt. Da es sich um ein 
instationdres Problem handelt, werden durch die 
zeitliche Anderung der gebundenen Wirbel freie 
Wirbel induziert, die mit der Strémung hinter der 
Tragflache abflieBen. Um spater die fiir dieses Pro- 
blem benétigten Ausdriicke bequem durch Superposi- 
tion bilden zu kénnen, wollen wir zunachst das Ge- 
schwindigkeitspotential ©, eines an der Stelle (€, —h) 
gelegenen Wirbels a y(é, t) ableiten. 

An der Wasseroberflache ist der Druck iiberall konstant. Aus der Bernoullischen Gleichung er- 
halt man dann unter der Voraussetzung, dab das Quadrat aller durch den Tragfliigel auftretenden 
Stérgeschwindigkeiten klein ist gegen u2 und mit ®, als Geschwindigkeitspotential des Wellen- 
anteils der Anstrémung 


0 a 
a (Pit D,) + uy = (®, + ®,) +gY=c). 


1 W. H. Isay, Ing.-Arch. 27 (1959) S. 295 (Arbeit I). 
* P. Kaplan, A hydrodynamic theory for the forces on hydrofoils in unsteady motion, Dr. Sc. Thesis, 


Stevens Inst. of Technology Hoboken, New Jersey, 1955. (Herrn Prof. Dr.-Ins. G. Weinbl 
danke ich fiir die freundliche Uberlassung dieser Arbeit), ‘ 2 iin one ae eee eee 


P. Kaplan, The drag of hydrofoils in unsteady motion, Stevens Inst. of Techn, ETT Report 617. 


P. Kaplan, The forces acting on hydrofoils in unsteady motion; Proc, 9th Int. Congress of Applied Mechanics 
Briissel 1956. : 


P. Kaplan, The Waves generated by the forward motion of oscillat 1 istributi “a 
western Conference on Fluid Mechanics 1957. Se ee 


P. Kaplan, The forces and moments acting on a tandem hydrofoil i 
ETT Rese ee. § on a tandem hydrofoil system in waves, Stevens Inst. of Techn. 


3 Wen-Hwa-Chu und Abramson, J. ship research 3 (1959) S. 50. 
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Dabei gibt Y(x, t) die genaue Form der Wasseroberflache an. Weit vor dem Fliigel (fiir x > — 00) 
ubt dieser keinen Einflu8 mehr aus; es bleibt dort also 


0D. oD 
Bye) =e) Sn. 


Da bei einer wellenférmigen Anstrémung mit der Wellenlange 1 fiir jeden Zeitpunkt t 
x+l 
lim f Y(x,t)dx =0 


x—>— OC xX 


sein muB, folgt c(t) = 0. Damit erhalten wir fiir die Form der freien Wasseroberflache 
Ug @ 
Y(mt)=—4 FO, +6) —" 2 0,44). (1) 


Die kinematische Strémungsrandbedingung an der freien. Wasseroberflache lautet 


dV AE OY. OYe neo 
de et M0 Ge = By (i +O), 


und wegen (1) wird daraus 

C) a2 oe? 2 
8 x (M1 + ®,) setG aes, (Onc ey) + 2 => (M1 + ®,) +55 (GP, + D,) =0, (y=9). (2) 
Da die Anstrémungswellenbewegung unabhangig von dem in der Flissigkeit befindlichen Trag- 


fliigel der Relation (2) geniigen muB, gilt auch fir @, allein die Bedingung 


a®@,, 2 &D, PD, , PD, eid ; 
cay to Ox? stot Uy aa ag (y =0). (2’) 


& 


Die Wellenanstrémung habe die Kreisfrequenz w; wir setzen sie wie iblich als ebene Welle an 
mit dem Potential 


@, = Be’ ei(ot—Ax) 4 B ety e—ot—Az) 


Die Bedingung (2), (2’) liefert dann zwei mégliche Wellenzahlen /, namlich 


ay = Fal» Ja = a Ma (4 > 9, M2 >), 


wobei 
(SS | Cos eee Rie Ns ie eee Se (3) 
b= Sear or coy by = aes amen aa x wag 
ist. Damit erhalten wir die Wellenanstrémung unseres Problemes in der Form 
Se-inn . . 4 utiann 
D = B, eo es + B, ei?! e%o 
ES iy—inus — | Siytixr (4) 
+ Be ?tewo + B,e—'?* e%o _ 
oD aD 
My = > tis yy (ui< uw » vi<up) « | 


Wir kommen nun zu der Berechnung des Geschwindigkeitspotentials ©, eines an der Stelle 
(&, — h) liegenden Wirbels der Starke a y(€, t). Wir setzen an: 


VE, t) = d(€) +6,(8) ef ?* + (6) Fe", (5) 


@, = 5" Hawi? sin x —£) FS” eg sin Aw —£) 2 4 OF + OY. (6) 
0 0 


Die beiden ersten Glieder in (6) stellen das auch in der stationaren Theorie auftretende Potential 
eines Wirbels in unendlich ausgedehnter Flissigkeit und des an der Ebene y = 0 gespiegelten 
Wirbels dar. Das Zusatzpotential ®* ist so zu bestimmen, daf es zusammen mit den beiden ersten 
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Gliedern der Relation (2’) geniigt. Dieses fiihrt nach elementarer Rechnung auf die Bedingungs- 


gleichung 
oc 
ody, 2 PDF Dy , POY g Bees 
By ute remietee rr mare = ay(&, t) = ine sinA(x — &)d/. (7) 


0 


Der Zusatzanteil ®** mu dann fiir sich der Gleichung (2’) geniigen, und er ist weiter so zu bestim- 


men, daB 


Seed pre shoes >) Wee | 
fim (2! =0, lim 57 =0, lim ar eae (8) 


t= —— Co) Ox x—>— CO > — CO 
ist; denn der Wirbel iibt ja weit voraus keinen Stérungseinflu8 mehr aus. 
Der Ansatz 
Dy = eft f eby (04-9 GA, €) +e ME—9 GA, €)) dh 
+ e—iat f ely (e— te —4) GA, E) + ef —4) G,(, €)) dA + f ey G,(A, €) sin A(x — €) dd 
f 


liefert gem4B Gleichung (7) nach elementarer Rechnung 


(uz 72 + 2uywA—gd + w*) G, = axultjem*, 
(uz 22 —2 uy oA—gA + 0%) 6, =— 5 0g 5,6) e**, 
(uj 48 — g A) Gy = ——F 8,(6) e-**. 


Wir setzen noch 2 = s und erhalten nach (5), (6) das gesuchte Wirbelpotential 


oe = (So(€) + 6,(&) ef °* + 6,(é) es) fleets 26 —9)) sin A(x — £) 7 
3 J 
[ese Bore tas sin Beeb) 
ee ae we—pt+lpat @ ae a a ee 
Leis Seen pats eae pr ay ia eb 
Se ww 2h Q+ 2 ae ay ery ER, |e (9) 
— dango foe? sae du \  —-Se-» g 
7% % Te =?) ig ay “0 ee ae 
0 
es 
9 1 ears e aie 0 e ; 


SS reas 
1 (, ng ho) fae 8 


Higa) tines E28) 
up e ub 
fy — [lp 


—+6,(6) e—iot 


J 


Um die Entwicklung der Theorie und ihre Ubersichtlichkeit nicht durch Zwischenbetrachtungen 
a Sey ce haben bale os E ormel (9) gleich auch das Ergebnis fiir ®** (das sind die Glieder 
ne Integral) mit angegeben. ®O** ceniigt der Relati aS i a 
Sogn eames 2 0 coon er Relation (2); denn y,, “, sind gemaB (3) Wurzeln 
Damit erfillt auch ®, die Randbedin ‘ys ii i 
ge 5 gung (2’); ®, geniigt ferner der Bedingung (8) und hat 
et alle geforderten Eigenschaften. Um das leprtere zu erkennen, bendtigt a ae wichtige 
ntegralformeln, die auch in der ganzen weiteren Theorie eine entscheidende Rolle spielen, so daB 


* 
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wir fiir sie neue Abkiirzungen einfiihren wollen. Es gilt 
co 
—He by) —inS@—-) dy 
L(x, y, §, fy) = : : v ; why 


x [Ci(us go) +iSi(14 i («—8) Lit e 


, (10) 
—4 5 b—y) —im4 w—8 Me j 
=—e e “9 [ost72-+ Inti I [(x—&)? + (h—y) ]-+i aretg -—* 
0 
co il Ps ; : 
ace) aie +a—yrl, (+ firx>&, —firx<&), 
—1 = (h—y) —in,-% («8 
lim T,(«, y, &, yh) = Fine “0 e af) 4 
X%¥—>-+ © 
g g g 
BP a a ee 8) din Q- (h—y) 
I(x, y,&)= | e 0 e “0 Begarce. 0 x 
: = hy) 
40 
19S (e-8 +122 @—9 a 
CO) © Ce s,(2 f =—£) |= 2° = 
x4 ¢ ( me ))+i (x ) ri #4 [Fe e—af (11) 
0 ah 
L(y apie 
ae oe oa8 | : josst72 + In SE 4 3 [(a— €)? + (h—y)?] —i aretg ;—— 
Vets a Gx—ie—ntynt, (—firx >€, + fire <é), 
fn (x, v5 é). = 0. 
x—>-+ oo 
—u 5 b—y) in 5 @—§ a 
I(x, y, 6) = e 0 e 0 pratip 
; 0 
(a 4+ ip)(h—y—ix+ié 
ent + if)(h—y + (erty = In & (0? + B) -++In [(x — é)? + (h— y)?] (12) 
— iaretg = — priate lt Sa =f = (a +1 pGx—ig—hty)]|, 
lim I3( 98). = 0, (e207 Pp 9), 
(edrees 2 
hg? ee du 
L(x, y, €) = e pta—ifp 
0 
go eos “lo 5772 +7 & (t+ ae In [(x—€)2 + (h—y)?] (13) 
— i aretg 7 — Set ae nl {ua 5 («x —if) (iz—ig—h ty) . 
n=l 
0). 
lim I,(x,y,é) =0, = (a=0, B>9) J 


x—>-+ 00 
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Fiir den Beweis der Integralformeln (10), (11), (12), (13) wird auf Ziff. 10 verwiesen. Die GréBen 
« und f sind in folgender Weise definiert : ; 


mM. | 
1; 1 
0 Eee ge ee 


Wir beschranken uns dabei auf den Bereich Q > 1/4; dieser ist fiir die bei Tragfliigelbooten vor- 
liegenden hohen Froudeschen Zahlen allein von Interesse. 

Mit Hilfe der Integralformeln (10), (12), (13) ist es nicht schwer, zu zeigen, dafs das Wirbel- 
potential (9) auch der Bedingung (8) geniigt. Die zeitunabhangigen Glieder in (9) stimmen natiirlich 
genau mit den bekannten Ausdriicken aus der stationaren Theorie iiberein (vgl. Arbeit I, Ziff. 2). 


3. Das Geschwindigkeitsfeld der gebundenen Wirbel. Mit dem Potential @,, eines an der Stelle 
(£, —h) liegenden Wirbels aus Formel (9) kann nun leicht das Geschwindigkeitsfeld u,, v, der 
gebundenen Wirbel des Unterwassertragfliigels erhalten werden, und zwar ist 


(14) 


a@,, dé _ f a, dé 


u,= | —-— a) 


v ax a” Y dy a’ 


Wir iibergehen die elementare Zwischenrechnung und geben gleich das Ergebnis an: 


teak ; iwt S (£\ p—ioat yth y—h 
ge eas J LONE) oi) Se le—s Tot t+ eH tO—m| 
Lg ; a aaa eae g —% (h—y) . nog 
ini dae ua | 6o(§) [Ls Wes 1) ae I(x, oe) & 1)] dé =F uz ie re [ 0 Sane (x — &) dé 


1 Soar - i 1 aes - a 1 
ee gE. ve | 6; (£) lz(1—#4) I,(x, y, &) iat i al Tale 9, 8) + = S iz, libs Yo Et) 


+ lt heb u)| d&g Set [5G [21 +44) hex 


Pass 
1 ats Vea a ea i ——— yb Seas ee 
=p ce (1 a I(x, y, &) ae E(%29; €, 4) BEE TG y.6: #)| dé 
: a po aes g A 8 
+ 8 yias MR Dees ie ca) was So or med Pe tat 1) 
ae : uz ‘ I oe Y e 0 oh eee 0 e 0 dé 
— . a j= agi a & g : zg 
Peay ape ea ee ee ee eg eg OT) et) 
2s cee i (6) fae pe tg Ue Ide, (15) 
1 P eeu 
nfm yt) =—yz | Ld &) +6,(é) et +6 Sap Ane x—eEé ee 
te 3) o( ) 1 ) Sig 1(€) e€ ] Ge) Liye kt @— Se dé 


is, eu Ee Soe nent 
Soe [ae [L(x ¥, é; 18 ine L(x, y, &; 1)] dé — Fe “o ” f exe ws £ (=e) dé 


bie r 1 5G: 1 : 
De uz © fae [z(1-+5) T(x, y, &) +3(1 +15) I(x, y, &) ectnisre I(x, y, §, uy) 


1 SCE pert M&W ee oN Rosana pe 
awe (1—i 5) I,(x, y, Speier L(x, 9,6, a sare Ty.¥.5,49)| dé 


XXIX. Band 1960 W.H. Isay: Zur Theorie der Unterwassertragfliigel bei wellenformiger Anstrémung 165 


a g ais g ens < 
; Hig (h—y) —i—— («—£) ny —ts— (h—y) —i pn, (*«—£) 
s Looe ge 2 lu 2 2 
—p et 16 | —e€ 0 e “o ag 40 e “o d 
2 u2 : ls) My— Pg a : 
mons (F 


& * & c 
My ita EY ae 


e & eed, 5 
g reese oo | Me —Ha 7) (h—y) tpg 2 (x | 
— ~e 6,(é e e + e 0 e “o dé. (16 
2 us | a( ) Mi— He He— fy ¢ ( ) 


4. Das Geschwindigkeitsfeld der freien Wirbel. Durch die zeitliche Anderung der gebundenen 
Wirbel werden freie Wirbel induziert, die mit der Anstrémung uy hinter dem Unterwassertragfliigel 
abflieBen. Der kleine wellenformige Anteil der Anstrémung und die gegenseitige Beeinflussung der 
freien Wirbel wird dabei vernachlassigt. Ein an der Stelle (&, — h) liegender gebundener Wirbel 
der Starke ay(&,t) induziert eine kontinuierliche Folge von freien Wirbeln der Starke 


ae Re ee 8) dee 
Zea ret— = ‘a= i 4,6 Te = eae Fe 
0 


Ug 
(ES X<oo). 

Wir wollen dabei annehmen, daf alle diese freien Wirbel den gleichen Abstand h von der Wasser- 
oberflache behalten; dieses ist natiirlich nur in der Nahe der Tragfliche naherungsweise richtig. 
Da jedoch die in der Nahe der Tragflache befindlichen freien Wirbel den gréBten Einflu8 ausiiben, 
erscheint es zulassig, die obige Annahme zu machen. Um das Geschwindigkeitsfeld bzw. das Poten- 
tial dieser freien Wirbel zu berechnen, kénnen wir von dem Geschwindigkeitspotential eines Wirbels 
nach Formel (9) ausgehen. Dabei ist natiirlich sinngemaB a y(€, t) zu ersetzen durch 


ow 


Lor pees (KE) 


aes) : 
~i24l@. s e'ot__ 6 (Ee Ee) dX , 
0 
und auBerdem ist iber X von € bis 00 zu integrieren. Hierfiir benétigt man eine Reihe komplizierter 
Integrale, fiir die wir teilweise wieder neue Abkirzungen einfihren. Es ist 


” o o “oh ah 
nen 2)= | Me te eee eer Ae: oe 
g 
oh 
eS aie. wh ah a—€. wm, Wh e= | 
——te Et (24) + 260) 2 arctg*=* + iin 24 tn (1 + i 
o oh 
7 Saga) fea <1 1 /o\r : eae Bs 
ae G a [(h +ix—ié)*—h"] (17) 
i SA ge oe 2) il @ \n ja rm i én ( h)"] 
ge es SBP piesa, 
I; (x, 0, é) — 0 9 
sD i e—) 
lim J,(2, hk, €)=ae " e ™ , lim I,(x,h,)=0, 
CO) eeyax 2 ee 8, = or ot) | 
I,(x, h, €) = e ; @—XpP+R ee: ume ate AEN 
& 
oh 
— ain rh — &)? 
=e. Ga eniee (2) —21Gin S* arctg "5 * + Cof < In(1 pe 
0 
o wh o& 
fe fo, sp ge (18) 
+e e Paella} [(h bix—idr—hr] 
1 26-8 cy 1 fo (ht ie te —(— Ay", | 
a5, © e — n-n!\ up Ro 
oh 


lim I,(x,h,é)=ize “e ™ , lim I,(x,h,é) =0, 


z—>0o rE 
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es ewe fe ee > 
T(x, h, &, ty) = ba e ue 
3 gh 
M2 
yeaa X) 
— i> (#—X) 1 g x -g My g x 7 15 e* d? 
xX 4e 0 ie: (x— )) +i (M3 (x— +ig|+ a4 ie @— x 
u 
‘ 0 
eae i areas [elas 2 
PAG) : 
— ait * ao(® @—9)+15(2 @—9) +14 
as (19) 
& ide 
—i— x —6&) -— (uy — 2) 3 (M4 ) h 
My —_ 
AE é) 
bat) 
aay elke ae wighy? gh) 
sie e i (1 +( ua za) — ib arets (mF =] 
v 
oS 1 gh ata te i t)“11 d 
55 wat Sg tit — (7) v9 
n=1 
wo = h .o 118 z 
—i?@—s Bw—a -i2@-H —i4e—-§ 
bese h, &, fy) =o 5 : Up eXo 2a Ae uz ), 
fae 2, €>2,). — 0; 
i2 (x—f —u SF —y) in Se) 
e 2 e 49 en 8 g& 
W—Met2uQ+ 2 age ee 
— 0 
iG—O ay) é 
so ae T,(x, y> é) ys I(x, y, é) = U I(x, ¥, 9) s 3 (20) 
i —i 2 (x—9 —~Se—y in Se—x) 
e 2 e bef) a) g 
ve e—pt2pQat 2 ar 
: po 1 —ip—Q+1 i 
er, - I(x, ¥58) —"—Fo * Ley.) — 3H hts). (21) 


Fir den Beweis der Integralformeln (17) bis (21) wird auf Ziff. 10 verwiesen.} 
Mit Hilfe dieser und der friheren Integralformeln (10) bis (13) erhalt man nach einigen leicht 
nachzurechnenden Umformungen das Geschwindigkeitsfeld der freien Wirbel: 


uf(% ¥, Var” [AOI eee (AO TIR 


0 


1 Das in Formel (19) noch auftretende Integral mu8 numerisch ausgewertet werden. Dabei kann der fiir 
t = 0 logarithmisch singulaére Anteil des Integranden etwa in dem Bereich 0 <1 < 0, 02 direkt integriert 
werden. 
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mit 


Oe Gy, 2) hee ye) Fay, 6) TPE ey 8) 


2B 2B 
ae? 5) 
7 0 
get TNO) eee USO) arenes I(x, h — y, €, pg) 
Ay o g 22, 
— 1 fy Sey) 1 ( ral é) bala) a 
=3 ie A 2 e mee ¥ (sun, £)4 
0 (five S)s 
mam MO 1 (19 fos) 
me esr aa aN ae Sri ress fir x >é), 
0 (fiir x << &), 
Wee set 1 tot (TRE 
mr Y= sa Gee [AEC Ad tes eet [HOC a 
mit 
. : ip — 2 Sayers 
fo} st Le h—y,8) +i Tele h+y,8)— = FE * Bey, 8) +8 19,8) 
= t~) — oF hy) 
i I(x, y, Se eT (xh —y- £ Cis ee wal 
ae ol Bf é) == fa — Be (x nA e fy) + fg I(x, h Y> E, [s) (23) 
1 & é . oO : wes 
lang aa) i— («—é) i—> («—&) 
I ne 1 | Up 7A °. 
tae : gO aiaes (it te) 
0 (fiir x < €), 
L2§ - oO 5 ~ yo & 5 
— 5 (k—¥) mer Ae 19) Set — ¢) 
Ef u2 1 ( Ug u os 
agen = 2— Q St (ft 0 5), 
0 (fur 4=< 6), 


5. Randhedingung am Fliigel und Aufstellung der Integralgleichung. Wie iblich wird die 
Zirkulationsverteilung y(&, t) aus der Randbedingung am Fliigel bestimmt; diese besagt, daB die 
Strémung relativ zum Fliigelprofil tangential verlaufen mu. Nun ist zu beachten, daB durch die 
wellenférmige Anstrémung eine Anderung der Auftriebskraft und damit eine Tauchbewegung des 
Unterwassertragfliigels bewirkt wird. Dadurch fahrt der Fligel nicht mehr in der konstanten 
Tiefe y = — h, sondern es gilt y = 7(t). Da eine solche Tauchbewegung ein stabiler Schwingungs- 
vorgang ist, und da ferner die Amplituden der wellenfo6rmigen Anstrémung klein sind nach Voraus- 
setzung (s. Ziff. 2), wollen wir in unseren Formeln doch weiter eine konstante Tauchtiefe annehmen: 
n ~ —h, mit — h als zeitlichem Mittelwert von 7. Dagegen mu8 die Tauchgeschwindigkeit dy/dt 
auf jeden Fall in der Randbedingung am Fliigel beriicksichtigt werden. Der zeitliche Mittelwert der 
Auftriebskraft K, halt gerade der der Schwerkraft der Masse oberhalb des Fliigels das Gleichgewicht, 
wenn wir annehmen, daB das Boot ausgetaucht ist, und der Fliigel keine nennenswerte Masse hat. 
Wenn man von Dampfungseinfliissen des Wassers und anderen Stérungen absieht, ergibt sich dann 
durch die instationare Anstrémung im einfachsten Fall eine Tauchbewegung der Form 


d2 
Mo" ——Mg+K,(t). 


Fiir die hier durchgefiihrte Naherungsbetrachtung ist es ausreichend, zu setzen 
Ky(t) © @ up T(t) = @ Up dp + OW Ae?! + 0 up A, e—*! 
Ay = foé) dé, A, = f0,(6) dE. (24) 


mit 


1 Es sei denn, man beschrankt sich auf die Untersuchung von Fliigeln, die etwa in einem Priifstand fest 
montiert sind. Dieser Fall ist in der hier entwickelten Theorie mit enthalten; man braucht nur die Tauchge- 
_ geschwindigkeit wegzulassen. - 
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Mit M g = ¢ uy Ay und unter Beriicksichtigung der Periodizitatsbedingung 
22/w 
9 f 
a —nlt+o) gg | Wes 
0) 


liefert die Integration der Bewegungsgleichung 


3 A 
n(t) = —h — 8, eit ae he 
7 (25) 
dy => 8 Ay giot 1 8 At iat 
dt @ Ao w Ay 


Ist f(x) die Profilskelettlinie des Unterwassertragfliigels, so lautet die Randbedingung 


v(x, — h, t) + v,(x,—h, t) + ve(x, — h, i 


J AA u(x, — h, t) + uy(x, — h, t) + us(x, —h, 1) + Th 
Durch Einsetzen von u,, v,, Uf, Vp» Uy, U4, dy/dt gemaB (4), (15), (16), (22), (23), (25) ergeben sich 
durch Koeffizientenvergleich in 1, e'” ‘zwei Integralgleichungen fiir die Zirkulationsfunktionen 6,(&) 


und 6,(é). Auf die dritte Gleichung fiir 6, kann dabei verzichtet werden; denn mit 6, ist ja auch 6, 
bekannt. Bedenkt man noch, daB I,(x, 0, €) fir x = & singular wird wie In |x — é| (Enteptalecsinuay 
so nehmen die beiden Integralgleichungen folgende Form an: 


— 2 uy f(x) = + fee sagt Kl) as (-#s; <4), (26) 


jeg 


8 ; ans Tae 5 & oy pe —in Te -§ A, 
He oie Behe Rete oe +2 uae Ba (lL +is'@))e 0@ ° +215 = 


=% [50 io z+ h(x, &) — "eocem «fa (-ss : S«). (27) 


Die Funktionen K(x, £) und k(x, &) sind im ganzen Bereich in beiden Variablen stetig; es zeigt sich, 


daB die Integralgleichung (26) genau mit der Gleichung (I, 12) der stationaren Theorie iiberein- 
stimmt. Denn 


K(x, §) = Sie. tik sf@e “sin (@—8) +2050 Mp cos ig (x — &) 


2hg = 


= = (f'() — i) Le, — Ah, &, 1) — = (f(x) + i) L(x, —h, é, 1) (28) 


ist genau mit dem Ausdrucke (I, 13) identisch. 
Der Kernanteil k(x, €) wird leider sehr umfangreich: 


—2hf’ 1 
be) =" oepan 2g @ tS) [[1—i 5) he—hs) + (1 +15) Les — 1,9) 
Doe eae Iilscc obs & ne tee es 7) 
\ ah Bs E 2hg "8 
pen | 2m ote in het nt hee 
SK ee i a 


26 


ahg 2hg 
2 


Ty(x, 2h, &, uy) + 


+7, 1—if'@) a Lis, = hy 6) 2 ee eee | 


cae “al +i f’(x)) Pac I,(x, 2h, &, Be 
142 [709 4.282) — feet hs) — Hes even J 


| 
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o ee by —m as n ( i> @—g p28 Gee 
ee emcee ea Stig. ts (fiir 2 > 2), 
0 (fix x < é), 
a) : py Tee “8 a ( i= (x —8) i M5 Oe &) 
—2(+if@O)lgom? tale * gerne ieee 
0 (fur x < &) 


(29) 

6. Lésung der Integralgleichung. Wir werden hier nur die Auflésung der etwas komplizierteren 
Integralgleichung (27) behandeln. Denn die Auflésung der Integralgleichung (26) bzw. (I, 12) wurde 
bereits in Arbeit I durchgefiihrt und stellt auBerdem im Prinzip nur einen Sonderfall der Lésung 
von (27) dar. Da zunachst also die Gleichung (26) als gelést angesehen werden kann, ist’auch die 


Gesamtzirkulation A, (24) aus der AbfluBbedingung an der Profilhinterkante bekannt. Wir setzen 
nun 


%—=—acoss, £&=—acoso, 6,(é)= : P(c) (0s ¢ sa. (30) 


~ sing 


Mit Hilfe einer harmonischen Analyse mit der Schrittweite 7/6 stellen wir die stetigen Funktionen 
durch Fourierpolynome dar: . 


By 8s gre ge esa 
ay Se eae ae 2S ir — Pata ae fi ‘ 
deins (1+ fe) £ [Bian Be + Bapge a cae | = _ > F,sinds, (31) 
A=l 
5 6 
asins:-k(s,o)= > >» b,,sinAscosvo. (32) 
A 20 : 


Fiir die Berechnung der Fourierkoeffizienten b,, vergleiche man Formel (I, 18). Damit geht (27) 
in die Form iiber 
pe > F Aes 9%, See Mm Sia 
bi ene sg @ Ay 


erie a 5 6 
== | P(a)( sin s Ss SS bj, sin 2s cos. —i S* sins In [eos s — cos 0) do. (33) 
yi 0 


cos 6 — Cos s Sie 


0 
Fiir die Behandlung der Integralgleichung (33) ziehen wir die bekannte Darstellung’ 


frlcos:s*=cos. og) == Ina 2 ~ cosy s cosa (34) 
p= 
heran und machen den Lésungsansatz 
Hs Ps gies Ey A, 
P(o) = <2 + > Sip crs mit “P= 2 (35) 
|x Tyal a \x 
Damit ergibt sich zur Bestimmung der Lésungskoeffizienten P,(v = 1, ... M) das folgende lineare 
Gleichungssystem : j 
ee 26 .oal sy mag fl A, .oa 
ice irae ree =F,+ y2mi i Fe tu Up n2), 
A, 


rs Ee aT Ree ee fee en ee tag 


efuell re) 76, eee ie (6) Lee) 68 OT ee 


In allen praktischen Beispielen diirfte es ausreichend genau sein, bereits P, = 0 zu setzen, also 
M = 6, so daf sich das an sich unendliche Gleichungssystem (36) auf ein System von sechs Gleichun- 
gen fiir die sechs Koeffizienten P, reduziert. Fast immer wird w a/uy< 1 sein. Die Gesamtzirku- 


; 1 W. Schmeidler, Integralgleichungen S. 67. Leipzig 1950. 
12 
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lation A, wird wie iitblich aus der AbfluBbedingung an der Profilhinterkante bestimmt, d. h. aus der 
Gleichung P(x) = 0. Damit ist die Integralgleichung (27) bzw. (33) gelést. 


7. Berechnung der Fliigelkrafte. Die momentanen resultierenden Fligelkrafte in x- und y- 
Richtung K, und K, werden prinzipiell genau wie in Arbeit I, Ziff. 5 nach dem Kutta-Joukowski- 
schen Satz bereclnee! Dazu denken wir uns den Fligel durch einen im Druckmittelpunkt (d. h. 
etwa im 1/4-Punkt) liegenden Punktwirbel gleicher Gesamtzirkulation ersetzt: 


I) = f v(t) dé. 
Dann wird mit der Wasserdichte 0 
K, =e co + uy, (= $3 —hat) + us ed —h, t) + “(3 = shyt) I(t), 
itso at Grapes: t)+ vt pg goa + o(—5>—h 4] T(t). 


h, t)sind aus den allgemeinen 


(37) 


Die Geschwindigkeitskomponenten u; = > ,—A, i), vy i ce 


Ausdriicken (15), (16) ohne weiteres zu bilden, indem statt der Wirbeldichte y(€, t) ein bei € = — a/2 
gelegener Punktwirbel /'(t) angesetzt wird. Der Stern bedeutet, daB natiirlich die vom gebundenen 
Wirbel selbst stammenden Anteile wegzulassen sind. 


Die Berechnung von uy (- + »—h, ) nach (22) macht ebenfalls keinerlei Schwierigkeiten, nur 


bei der Bildung von vs |— --, —h, t) nach (23) bedarf der von dem Glied I,{— > = (UF a) stammende 


Anteil einer besonderen eee Wir haben zu untersuchen, welcher Wert physikalisch sinn- 


gemaB fir den formal unendlich werdenden Grenzwert lim I, (* 0, rd einzusetzen ist. Da 


x—>—a/2 


ua s}e 50a le+s) +49 +08 +8) 


ie 


a 


ist, wird 


lim I(x, Tees 3)= =i 5+ 05772 + In 22 tin Th 


x—>—a/2 Uo x—>—a/2 
Dieser unendliche Grenzwert beruht darauf, da8 wir bei der Kraftberechnung von der Wirbeldichte 
zum Punktwirbel iibergegangen sind, d. h. wir haben den eigentlichen Ausdruck 


J y(% t) vs (x, — h, t) dx 


ersetzt durch 
4(—F5—h ) I(t). 


Dieses ist ja bei stetigen Geschwindigkeitsanteilen eine gut brauchbare Naherung. Es liegt daher 
nahe, fiir den singularen Anteil zu schreiben 


=T9 fre t) In ress ad. 


Dabei ist es fiir diese Naherungsbetrachtung méglich, 


lim In 
x—>—a/2 


P 2 


It — 
y@r,t) =O /e—* 
zu setzen, d. h. nur den ersten Birnbaumschen Anteil zu nehmen. Dann wird wegen (34) 
. a SS x x | 
jim In) 5 Zz aS an alee a+x dx 
-ifa + cos s) In + —c0s s dsi= Sn. 
fh 


teeter nee & 
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Damit wird der fiir die Kraftherechnung gesuchte Wert 


L(+, Ores Z| = iF + 0,0772 + In 2. 

8. Die 1/4—3/4-Punkt-Methode. Um ohne die Auflésung der Integralgleichung naherungsweise 
die Gesamtzirkulation J’ des Fliigels zu bestimmen, hat sich in der stationaren Theorie des Unter- 
wassertragfliigels (Arbeit I) die bekannte 1/4—3/4-Punkt-Methode gut bewahrt. Es liegt daher 
nahe, sie auch bei dem vorliegenden instationaren Problem anzuwenden, um bei N aherungsrech- 
nungen die Auflésung der Integralgleichung zu umgehen. Die beiden Bestimmungsgleichungen fiir 
A, und A, ergeben sich unmittelbar aus (26) und (27); sie lauten 


tur (f)=4[t-+K(5-— 3) 


gh : ; a 
SPY He et g se ey 23 5 Tes ee 
2am B(1 +if'(2)) 6 “he 7 2a py & BL +if (S))e be 2h (38) 


9. Die Form der Wasseroberfliche. Zur Berechnung der Form Y(x, t) der freien Wasserober- 
flache dient Gleichung (1), die wir genauer in der Form schreiben 


1 ee pee saa ) ap ae 0, t) + u,(x, 0, t) + u4(x, 0, 1)]. (1) 


Y(x, t) =—— 

fo) 

Dabei kann der Unterwassertragfliigel durch einen in seinem 1/4-Punkt angebrachten Punktwirbel 

gleicher Gesamtzirkulation ersetzt werden. Aus Arbeit I ist bekannt, daB dieses Verfahren prak- 
tisch ausreichend genau ist. 

Der Ausdruck 0@,/dt in (1) bezieht sich auf die freien Wirbel. Unter Benutzung unserer Ergeb- 

nisse fiir D,, u,, us gemaB (9), (15), (22) erhalt man nach einigen elementaren Zwischenrechnungen 


: gx é od 
7 2 SSD: See lal va 
¥(x,1) =— e§#) (u, — Q) Bye 0 + (Uz — 2) Bye 
: ip, eo Ginee 
lg Oe (GLO Ee" 
apis My 1 Me 2 
ie Ao 1,[20, ¢.1)+4(x0, $,1)—206 “0 sin (e+2)| 
Tt 
KA ORE ee Gere 
ee ete t } 
mit ; : f : 
- & 1 
f y= z(1—#g) B(x z)+3(1 ! iB) Len 9-3) + L(0 ol) 
gonad) 
M->G My x + 
M,— 2 ee: eee —— fe 7) u 2 
mee a 1,(x,0, Be) Bo acai * 0e ss 
h 
— Pa 2 x + —r ae) 
Fa (a 
gh J 
PO ator a My— 2 
neers + 1,(s0, 5H) H2— My 
h o a 18 is js a 
Gua eal eee ‘ p+3)) (rir 2 >—$), 
+ Pe tee 122 
0 (fi zi =F), 
gh eeepc x z —t vet ee 
in Q seat = 5: 7) e E a 2) (fir « >—$), 
+ Mo — fs 
0 = See 


12* 
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Mit Hilfe der Integralformeln (10) bis (13) und (17) bis (21) erkennt man leicht, daB lim Y(z, t) 


Sar eae, 0) 


nur durch die wellenférmige Anstrémung bestimmt ist, wie es auch sein muB. Der Wert von lim 
x—> CO 


Y(x, t) ist ebenfalls ohne weiteres aus (39) abzulesen; jedoch kommt d’esem keine physikalische 
Realitat zu, denn die freien Wirbel bleiben ja in Wirklichkeit nicht in unveranderter Starke hinter 
dem Fligel erhalten, wie wir es hier voraussetzten, sondern sie zerfallen allmahlich. Gleichung 
(39) gilt deshalb fiir x > 0 nur in der Umgebung des Tragfliigels. Fir einige grundsatzliche 
kritische Betrachtungen zu der Formel fiir die Form der freien Wasseroberflache kann auf 
Arbeit I, Ziff. 6 verwiesen werden}. 


10. Beweis der Integralformeln. Wir haben nunmehr den Beweis fiir die Integralformeln (10) 


bis (13) und (17) bis (21) nachzutragen. 
a) Beweis von (10), (11). Wir fahren den Beweis fiir Formel (11); ganz analog ist er fiir For- 


: vs es & : 
mel (10). Im Folgenden sind A und H Abkirzungen fir uw (« — &) und we (h—y). Es ist 


| Pe — eH 1941C (OA) +iS8(QA) =i Z| | 


; B+ 2 
— e2H—iQA | (1 — ¢— 2H) nar et, A>? 
ro t-A=0 
1 
mit 
TE ee [cet CL OA; = C(O A) 
und eet | 
: d 
s(@ 4) = { sine = ; Si 0 A) = BAe | 
oy F | 
Aus | 
[ rortiedear = EERE EA: 
- O—iQA 
folgt durch Integration iiber } 
r d ro O Q 
2 gti It) oe Ai eee 9 OTIiQA 
Je e 2 yer s Beg RO re ers 


0 


Damit ist der erste Teil von Formel (11) bewiesen. Es gilt nun 


2H OLIQA Aue QH 
iQA —? E iif Pt eA (O14 do 
Je reoa =| ame Ot | Cet hn oe 
iT) 0 
1, 44H. H Cos de ee 
= n—],— +i aretg 5 + S' — [§ Q4—QHY— (QA). 
me 


Hieraus ergibt sich die zu beweisende Formel (11), wenn man noch bedenkt, daB 


H I A A ) 
arctg — = ++ — — arct oe 


— und ferner fiir C,(2 A) und SQ A) die bekannten Reihenentwicklungen einfihrt | 


(QA)? | (Q A)! 
2.21 1 SS Te 


C(QA) = 0,5772 + S in. 4a — 


+ Der zeitunabhangige Anteil von (39) stimmt -natiirlich genau mit der Oberfachenform der stationaren 


2 (2.4) : 
S(2 A) = QA—Sat+eee. , 
Theorie nach Gleichung (I, 1), (I, 34) tiberein. — 
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b) Beweis von (12), (13). Wir fiihren den Beweis fiir Formel (12); ganz analog ist er fiir (13). 
Dabei kénnen wir A > 0 voraussetzen, da der Fall A < 0 durch den konjugiert komplexen Wert 
der Formel (13) erfaBt ist. Es gilt 


foo) fore) o 

i a “ e Htt+iAr e Hrt+iAr 
a aE aa 

A eu Ft a! ie ’ t+ 4B F t+iB 


Die Integration iiber den Viertelkreis des ersten Quadranten der (t + i 1’)-Ebene liefert nach der 
Residuenmethode fiir unendlichen Kreisradius 


co 


fore) 

e Hrt+iAr e tH’ —Ar’ 

——_— dt = dt’ 
Tee ; ve +B 

0 0 


— — ef 4+ipH o.s772 +5 In (42 + HP) +> Inf? + i arctg + 5 


! 
=“, n-n! 


—iBH—Bf A)", 


wobei die letzte Umformung auf (11) beruht. 


Ferner ist 
f Hr +iae iHp+ AB r : ig ip 
== (i A — H) (rt + 
| ae. v a hl 
1 1 2 (i A— Hn ey 
= ob 44101 In (at +B?) —y In B — i aretg 5 Ye: = [(« +78)" — (iB) I. 
Also wird 
fees e =, (H —i A) (« + if) 10.5772 Ages F (A? + H?) 
weatip (at ee Pas : 2 


+ In (a? + 2) + i arctg & —i aretg F =e (A — Hy" («+1 6” 
= — e#—i4)(« + #8) F(A — H) (a +7 8)), 


wie man unter Benutzung einer bekannten! Darstellung fiir die E,-Funktion mit komplexem Argu- 
ment erkennt. Es gilt nunt 


. (i A — H) (a + if) 
pp er UN) a, Um G@ A—H)(@ +38)’ 


und damit ist Formel (12) vollstandig bewiesen. 


c) Beweis von (17), (18). Wir fiihren den Beweis fiir Formel (17); ganz analog ist er fiir For- 
mel (18). Es ist 


foe) oo h 
=) 
.o 3 a) ae 
rae At *) h dX i Bae oe e o de 
Se ie Fl 
E Seems 
h 


00 A/H 
—inHr 
ariel fe ac +{ ei QHr o ' (gilt auch fir A< 0). 
T 
0 


Die Integration iiber den Viertelkreis des ersten Quadranten der (t ++ it’)-Ebene liefert nach der 
Residuenmethode bei unendlichem Kreisradius 


foe) co y. 
tQHr — 2 1 
| aes ee tae 


0 


1 Vgl. z. B. Lésch-Schoblik, Die Fakultat und verwandte Funktionen S. 118—122, 1951. 
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also wird? 

me H p — 2Hr PO HY i 2 db 
{OH i — 1 bsp ale mee = 0 ei — oT 
e 28 Oo OH te ey ee a dr’ = —e 2&4 +—e [e 9 

|Gare-¥ ae ae iy , eeely 2 2 ce 


— + eH : 3a e- 91 Es(Q H) + 4 °F E(—QH). 
H 


Weiter ist 
A/H A/H A/H 


On : 
[ eon dt cae {a dpa a | GOR ee 
6 


wisp I 21 


A|H aS 


QH . QH $ dt 
fe ee pees iQH(e+i)__] 
atte | ie 25 [ ¢ \aee 
0 


ee ee AB Pt ee a | ; : . 
cee! 2 DawWlOHt+i@d) — (2 H)") 


A 
= Co) 2 H- arctg it 


oo 


QH+i QAy— Cle. 


Damit ist Formel (17) bewiesen; die Richtigkeit des angegebenen Grenzwertes fiir x —> -L0O ent-' 


nimmt man leicht aus der obigen Herleitung. 
d) Beweis von (19). Es ist 


Rae ey i Aleit i 


C(t) — 1 S(t) +t > dt 
ae (x 
yee 
r ee ene e O+it 
= | si | Pegs 2 UO dt 
— "5 @—8) : 
0 
Zunachst ist (A beliebig) 
oul tek Ve i) eS ae eee 
io = Je) 1840) +i Fae = — | ma a [ <woe 
—mA =A Tt 
: Q r : im ra . 2 
a i py Boe face ify a eae dt 
aloe Deel aati c leery ee 
id —,A —p,A 
at i fy = 10) A S oer Ps 4 Bb 1 My . - 
ear ied C(uy A) +i Sy A) +i | Hs CMA) +i S(2.4) +i 
Weiter haben wir 
=f on WH ee 2 A MA 2 4H 
ee ek i ir GES 
ie Je pepe ddr = fe| eye woos rao [er i eae ae WO de. 
WA oo 9 $b o=o0 


1 Die beiden reellen Funktionen E*(x) = — | e? ae und E,(—x) =—, | Cae ae 


bel Jabuke-Kinde; Tafel hohe er ae unlaiouew = 


dah aca ala a ea, 


PUR ie xP coms, 
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Nach Formel (13) wird nun 


co 


ete} oo 22 Q H 
\ es 
fe ste fo tigate sm on Bt te 8 fey 
0 0 u 
oe 
Ie Lae Sine i E* 28 
My 
und daraus folgt weiter 
4H co My 1 Q My 5 
yi Se Sees. 
af e.’ motte aa aal a i) wife z) [ et a0 
o=0 0 _ Re 
My 
Bj oe) “oH 4H 
: foe 3 
ee oi Hyu,—HaQ Uy My eeady i fy dé 
se (LE ==5)| a Mok a 
| Be ii 0 4 Je x See : o 
2d 0 


i ln BO — i HO BS(O H) +i BM (uy HW). 


SchlieBlich ist noch zu berechnen (A beliebig) 


fy 


[yA 2 fy H [4 A 2 ) fy 
ya o—it ee ey o—it dd 
by o Set ly = oO a Res. 
[- | e fp war =| dt as (e? + Noth 
Ss 0 0 . 0 


2Q 
eA ceyey hA ae ; H <1 (u, H oe 
Feta tyme Siac 
0 


n-n! 
n=l 


Durch Zusammenfassung der einzelnen Integralanteile ergibt sich die zu beweisende Formel (19). 
Auch die in (19) angegebenen Werte fiir lim J, sind aus deni in der obigen Herleitung enthaltenen 
Formeln ohne Schwierigkeit zu entnehmen. 


e) Beweis von (20), (21). Wir fihren den Beweis fiir Formel (20); ganz analog ist er fiir For- 
mel os Wir schreiben zundchst unter Vertauschung der Integrationsreihenfolge 


ey 72 x —é 
f 29 a ie ae en the a eral 
. bat) g Laie: ss = Ze = 
ae are aay erie i} @+2)(W—n+2uQ+ O) 
5a : 
BS i X)—x 
te et C8) Bee" a ae 
eee (u+ Q)(W— w+ 2H Q+ 2%)" 
0 
Da 
bed —e@-ip—2 il a Ope 1 <a 1 


ist, folgt nach (11), (12), (13), daB das X, enthaltende Integral fiir X)—> co verschwindet, und es 
ergibt sich genau Formel (20). 


(Eingegangen am 20, Juli 1959.) 


Aus dem Institut fiir Angewandte Mathematik und Mechanik 
der Deutschen Akademie der Wissenschaften zu Berlin 


Anschrift des Verfassers: Dr.-Ing. habil. W. H. Isay, Berlin-Dahlem, Schweinfurthstr. 90. 
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Beitrag zur Berechnung des Warmeiiberganges und Druckabfalles 
laminarer Einlaufstrémungen. * 


Von Karl Stephan 


1. Einleitung. Die theoretische Grundlage von Problemen der konvektiven Warmeiibertragung 
sind bekanntlich die Fouriersche Energiegleichung! und die Navier-Stokessche Bewegungsgleichung* 
aus der sich fiir sehr kleine Zahigkeiten die wesentlich einfachere Prandtlsche Grenzschichtgleichung? 
herleiten laBt. Im Fall temperaturabhangiger Stoffwerte kénnen diese Gleichungen nicht unab- 
hangig voneinander gelést werden, weil das Strémungsfeld auf das Temperaturfeld einwirkt und 
umgekehrt eine Riickwirkung vom Temperaturfeld auf das Strémungsfeld besteht. Der mathemati- 
schen Behandlung des Problems stehen betrachtliche Schwierigkeiten entgegen. Aus diesem 
Grunde ist bisher noch keine umfassende Lésung iiber den Einflu8 der Temperaturabhangigkeit der 
Stoffwerte bekannt geworden. 

Sofern jedoch die Temperaturunterschiede innerhalb des strémenden Mediums nicht sehr groB 
sind, kann die Zahigkeit konstant gesetzt werden. Das ist bei einer Reihe von Stoffen recht gut 
erfiillt.4 Ist ferner die Dichte konstant, dann entfallt die Riickwirkung des Temperaturfeldes auf 
das Strémungsfeld. Man hat also zuerst die fir das Strémungsfeld maBgebenden Gleichungen 
zu lésen und kann ausgehend von dieser Lésung das Temperaturfeld berechnen. 

In der folgenden Untersuchung wird auf Grund einer solchen Rechnung die im Einlauf von Roh- 
ren und von ebenen Spalten iibertragene Warmemenge bestimmt. Gleichzeitig wird der Druckab- 
fall ermittelt. Es wird Laminarstrémung vorausgesetzt. Neben der Zahigkeit sollen auch alle tib- 
rigen Stoffwerte nicht von der Temperatur abhangig sein; die Reibungswarme und die Warmeleitung 
in Strémungsrichtung werden vernachlassigt. Es wird angenommen, eine Fliissigkeit stréme aus 
einem gréBeren Raum in ein Rohr oder in einen ebenen Spalt. Dann ist die Geschwindigkeit aller 
Flissigkeitsteilchen tiber den Eintrittsquerschnitt des Kanals konstant. AuBerdem sollen die 
Temperatur der Kanalwand und die der Fliissigkeit im Eintrittsquerschnitt konstant, aber vonein- 
ander verschieden sein. 


2. Physikalische Beschreibung des Problems. Auf Grund der Vorstellungen der Grenzschicht- 
theorie gelangt man zu einer physikalisch anschaulichen Deutung der Vorgange im Einlauf. Tritt 
die Fliissigkeit in den Kanal ein, so werden die Teilchen in unmittelbarer Wandnahe unter dem 
EinfluB der Reibungskrafte verzigert. Aus Kontinuitatsgriinden erfahrt die praktisch reibungsfreie 
Strémung in der Mitte des Kanals, die sogenannte Kernstrémung, eine Beschleunigung, welcher 
nach der Bernoullischen Gleichung ein Druckabfall entspricht. Stromabwarts werden immer 
mehr Teilchen von Reibungskraften erfaBt. Die Grenzschichtdicke sowie die Geschwindigkeit der 
Kernstrémung nehmen zu. Das Geschwindigkeitsprofil geht schlieBlich asymptotisch in die 
Poiseuillesche Parabel iiber. 

Entsprechend der Strémungsgrenzschicht bildet sich eine Temperaturgrenzschicht aus, die 
ebenfalls mit der Lauflange anwachst. In hinreichender Entfernung von Eintrittsquerschnitt andert 
sich auch die Temperatur der Teilchen in der Mitte des Kanals. Das Temperaturprofil geht asym- 
ptotisch in die konstante Temperatur der Kanalwand iiber. 

Man kann den Einlaufvorgang als einen AusgleichsprozeB auffassen. Dieser fihrt zu einem 
Gleichgewicht, das erreicht ist, wenn das Geschwindigkeitsprofil in die Poiseuillesche Parabel tiber- 
gegangen ist, und wenn sich die Temperatur der Flissigkeit nicht mehr von der der Kanalwand 
unterscheidet. Das Gleichgewicht wird um so rascher angestrebt, je weiter man von ihm entfernt ist. 
Ganz vorn im Einlauf wird daher mehr Warme je Flacheneinheit iibertragen als weiter hinten. 
Ebenso sind die Wandschubspannung und damit auch’der Druckabfall in der Nahe des Eintritts- 


querschnitts gréBer als weiter stromabwarts. Der bessere Warmeiibergang von Einlaufstrémungen 
ist also mit einem héheren Druckabfall verbunden. 


* Auszug aus der von der Fakultat fiir Maschinenwesen der Technischen Hochschule Karlsruhe genehmigten 
Dissertation des Verfassers. 

1 J. B. Fourier, Théorie analytique de la chaleur, Paris 1822. 

* M. Navier, Mem. d. L’Acad. d. Sci. 6, (1827); G. G. Stokes, Trans. Cambr. Phil. Soc. 8, (1845). 

° L. Prandtl, Verh. d. III Intern. Math. Kongr. Heidelberg 1904, 

* G. H. Géttner, Erdél und Kohle 3 (1950), S. 598. 
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Trotz der einfachen physikalischen Modellvorstellung des Kinlaufvorganges existieren indessen 
nur wenige Untersuchungen, die sich mit der Uberlagerung von hydrodynamischer und thermischer 
Kinlaufstrémung befassen. Es wurden lediglich eine Reihe von Naherungslésungen bekannt!?:3.4:5, 
die vorwiegend in einiger Entfernung vom Eintrittsquerschnitt und fiir solche Stoffe gultig sind, 
deren Prandtl-Zahl Pr ~ 1 ist. Die Prandtl-Zahl ist hierbei definiert durch Pr — fue,/A mit der 
dynamischen Zahigkeit uw, der spezifischen Warme bei konstantem Druck c, und der Warmeleit- 
zahl 4. Neuerdings hat W. Geck® ein Differenzenverfahren zur Lésung der Prandtlschen Grenz- 
schichtgleichung und der Fourierschen Energiegleichung angegeben. Er berechnete den Warme- 
tibergang und den Druckabfall im Einlauf luftdurchstrémter ebener Spalte. Fiir Luft ist Pr ~ 0,75. 
Das ziemlich miihsame Differenzenverfahren erfordert jedoch einen kaum zu bewaltigenden 


Rechenaufwand fiir Pr > 1, weil man dann eine auSerordentlich kleine Schrittweite wahlen muB. 


3. Mathematische Formulierung des Problems. Es werden folgende Gréfen eingefihrt: 


x = Koordinate in Strémungsrichtung, gerechnet vom Eintrittsquerschnitt an; 

y = Koordinate senkrecht zur Strémungsrichtung, gerechnet von der Kanalwand an; 
u = Geschwindigkeit in Strémungsrichtung; 

Us = Geschwindigkeit der Potentialstrémung; 

v = Geschwindigkeit senkrecht zur Strémungsrichtung; 


U,)= konstante Geschwindigkeit im Eintrittsquerschnitt ; 
= halbe Spaltbreite oder Rohrradius; 

= Dichte; 

= kinematische Zahigkeit ; 

= Temperatur eines Fliissigkeitsteilchens ; 

t,, = konstante Temperatur der Kanalwand; 

t) = konstante Temperatur im Eintrittsquerschnitt ; 

p = Druck eines Fliissigkeitsteilchens. 


eo) 


Fiir die Rechnungen erweist es sich als zweckmaBig, dimensionslose Gréfen einzufiihren. Wir 


definieren die Reynoldsche Zah] durch 


Recess 
Vv 
und es sei 
u ug v 
——— = — Vie=— Re; 
PERL da ai U, 
st ataee = (eS 
X= Re’ Sarg, ty — tw 


Die exakte Beschreibung des Warmeiiberganges erfordert die Berechnung des Temperaturfeldes, 
die wiederum die Kenntnis des Strémungsfeldes voraussetzt. Fiir die Geschwindigkeitsverteilung 
der ebenen Strémung ist die Prandtlsche Grenzschichtgleichung maBgebend, die unter der Voraus- 
setzung kleiner Zahigkeiten gilt. Sie lautet mit den oben definierten dimensionslosen Gréfen 


aU BU. =: dls *6°U : 
Og gy 8 axe Pays (1) 
Ferner eat die Kontinuitatsgleichung 
tan? (2) 
OX GON 


erfiillt sein. Die Temperaturverteilung der ebenen Strémung berechnet man aus der Fourierschen 
Energiegleichung 
oT Fy re BC § 


3 
: koxe ok at ate 612 8) 


1 §, Sugawara, T. Sato, Mem. Fac. Engn. Kyoto Univ. 16 (1954), S. 166. 
2 K. Elser, Allg. Warmetechn. 3 (1952), S. 30. 

3 H. Hahnemann, L. Ehret, ZMB-Forsch. Ber. Nr. 1751 (1943). 

4 W. M. Kays, Trans. ASME 77 (1955), S. 1265. 

5 J. A. Businger, Ingenieur 64 (1952), S. 17. 

6 W.Geck, Diss. T. H. Karlsruhe 1953. 
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Hierin sind U = U(X, Y), V = V(X, Y), Us = U(X), T = T(X, Y). Die Randbedingungen 
lauten 


X>0, Y=0 wud Y=2:U=V=0, T=0, 
Xo 0 pes SO fe Ler feds 


Wie H. Schlichting! gezeigt hat, gelten die Gleichungen (1) bis (3) nicht nur fiir die ebene Stré- 
mung. Sie bleiben auch noch fiir gekriimmte Wande giiltig, solange die Grenzschichtdicke gegen- 
iiber dem Kriimmungsradius der Kanalwand klein ist. Diese Voraussetzung trifft vorn im Einlauf 
von Rohren und auch von Ringspalten zu. Dort darf man also zur Berechnung der Grenzschichten 
ebenfalls das System der Gleichungen (1) bis (3) benutzen. 


4, Lésung fiir das Strémungsgeschwindigkeitsfeld. Die Integration der Gleichungen (1) und (2) 
erfolgt nach bekannten Methoden der Grenzschichttheorie®**, Diese bestehen darin, da man 
eine Stromfunktion Y/(X, Y) einfihrt, welche die Kontinuitatsgleichung erfillt. Es soll also 
gelten 
(X,Y) __ 


= O10. 6 und =k V(X, Y). 


Fiir die Stromfunktion macht man einen Reihenansatz 


P(X, ) = foln) XU? + uy film) X + us fo(n) X97 + --- (4) . 
it 
m1 - 
Jims 5 Cee 


In 7 ist die Koordinate Y durch den Wert 2 X1/? dividiert. Dieser ist der Grenzschichtdicke 6,(X) 
der langsangestrémten ebenen Platte proportional. Die f;(7) (i = 0, 1, 2, ...) sind die sogenannten 
Koeffizientenfunktionen des Geschwindigkeitsfeldes. Die u; sind noch naher zu bestimmende Kon- 


stanten. 
Bildet man in Gleichung (4) 


rae UX Y) = [foln) + uy film) XH? + we fa(n) X +--+], (5) 


so erhalt man eine Aussage iiber die Potentialstrémung U,(X). 

Fiir hinreichend grofe Werte 7 muB U(Y, Y) in die Potentialstrémung U,(X)' iibergehen. 
Gleichung (5) mu also fiir 7->0o unabhangig von Y werden. Das ist nur méglich, wenn die 
fin > 0c) (i = 0, 1, 2,...] konstante Werte annehmen. Es ist zweckmaBig filn—> co) =2 zu 
setzen, weil dann der Faktor 1/2 in Gleichung (5) wegfallt. Damit lautet der Ansatz fiir die Potential- 
stroémung 


U(X) yp RO ee yg RE (6) 


Zur-Lésung der Prandtlschen Grenzschichtgleichung (1) setzt man die Ansatze (4) und (6) in 
Gleichung (1) ein und nimmt die erforderlichen Differentiationen vor. Setzt man die Glieder gleich 
hoher Potenzen von X gleich null, so erhalt man ein System von unendlich vielen gewohnlichen 
Differentialgleichungen. Dieses enthalt noch die Konstanten u;, von denen man sich befreien kann, 


wenn man entsprechend dem Vorgehen von Howarth® weitere Koeffizientenfunktionen einfiihrt. 
Es sei 


faln) = 820) + 2 hgln) 


Uy U2 


fa(n) = 83(n) + 


us 
h(n) + = h(n) 


Ug 


u uj 


Fale) = eat) + Sn) + hale) + E28) +h mln) 


U, u, 


1 H. Schlichting, Grenzschichttheorie. S. 106—107, Karlsruhe 1958. 
? H. Blasius, Z. Math. Phys. 56 (1908), S. 1. 

° L. Howarth, ARC Reports and Memoranda 1632 (1935). 

4 X. Hiemenz, Diss. Gottingen 1911. 

° L. Howarth, Proc. Roy. Soc. Lond. A 164 (1938, S. 453. 
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An die Stelle der Prandtlschen Grenzschichtgleichung (1) tritt dann das folgende 


: | rekursiv zu 
lésende Gleichungssystem mit der unabhangigen Veranderlichen (ae 


fv +fih =. (7) 
A +hifo—fifo t2Ah =—4, (8) 
By +85 fy—2efo +3 afi =—8, 
ee ee ©) 
8 +83 fo—3esfo+4esfo =—12, 
hy’ + hy fo— 3h fy + 4hg fo =—12 +3 fig; —3f% a —2f, 2%, (10) 
key’ + hy fo — 3k fy + 4 hg fy = 3 fi hy — 3. hg — 2 f, WY, 
G+ Siho—4eifo +5 efo =—16, ] 
Wy + fo—A Ry fy + 5 hy fe =—16— 2 f, + 4g, —Aaaft | 
ky + kg fo—4 ky fo + 5 ka fo =—8—3 8) +282, \ (11) 
Wt Me fo— 4s fo +5 le fo = — 2M fy + 4M fi — A fy — 315 fe + 4 — Bh ge | 
mg + my fo— 4 mi fy + 5m fo = —2 hy fy + 4 ky fy — 4 hey fy — 3 hig hy + 2 AZ, | 


Das System dieser Gleichungen ist fiir die folgenden Randbedingungen zu lésen: 


1. Wegen Y (X, Y = 0) = 0 verschwinden alle Koeffizientenfunktionen an der Stelle 7 = 0. 


2. Wegen U (X, Y = 0) = 0 und V(X, Y = 0) = 0 verschwinden die ersten Ableitungen der 
Koeffizientenfunktionen an der Stelle 7 = 0. 


3. Fiir sehr groBe Werte 7 geht die Grenzschicht in die Potentialstrémung iiber. Es gelte fiir 
1 —> OO 


~ 


fo =fi = 8 =8 =& =: =2, 


Gleichung (7) tritt bei der Berechnung der Grenzschicht einer langsangestrémten ebenen Platte ' 
auf. Sie wurde zuerst von H. Blasius! durch Potenzreihenansdtze und spater von C. Tépfer? und 
L. Howarth? nach dem Runge-Kutta-Verfahren numerisch geliést. Nach dem gleichen Verfahren 
berechnete L. Howarth ferner die Koeffizientenfunktionen g,, wahrend B. Punnis* die Gleichungen 
fiir f, und h, sowie die Gleichungen (10) nach dem Verfahren von Adams-Stérmer vertafelte. 


Da die Koeffizientenfunktionen bei der Berechnung des Temperaturfeldes von Bedeutung sein 
werden, kommt es darauf an, sie méglichst genau zu kennen. In der vorliegenden Arbeit wurden 
daher samtliche Koeffizientenfunktionen der Gleichungen (7) bis (10) erneut berechnet. Die numeri- 
sche Rechnung wurde von der elektronischen Rechenanlage der Farbwerke Hoechst nach dem 
Runge-Kutta-Verfahren ausgefiihrt. Die Zahlenergebnisse sind in der Dissertation des Verfassers® 
aufgefiihrt. 


Die Koeffizientenfunktionen des Systems der Gleichungen (7) u. f. haben universellen Charakter. 
Sie sind unabhangig von der Geometrie des durchstrémten Kanals. Diese wird erst durch die Kon- 
stanten u; beriicksichtigt. B. Punnis* hat in seiner Arbeit ausfihrlich beschrieben, wie die u; 
berechnet werden. An dieser Stelle seien daher lediglich die Zahlenwerte fiir die u; mitgeteilt. Sie 
weichen etwas von den Werten ab, die B. Punnis* ermittelte. 


| Mt ; | & | e | ss 


ebener Spalt 1,7208 — 22,2494 19,823 — 208,41 
Rohr 3,4416 — 13,3347 | 206,252 | —4709,92 


1 siehe FuBnote 2 von S. 178. 


2 C, Tépfer, Z. Math. Phys. 60 (1912), S. 397. 
3 siehe FuBnote 3 von S. 178. 

4 B. Punnis, Diss. Gottingen 1947. 

5 K, Stephan, Diss. T. H. Karlsruhe 1959. 
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5. Lésung fiir das Temperaturfeld. Es ist naheliegend, die Methoden, welche sich bei der Be- 
rechnung der Strémungsgrenzschicht als fruchtbar erwiesen haben, auf die Berechnung der Tem- 
peraturgrenzschicht zu tibertragen. Man macht demnach fiir die Temperaturverteilung den Ansatz 


T(X.'n) = Tog) 2 Ti) Ee) ee a) eee (12) 


Darin sind die T;(7) (i = 0, 1, 2...) die Koeffizientenfunktionen des Temperaturfeldes. Ahnlich 
wie beim Strémungsfeld werden die Koeffizientenfunktionen in folgender Weise weiter aufgespalten: 


T,(y) = u, H,(7) , 

T,() = Uy Go(y) + uj H,(n) » 
T3() = ug G3(n) + uy Ue H3(y) + ui K;(7) » (13) 
T4(y) = ug Gy(n) + wy ug Hy(n) + ut Ky(m) + uj uy L4(y) + uy M,(n) 


Nach Einsetzen von (12) mit (13) und von (4) in die Fouriersche Energiegleichung (3) erhalt 
man durch Koeffizientenvergleich das folgende rekursiv zu lésende System gewdhnlicher linearer 
Differentialgleichungen, in das nun auch die Koeffizientenfunktionen des Geschwindigkeitsfeldes 


eingehen: 


1 at / 
Pr fis =P fo To = 0 2 () 
a 


1 ZA is a, 
—pp Hi + fo th —fo Hi = —2f, TM, (15) 


1 fr f - / 
p, & + fo Go —2 fo Gg =—3 Th, 


1 Vt 7 7 te 7 la 
~p, He + fo Hy — 2 fo Hp = —2 Wf, + H, fy —3 The, 


(16) 


i I} la / 7 

—p, &s + fo G3 —3 fo Gs =—4 Tog,» 
1 wt / , / 77 / 77 y ; 
~p; As + fo Hs — 3 fo Hs = —2 fi Gy + 2 fi Go—3 Ay 8. + Hy g2—4 Ty hg, (17) 


1 aa 2 , rs a Ya 7 , , 
-p; Ks + fo Ks — 3 fo Ks = —2 f, Hy +2 f, Hy —3 Hy hy + Hy hy—4 Ty kg, 


1 AY va Ta - 

—p, & + fo Gi — 4606, =—5 Tog; 

1 4 t / / vA 4 / 7 
—p, Us + fot, —4f, o, = — 27,6, + 3f,6,—4 Hj, g, + H, g,—5 Ty hy, 


1 ZA 7 f La 9 / 7 
—p; Ka + fo Ka —4 fo Ky = —3 9,6, + 28,6,—5 Th hy , 


1 wt io 7 7 7 iy / 1 
pp La + fobs —4 fo Ly = —2f, Hy + 3 f, Hy — 3 g, Hy — 36, hy ee 


+ 2g, H, + 2h, G,—4 Hy hs + hy H,—5 Th], 
1 Mt v va 7 LA / , 

p, Ms + SoMa — 4 fo My = —2f, Ky + 3f, K; — 3h, H, + 2h, Hy 

4 Wk FR 5 Tm, , 


Dieses Gleichungssystem ist fiir die folgenden Randbedingungen zu lésen: 

1) An der Stelle Y = 0, X > 0 muB T(X, Y) verschwinden. Das ist méglich, wenn alle Koef- 
fizientenfunktionen an der Stelle 7 = 0 verschwinden. 

2) Ander Stelle X = Oist T(X, Y) = 1. Diese Bedingung ist erfiillt, wenn fiir 4 —>0o T,(y) =1 
wird und alle iibrigen Koeffizientenfunktionen null werden. 


Gleichung (14) wurde erstmals von E. Pohlhausen! angegeben. Sie ergibt sich bei der Berech- 
nung der Temperaturgrenzschicht einer langsangestrémten ebenen Platte. 


1 E. Pohlhausen, Z. angew. Math. Mech. 1 (1921), S. 115. 


XXIX. Band 1960 K. Stephan: Warmeiibergang und Druckabfall laminarer Einlaufstromungen 181 


Wie aus dem System der Gleichungen (14) u. f. hervorgeht, gehért zu dem Strémungsfeld eine 
einparametrige Schar von Temperaturverteilungen mit der Prandtl-Zahl als Parameter. Um das 
Temperaturfeld iiber einen méglichst weiten Bereich der Prandtl-Zahl erfassen zu kénnen, wurden 
die Gleichungen (14) bis (17) fir 

ee == O51; ks 10 -100% 1000 


gelést. Fiir kleinere Werte als Pr = 0,1, wie sie zum Beispiel bei fliissigen Metallen vorkommen, 
darf die Warmeleitung in Strémungsrichtung im allgemeinen nicht mehr vernachlassigt werden. 
Diese ist jedoch in den Gleichungen (14) u. f. nicht beriicksichtigt. Es hat daher wenig Sinn, das 
Gleichungssystem auch noch fiir Werte Pr < 0,1 zu lésen. 

Fir Pr =0,1; 1; 10 wurden die Lésungen nach dem Runge-Kutta-Verfahren ermittelt. Da 
dieses sehr miihsam ist, wurde fiir Pr = 100 und Pr — 1000 ein Naherungsverfahren benutzt, das 
sehr gute Ergebnisse liefert. Man geht dazu von der Tatsache aus, daB fiir groBe Werte Pr die 
Temperaturgrenzschicht sehr diinn gegeniiber der Strémungsgrenzschicht ist. Die Temperatur 
T(X, Y) = 1 wird daher schon fiir sehr kleine Werte 7 erreicht, die innerhalb der Strémungs- 
grenzschicht liegen. Es ist daher naheliegend in den Gleichungen (14) u. f. die Koeffizientenfunk- 
tionen des Geschwindigkeitsfeldes durch die ersten Glieder ihrer Taylor-Entwicklung, Entwick- 
lungszentrum 7, = 0, zu ersetzen. Dann kann man die Gleichungen (14) u. f. durch Reihenansatze 
lésen. Auf die Wiedergabe der umfangreichen Zahlentafeln sei an dieser Stelle verzichtet. Sie sind 
in der Dissertation des Verfassers! angegeben. 

Zu einer Abschatzung des Konvergenzintervalles 0< X < X, der Temperaturverteilung 
I(X, Y) und der Geschwindigkeitsverteilung U(X, Y) gelangt man, wenn man zu U(X, Y) und 
zu T(X, Y) Oberreihen bildet und deren Konvergenzradius abschatzt. Eine Oberreihe zu U(X, Y) 
ist bekanntlich die Reihe U;(X). Eine Oberreihe zu T(X, Y) erhalt man, wenn man in dem Ansatz 
(12) die Koeffizientenfunktionen T;(7) durch ihre absoluten GréBtwerte ersetzt, die sich aus der 
numerischen Rechnung ergeben. Eine Abschatzung des Konvergenzintervalles dieser Oberreihen 
ergibt den Konvergenzradius 


X S5- 10 (19) 
mit 
v x x 
4s a UD ERD Rep 
v 


worin D der hydrodynamisch gleichwertige Durchmesser ist, der beim Rohr 2 s und beim ebenen 
Spalt 4s betragt. Der Wert (19) ist fiir das Strémungs- und Temperaturfeld des Rohres und des 
ebenen Spaltes giiltig. An der Stelle X, = 5-10 betragt die Geschwindigkeit der Potential- 
strémung etwa 1,15 - Up. 


6. Asymptotische Lésung fiir das Geschwindigkeitsfeld. Die bisherige Losung konvergiert nur 
fiir kleine Werte der Lauflange X < X,. Zur Berechnung der Geschwindigkeits- und der Tempera- 
turverteilung fiir groBe Werte X ‘wurde eine Entwicklung von der ausgebildeten Poiseuilleschen 
Parabel und von der konstanten Endtemperatur aus angesetzt. Ausgehend von einer 4lteren 
Arbeit von J. Boussinesq? fiihrte bereits H. Schlichting’ eine solche Rechnung fiir das Geschwindig- 
keitsfeld der Strémung im ebenen Spalt aus. B. Punnis* berechnete nach den gleichen Methoden 
das Geschwindigkeitsfeld der Rohrstrémung. Er gibt den Ausdruck 


U(X, R) =2[(— RY) + eA RR) + oF eo PAXRAR) | 
4 che 3X RR) + cy AX RY(R) +..-] (20) 


an mit 


c, = 2,242, tc =—0,595 und 4, = 16,0397, A, = 47,844, 


der in der vorliegenden Arbeit itbernommen wurde. In Gleichung (20) ist R = 1 — Y die dimen- 
sionslose Koordinate senkrecht zur Strémungsrichtung, die von der Rohrachse aus gerechnet wird. 


Die Funktionen R(R) (n = 1, 2, 3) und R¥’(R) wurden von B. Punnis* vertafelt. 


1 siehe FuBnote 5 von S.179. 

2 J. Boussinesq, Journ. Math. pur. et appl. 13 (1868), S. 377. 
3 H, Schlichting, Z. angew. Math. Mech. 14 (1934), S. 368. 

4 siehe FuBnote 4, von S 179. 
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Fiir die vorliegende Arbeit erwies es sich als zweckmafig, die Rechnungen von H. Schlichting! 


zu erweitern. Es sei jedoch nur das Ergebnis der Rechnungen mitgeteilt. Die Geschwindigkeits- 
verteilung im ebenen Spalt wurde ermittelt zu 


U(X, R) = = (1 — R2) + ¢, e74*S,’(R) + c2 e24XS,(R) +... (21) 


mit?¢— — 0,2881, c. = 0,083 und yy = 18,7441. Hierin ist R = 1— Y wiederum dimensions- 
lose Koordinate senkrecht zur Strémungsrichtung, die von der Symmetrieebene des ebenen Spaltes 
aus gerechnet wird. Die Funktionen S,/(R) und S,/(R) sind in der Dissertation der Verfassers? ver- 
tafelt. 


7. Asymptotische Lésung fiir das Temperaturfeld. Zur Ermittlung des Temperaturfeldes der 
Rohrstrémung mu8 man fiir groBe Werte X von der rotationssymmetrischen Form der Fourier- 
schen Energiegleichung ausgehen, weil fiir groBe Werte X die. Grenzschichtdicke nicht mehr ver- 
nachlassigbar klein ist gegenitiber dem Kriimmungsradius der Rohrwand. Die Fouriersche Energie- 
gleichung lautet dann mit den in der vorliegenden Arbeit benutzten dimensionslosen GréBen 


oT oT Ih GE held 
axe igs (aR TR eR) ia 
Fir die Temperaturverteilung macht man den Ansatz 
T(X, BR): => see es (23) 
n=0 


worin die yz, noch naher zu bestimmende Konstanten sind. Den Ansatz (23) sowie Gleichung (20) 
setzt man in die Fouriersche Energiegleichung (22) ein und nimmt die entsprechenden Differen- 
tiationen vor. Durch Koeffizientenvergleich erhalt man fiir die Glieder, welche den Faktor e “* 
enthalten, den Ausdruck 


36 + 28) +2 uy Pr(1— RF’) 8, =0, (24) 


der fiir die Randbedingungen #, (R = + 1) = 0 zu lésen ist. Gleichung (24) tritt bei der Berech- 
nung einer thermischen Kinlaufstrémung auf, deren Geschwindigkeitsprofil eine Poiseuillesche 
Parabel ist. Sie wurde zuerst von L. Graetz? und spater von W. Nufelt* gelést. Als Lésung ergibt 
sich eine Eigenfunktion #,(6;, R) mit 


6? =2 Pruy;=2Prn, (i =0,1,2...). 


Die Eigenwerte 8; wurden auBer von L. Graetz? und W. Nufelt® noch von einigen anderen Autoren 
ermittelt 5 & 7, 
Sehr genaue Rechnungen stammen von J. Schenk und J. M. Dumoré’. Sie errechneten die 
Werte 
Bo = 2.70437), B, = 6,6790 , 
pa LOG SS = Bo = 146TH. 


Fir groBe Eigenwerte f; gilt die Naherungsformel von H. A. Lauwerier® sowie von J. Sellars, 


M. Tribus und J. Klein® 
B; = 4 u a. z, ° 


Kin geschlossener Ausdruck fiir die Eigenfunktionen #,(f,, R) la8t sich nicht bestimmen. Man kann 
sie jedoch in Form einer Potenzreihe angeben. Paft man diese den Randbedingungen an, so ergibt 
sich 


Bo(B:. R) = B; oe (B;. R) : 


1 siehe FuBnote 3 von S. 181. 
2 siehe FuBnote 5 von S. 179. 
3 L. Graetz, Ann. Phys. 18 (1883) S. 79 u. 25 (1885), S. 337. 
page ela 54 (1910) S. 1154. 
. Sellars. M, Tribus u. J. Klein, WADC Technical Report 54—225, Michi iv. 
6 Th. M. Drew, Amer. Inst. Chem. Eng. 26 (1931), S. 26. a 
7 J. Schenk u. J. M. Dumore, Appl. Sci. Res. 4 (1953), S. 39.. 
8 H. A, Lauwerier, Appl. Sci. Res. A 3 (1953), S. 58. 
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Die noch freien Konstanten B; sind so zu ermitteln, daB T(X = 0, R) = 1 erfiillt wird. Das ist 
gleichbedeutend mit 


3 9,8) = 
Zur Berechnung der Glieder #,(R), 3,(R), . . geht : man so vor: In der Fourierschen Energie- 
gleichung (22) in welche die Ansatze (20) und (23) eingesetzt wurden, tritt die Exponentenfunktion 
e—*X mit folgenden Exponenten k auf: 


Ho> ] 
fas Myo +A | 
Me> My +A, Mo + 2A, \ (25) 


M3 > fe +A, fy + 2A,, fo + 3A,, 


An Stelle von j, sind dabei die Werte /, (n = 1, 2,...) und an Stelle von py die Werte Li, 
(@ = 0, 1,2...) zu setzen. ; 
Uber die Bepodeneen (25) verfiigt man nun in der Weise, daB die in jeder Zeile stehenden Werte 
einander gleich gesetzt werden. Durch diese Festsetzung ist es méglich, die Fouriersche Energie- 
gleichung in ein System von unendlich vielen gewéhnlichen Differentialgleichungen zu verwandeln, 
aus dem sich die noch unbekannten Funktionen #,(R), 3,(R), ... berechnen lassen. Durch Koef- 
fizientenvergleich ergeben sich die Gleichungen ; 


Of +H +2 uy Pr(L—R)9,=—2 po PR +24 PrERH, (26) 
ET el ae ae ao a, 


—2 pp Pr Ri 8) + 4A, Pr— 7 Ra 8 (27) 


Die Randbedingungen sind 
Bye Cites tea by 0 


2) 39, (R) = 1. 


te CSeETe (26) und (27) ist R, =c, R, und R, = c? 2 R,. Schreibt man 0, = c, B; 0, und 


8, = c2 B;9,, s0 kann man die Konstanten c, und B; kiirzen. Die Funktionen 8, und #, lassen 
sich fnrcls Potenzreihenansatze ermitteln. Die noch freien Konstanten B; Losi man so, dah 
das Temperaturprofil T(X, Y) der Entwicklung von vorn im Sinne des Kleinsten Fehlerquadrates 
moglichst gut mit dem Temperaturprofil T(X, R) der Entwicklung von hinten ibereinstimmt. Es 


soll also an einer Stelle X = Gs zum Beispiel X, = =5-10~—* nach Gleichung (19), gelten 


oe i [1(X,, R) — T(X,, Y)]’ aR =0, (28) 


1 
iy, | (7B 2) — 11%, YI aR =o, 9) 


Mit den hieraus errechneten Werten Bp, B,, ... laBt sich die Temperaturverteilung anschreiben. 
Man erhalt fiir die Temperaturverteilung Ausdriicke von der Form 


1(X, R, Pr) = > C, (Pr) e—*»P9* ©, (R, Pr). (30) 
n=0 


- Die Konstanten C,(Pr), k,(Pr) sowie die Funktionen 0,(R, Pr) sind in der Dissertation des Ver- 


fassers! vertafelt. 


1 siehe FuBnote 5 von S.179. 
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Zur Berechnung der Temperaturverteilung im ebenen Spalt geht man genau so vor wie bei der 
Rohrstrémung. Man hat lediglich die Fouriersche Energiegleichung fiir die ebene Strémung (3) 
zu benutzen. Das Ergebnis laft sich wiederum in Form der Gleichung (30) angeben, deren Werte 
C,(Pr), k,(Pr) und @,(R, Pr) ebenfalls in der Dissertation des Verfassers! vertafelt sind. 

Mit diesen Ergebnissen kénnen das Geschwindigkeits- und das Temperaturfeld graphisch dar- 
gestellt werden. Abb. 1 zeigt das Geschwindigkeitsfeld des ebenen Spaltes, Abb. 2 das Temperatur- 
feld fir Pr = 1. Zum Vergleich sind in Abb. 2 noch einige Temperaturprofile fiir Pr 1 ein- 


gezeichnet. 


if 


10° 


36:0" 81-107 6:10" 


Abb. 1. Geschwindigkeitsverteilung im ebenen Spalt. 


36-00% 87-192 


Abb. 2, Temperaturverteilung im ebenen Spalt. Pr = 1 (0,1; 10). 


Damit ist die Berechnung des Geschwindigkeits- und die des Temperaturfeldes abgeschlossen. 
Als wichtiges Resultat ergibt sich: Ganz vorn im Einlauf werden die Geschwindigkeits- und die 
Temperaturverteilung hinreichend genau durch die Gleichungen der langsangestrémten ebenen 
Platte geschrieben. In groBer Entfernung vom Eintrittsquerschnitt, wo die Poiseuillesche Parabel 
praktisch ausgebildet ist, gehorcht die Temperaturverteilung der Graetz-Nufelt-Lésung. 


8. Berechnung des Druckabfalles. Die Kenntnis der Geschwindigkeitsverteilung gestattet die 
Berechnung des Druckabfalles. Solange man von einer Potentialstrémung sprechen kann, ergibt 
sich der Druckabfall aus der Bernoullischen Gleichung 


Bae de Sey BO, 2 ease (31) 


mit Ap = p(X = L)— p(X =0). 
Aus der Beziehung (6) fiir die Potentialstrémung ergibt sich 


i 


al ak by S ae. (32) 
ovo iS 
2 
mit X = 3 ti DD: Dabei ist D der hydrodynamisch gleichwertige Durchmesser. Es ergaben sich 


v 
folgende Werte a;: 


— 24,379 + 2043,5 — 8192 
= 59,299 + 2565,1 — 136515 


ebener Spalt -+ 13,766 
Rohr + 13,766 


Kine Abschatzung des Konvergenzradius der Potenzreihe (32) ergibt: 
Xi-G 52 10ct 
1's, FuBnote 5 von S. 179. 
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Fir groBe Werte der Lauflange geht das Geschwindigkeitsprofil asymptotisch in die Poiseuillesche 
Parabel tiber. Dementsprechend gehen auch die Werte fiir den Druckabfall in die Asymptote 


= aX+b (33) 


iiber mit der Konstanten a und der Konstanten b, die den zusatzlichen Druckabfall einer hydro- 
dynamisch nicht ausgebildeten Strémung gegeniiber der ausgebildeten Strémung angibt. 

Die numerische Rechnung ergibt fiir den ebenen Spalt a = 96, b = 0,64 und fiir das Rohr 
a = 64, b = 0,87. Den Wert b = 0,87 berechnete B. Punnis 1. 


9. Berechnung der Warmeiibergangszahlen. Um einen Vergleich mit anderen Warmeibergangs- 
problemen zu ermdglichen ist es iiblich, rein formal die sogenannte Nufelt-Zahl einzufiihren, die 
wir fiir die Prandtl-Zahlen 0,1; 1; 10; 100; 1000 berechnen kénnen. Fiir alle iibrigen Prandtl- 
Zahlen in dem Intervall 0,1< Pr < 1000 erhalt man die Nufelt-Zahlen durch Interpolation. 

Die mittlere Warmeiibergangszahl a(x), die ein MaB ist fiir die bis zu einer Stelle x = L iiber- 
tragene Warmemenge werde definiert durch 


(x) ty(%) f= Up Fo ey [ty — tm(x)] « 


Hierin ist i,,(x) die mittlere logarithmische Temperaturdifferenz 


bk x) at (t — tw) — (tm — tw) 


In 0 te 
pe 
Af say 

eens - al [ igs 
HE F 
i | 
10 1 wooo, | | |i ees || 
5 7 100 : 
is | [ 
: 10 | 
19° : - 
: Pr=01 ; T Tt la | 

fine lias | 
3 | | ie 
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Abb. 3. Mittlere Nufelt-Zahl 


__ mit der mittleren Flissigkeitstemperatur 


Ne 


ie 0 eer aac | t(x,r) u(x, r) dF(r) . (34) 
r=0 
Hier ist F die Flache des Kanalquerschnittes und V das in der Zeiteinheit durch den Querschnitt 
strémende Volum. Ferner ist f die zwischen x = 0 und x = L warmeibertragende Wandflache. 


x 1 Siehe FufSnote 4 von 8.179. . 
a . ; 13 
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Mit den bereits in Abschnitt 3 definierten dimensionslosen Groen ergibt sich aus Gleichung (33) 
die mittlere Nufelt-Zahl 


ma(%) = “97? = — Fw B(). a 


Da die Geschwindigkeits- und die Temperaturverteilungen bekannt sind, kann man ¢,,(x) nach 


Gleichung (34) und damit auch Tok ) und Nu(X) berechnen. Man gelangt zu folgendem Ergebnis: 
Die Nufelt-Zahl nimmt am Eintrittsquerschnitt unendlich groBe Werte an. Mit zunehmender 


Lauflange X nimmt sie ab und nahert sich fir X —0o einem Grenzwert Nuco- Dieser ist unab- 
hangig von der Prandtl-Zahl und betragt fiir die Strémung im ebenen Spalt Nu = 7,55 und fiir 


die Rohrstrémung Nu, = 3,66. Beide Werte sind von mehreren Autoren?? gut bestatigt. 
Abb. 3 zeigt den beschriebenen charakteristischen Verlauf fiir die Nufelt-Zahl. 


10. Zusammenfassung. In der vorliegenden Arbeit wurden Warmeiibergang und Druckabfall 
nichtausgebildeter Laminarstrémung in Rohren und ebenen Kandlen berechnet. Hierzu wurden 
exakte Methoden der Grenzschichttheorie, wie sie zur Berechnung der Geschwindigkeitsverteilung 
iiblich sind, auf die Berechnung der Temperaturverteilung tibertragen. Diese Methoden bestehen 
darin, daB man eine ,,Entwicklung von vorn“ ansetzt, die vom Eintrittsquerschnitt an giiltig ist, 
und eine ,,Entwicklung von hinten“, die in gréBerer Entfernung vom Eintrittsquerschnitt an gilt. 
In dem Bereich, wo weder die Entwicklung von vorn noch die Entwicklung von hinten giiltig sind, 
weil nicht geniigend Glieder zur Verfiigung stehen, erzielt man durch Interpolation brauchbare 
Ergebnisse. Die in der Arbeit beschriebenen Methoden sind auch geeignet, den Warmeiibergang 
und Druckabfall in Ringspalten zu berechnen. Das wird jedoch eine erhebliche Rechenarbeit 
erfordern, weil als weiterer Parameter noch das Verhaltnis der Radien des inneren und des auBeren 
Rohres auftritt. 


(Eingegangen am 24. Juli 1959) 


Anschrift des Verfassers: Dr.-Ing. K. Stephan, Karlsruhe, Lehrstuhl fiir Thermodynamik der Technischen 
Hochschule. 


1 H. Hahnemann, Warme- und Kaltetechn. 44 (1942 S. 167; F - 
2 K, Elser, Schweiz. Bauz. 69 (1951), S. 641. ee pe pe ee), 8 2 


3 E. Eckert, Introduction to the transfer of heat and mass, S. 100, New York 1950. 
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; Integralgleichungen gerader Kipptrager 
mit dimnwandigen, offenen und quasi-wélbfreien Profilen 


Von L. Kraus 


1. Einleitung. Fiir die von Chwalla! aufgestellte komplizierte Differentialgleichung fiir das 
Kippen des geraden Tragers mit diinnwandigem, einfachsymmetrischem offenem Profil liegen 
Lésungen in verschiedener Form vor. Meissner? berechnet mit dem Energieverfahren die Kipp- 
zahlen von Tragern auf zwei Stiitzen unter Gleichstreckenlast, mit einem Endmoment und mit zwei 
gegengleichen Endmomenten. In einer umfassenden Abhandlung liefert Pettersson? Lésungen 
der Kippdifferentialgleichungen mit Hilfe von Potenzreihen bei linearem Verlauf der Kippmomente; 
er behandelt den Trager auf zwei Stiitzen mit einer Einzellast in Tragermitte ohne und mit gegen- 
gleichen Endmomenten, ferner den Trager mit einem Einzelmoment am Tragerende, und symmetri- 
sche Durchlauftrager auf vier Stiitzen mit einer Einzellast in der Mitte des Mittelfeldes. AuBer der 
Energiemethode und der Lésung mit einem Potenzreihenansatz lat sich auch das Differenzen- 
verfahren verwenden. Hier soll der Versuch unternommen werden, das Kipp-Problem des einfach- 
symmetrischen Querschnittes mit Hilfe der fiir die numerische Rechnung bequemen Integral- 
gleichungen zu lésen. Fir Trager mit einfach-symmetrischen Profilen gelingt es jedoch nur bei be- 
bestimmten Querschnittsformen und Belastungen, die Integralgleichungen mit Hilfe von bekannten 
Funktionen auszudriicken. Es ist dies der Trager mit quasi-wélbfreiem Querschnitt (T-, Winkel- 
oder Kreuzprofil mit Wélbwiderstand Cy, = 0), wenn zugleich das Biegemoment einen linearen 
oder quadratischen Verlauf hat. 

Im folgenden werden die Integralgleichungen fiir solche Trager angegeben und fiir ein bestimm- 
tes T-Profil numerisch ausgewertet. Die berechneten Kippzahlen werden in der Gestalt von Dia- 
grammen und Tabellen mitgeteilt. 


2. Kippdifferentialgleichung. Die Kippdifferentialgleichung fiir einfachsymmetrische Profile4 
lautet in den dimensionslosen Koordinaten # und s: 


— <M gv + {[1 4+ O81 M(s)] OY +2 | Mr) + of P/BD a5|9 ao. (1) 


Es bedeuten Cy, den auf den Schubmittelpunkt bezogenen Wélbwiderstand [kpcm*], B die Biege- 
steifigkeit fiir die Minimumstragheitsachse [kpcm?], D die de St. Venantsche Drillsteifigkeit [kpcm?], 
Il die Tragerlange [cm], # den Drillwinkel und z die Langenkoordinate [cm], deren Achse die Trager- 
schwerachse ist. Weiterhin ist s = z/] die dimensionslose auf die Tragerlange 


bezogene Koordinate, h die Tragerhéhe [cm], v =@h/2 der Abstand [cm] der y 
Streckenlast p [kp/cem] vom Schubmittelpunkt, y,, der Abstand [cm] des PIs 
Schubmittelpunktes vom Schwerpunkt und y = fyrdF/J, eime dem Ale 

: F = 
Querschnitt eigentiimliche Lange (Abb. 1). Es sind ® = 2 (2 yy,—y)/h und s : 


B =h/B/D/ 21 dimensionslose QuerschnittskenngréBen. SchlieBlich stellt 
M, ein beliebiges Vergleichsbiegemoment [kp cm], M(s) den qualitativen, 


dimensionslosen Verlauf des Momentes und A = M, i/ /BD den dimensions- 44), 1. Das einfach-symme- 
losen Eigenwert dar. Striche bedeuten Ableitungen nach s. Seitliche Biege- Bee broth Coane 
momente, deren Vektoren mit der Achse fiir das kleinere Tragheitsmoment 
zusammenfallen, sind in der Differentialgleichung weggelassen. 
Um von der Differentialgleichung (1) zu Integralgleichungen zu gelangen, schlagen wir den- 
selben Weg ein wie friiher®. Doch zeigt es sich, daB die dabei auftretende homogene Gleichung 
cs , 
C 7U7 
Liu) =—peut + {1+ 0p4M()] v} =0 (2) 
1 FE. Chwalla, Sitz. Ber. Akad. Wiss. Wien, mathem.-naturw. Kl. Abt. Ila, 153 (1944), S. 25. 
2 F. Meissner, Diss. d. Deutschen Technischen Hochschule Brinn 1944 und Stahlbau 28 (1955). 
8 Q, Pettersson, Bull. Nr. 10 from The Div. of Building Statics and Struct. Eng. at the Royal Inst. of Techn., 
Stockholm. 
4 siehe FuBnote 1 und 3. 
5 L. Kraus, Ing.-Arch. 26 (1958), S. 1. 
13) 
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nur in besonderen Fallen geschlossen lésbar ist. Von der Lésung durch Potenzreihen soll abgesehen 
werden; denn dann ist es schon zweckmafiger, gleich die vollstandige Differentialgleichung (1) so 
zu behandeln. Setzen wir jedoch Cy, = 0 und nehmen wir fiir M(s) einen linearen bzw. quadra- 
tischen Verlauf an, so driicken sich die Lésungen von (2) 
durch In- bzw. Yt Tg- und arctg-Funktionen aus. Beziig- 
lich des Biegemomentes hat man damit die baustatisch 
interessanten Fille eingefangen. Was die Querschnitte an- 
betrifft, schlagen wir vor, sich einer zeichnerischen Nahe- 
rung zu bedienen. Mit Hilfe der im folgenden entwickelten 
Integralgleichungen lassen sich die Kippzahlen fiir T-Profile 
berechnen. Fiir den entsprechenden doppelt-symmetrischen 
I-Querschnitt sind sie entweder bekannt, oder sie kénnen 
auf dieselbe Art ermittelt werden. Tragt man diese Kipp- 
werte iiber dem Verhialtnis der Flanschbreiten auf und 
Abb. 2. Die Kippwerte iiber dem Verhiiltnis verbindet man die Punkte durch eine Kurve (Abb. 2), so 

Ra gg he erhalt man Naherungswerte fiir die sich dazwischen ein- 

ordnenden Profile. 


Die Integralgleichungen fiir das T-Profil lassen sich auch auf mechanischem Wege herleiten. 


Hierfiir betrachte man einen Ersatztrager mit der verdnderlichen Drillsteifigkeit D(s) = 
Dy {1 + « M(s)] unter der Wirkung eines Einzeltorsionsmomentes M =1 an der Stelle z =1t und 
berechne die Lésung fiir Is, t) in den beiden Bereichen 0 Xs Xt undt <s <1. Die zur expli- 
ziten Lésung von @ nétigen vier Bedingungen liefern die beiden dem Kipptrager zugeordneten 
Randbedingungen und zwei Ubergangshedingungen bei s = t (Gleichheit von # und eine Sprung- 
bedingung, die iiber die Gré®en der durch M in beiden Bereichen erzeugten Torsionsmomente Auf- 
schluB gibt). Die Randbedingungen lassen sich als Bedingungen fiir eine starre oder elastische Tor- 
sionseinspannung oder fiir ein freies Ende auffassen. Den mit j? und # behafteten Ausdruck in (1) 
denke man sich als eine tiber den Trager verdnderliche Torsionsstreckenbelastung m(s). Der 
Drillwinkel 3, den m(t) an der Stelle s erzeugt, hat die Form 


i 
Hs) = f Ws, t) m(t) dt. 
i 
Damit ist das Problem in der Gestalt einer Integralgleichung gegeben. 


3. Integralgleichungen bei linearem Verlauf des Biegemomentes. Wir setzen Cy, = p = 0 und 
M(z) = Mjs. Es ist M, das Biegemoment am Ende des Integrationsbereiches. Damit ergibt sich 
mit x = OBA aus (2) 


[1 +s) #) =—2?SO. (3) 
Nach dem Verfahren von Lagrange erhalten wir daraus die Integralgleichung 
Ae . 1+ xt 
00) =F [ae] P 00 @ +n tx) G+. (4) 
0 
Thre Ableitung nach s lautet 
4 : t? H(t) dt C,x 
aye oor 
3s) ee (5) 
: 0 


Mit (4) und (5) werden wir einige Lastfalle behandeln. 


a) Der Kragtrager. Es sei P die Kipplast am freien Ende des Tragers und v = 9 h/2 ihr Ab- 
stand vom Schubmittelpunkt. Da P im Flansch Zugspannungen erzeugt!, muB das Biegemoment 
das negative Vorzeichen erhalten. In (3), (4) und (5) ist also x durch —x zu ersetzen. 


1 siehe FuBnote 1 von S. 187. 
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Mit den Randbedingungen — #(0) =1 8 0 (0) und (1) = 0 erhalten wir dann nach einer 
kurzen Rechnung 


s 


Hs) = ae {/( + In(1—x?) 
x Lt gma—») Ff 


l—xs 
In i 


= a 


| t? D (t) dt 


== 


1 
+f t+gma—xs) 


naa (6) 


Bei umgedrehtem Profil wechselt 4 in (6) das Vorzeichen. Die Zuordnung der Vorzeichen ist in 


Abb. 3 dargestellt. 


éP 
7 System 
TN, M, y » 
2 Pp WINS WaT 
’ - Orillwinkelverlauf 
eriau. 
A>0  -A<O Be 
Abb. 3. Die Zuordnung der Vorzeichen Abb. 4. Systemskizze fiir den Trager Abb. 5. Der Trager auf zwei Stiitzen mit einer 
von 4 beim Kragtrager mit einer Ein- auf zwei Stiitzen mit Endmomenten. in der Mitte wirkenden Einzellast. 


zellast am freien Tragerende. 


b) Der Trager auf zwei Stiitzen mit Endmomenten. Das Biegemoment M lautet mit 
den Abkiirzungen 
_ MI 
M, — BD 


gemaB der Abb. 4 
M = M, [a + (l— a) s] = M, M(s). 


Statt (3) haben wir 


x (1—2) Bae # a : 
([) +7cees |} mec x) s] 0. (3a) 
Wenn wir x (1 — a)/(1 + a x) =x und 22/(1 + « x) = (A)? setzen, kénnen wir wieder (4) und (5) 
verwenden. Die Randbedingungen lauten (0) = 0(1) = 0. 

Daraus ergibt sich die Integralgleichung 


is) =Fia fl +3] IFas 


0 


i. {fase In 1 +* | [x + (1—a) P 0(0) dt 


1+% 
ae 


+ [in (1 + #s) |Inoe, fe +0—9er aad () 


Auch hier wechselt das Vorzeichen von A, wenn das Profil umgedreht wird. 


c) Der Trager auf zwei Stiitzen mit Einzellast. Gema® der Abb. 5 ist M(s) = s und 
LS Pr//BD . Die Lésung dieses Lastfalles liegt in geschlossener Form vor!. Um die Abwei- 
chungen von den strengen Werten zu ermitteln, wird dieser Fall mit den Methoden der Integral- 
gleichungen behandelt. Als Randbedingungen werden vorgeschrieben (0) = 0 und (1 —) X 
& (1) =-(« o/®) 0(1). Aus der letzten Gleichung geht hervor, da8 — ahnlich wie beim Rechteck- 
querschnitt — der Drillwinkel in Tragermitte keine horizontale Tangente hat. Fir @ = 0 ist die 
zweite Randbedingung mit der entsprechenden Randbedingung des Rechteck-Tragers identisch. 


1 siehe FuBnote 3 von S. 187. 
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Aus diesen Randbedingungen, sowie aus (4) und (5), bekommen wir die folgende Integral- 


gleichung: 
s 


Z (i +s) 
9) = spa —»—00 +H nd $4) i le (Lee Eres 


X [In (1 +2 0)]2.0(0) de + [|oo—a+e (nes | ae) 


-[In (1 +% s)]? A(t) a} . 1 (8) 


Das Vorzeichen von A und x wird negativ, falls P im Flansch Zugspannungen erzeugt. 


4, Integralgleichungen bei quadratischem Verlauf des Biegemomentes. Wiederum sei Cy = 0. 
Wir beschrinken uns auf rein parabolische Biegemomente M(s) = s? bzw. M(s) =s (1—s). 
Die Integralgleichung baut sich je nach der Wirkungsrichtung der auBeren Krafte aus {ry Tg- oder 
arctg-Ausdriicken auf. Auch wenn zusiatzlich Biegemomente an den Tragerenden wirken, 1aBt sich 
die Integralgleichung mit Hilfe der Partialbruchzerlegung finden. Wir geben zuerst die Integral- 
gleichungen fiir M(s) = s? an. Der Fall M(s) = s (1—s) baw. M(s) = 1— s? laBt sich aus dem 
ersten Fall ohne Schwierigkeit entwickeln. Fiir negative A und x [in (3b), (4b) und (5a) sind A und 
% dagegen positiv zu zaihlen|] ist also 


[l—xs)#”) =—2? *+obp/BD mp 0b. (3b) 


Daraus leitet sich mit der Abkiirzung m(s) = s+ + 0 B p /B D 1/M? die folgende Integralgleichung 
her: 


OR - | [Qe Tq (Vx t) — Me Tq (x s)] m(t) HC) dt + CMe Tq(Vxs)+ Cy. (4a) 


Differentiieren von (4a) nach s ergibt 


artes “ oe or 
eyes | TO) e+ G ES. (5a) 
Wechseln p, J und x das Vorzeichen, dann bekommen wir statt (3b) 


(1 +x) 0) =—}?|t— of pBD a5 0. (3c) 


Daraus erhalten wir mit der Abkiirzung n(s) = st —o B p /B D 1/M? die Integralgleichung 


s 


Hs) = le al [arctg (24 t) — arctg (Vx. s)] n(t) A(t) dt + C, arctg (Vx s) + (C,. (4b) 


Thre erste Ableitung nach s lautet 


s 


Fieetiee. n(t) Vx 
os) =— 2 eae aa+e  .. (5b) 
0 
a) Der Kragtrager. Wir behandeln den Lastfall Gleichstreckenlast. Es ist M (z) = — M,3? 


mit M, = pl?/2. Damit wird m(s) = st + 2 B/A und n(s) = st— 2 0 Ba. Randbedingungen sind 
#7 (0) =0 a B(1) = 0. Damit ergibt sich aus (4a) und (5a) [Zugspannungen im Flansch] die 
Integralgleichung 


Hs) = e | f [mo m(t) [Me Xa (Vx) — Wr Xa (x s)] H(t) de 


0 


= | m(t) [Ur Tg (Vx) — Wr Xo (y/x t)] a) a} (9) 


+ 
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Wirken im Flansch dagegen Druckspannungen, 


(4b) und (5b). Sie lautet 


Os) s2= Ve | | n(t) [arctg (x) — arctg (Vx s)] O(t) de 


0 


dann bekommen wir die Integralgleichung aus 


+ | n(t) [arctg (x) — aretg (Vx t)] Oe) de : 


(10) 


b) Der Trager auf zwei-Stitzen. Fir Gleichstreckenlast p = 2 q ergibt sich gemaB Abb. 6 
das von der Tragermitte aus gezihlte Biegemoment zu M(z) = M, (1 — s?) mit M, = ql?. Der 
Drillwinkel verschwindet an den Auflagern und ist zur 


Tragermitte symmetrisch. Bei positivem M,, A, x und Pp p=2g 

(Druckspannungen im Flansch) haben wir prenrsriisssr ane 
ee x ; = 22 ‘ 9 Be RENE ——- 

(2 i ta) 0 = ree — 3%)? aoa | O. (3d) 


Die Randbedingungen stimmen mit denen von Abschnitt 

4a iiberein. Die (3d) zugeordnete Integralgleichung ist 
mit (9) identisch, wenn dort 2/(1-+ x) statt 22 und 
z[(1 + x) statt x und m(t) = (1—#)? +20 B/A statt m(t) 2 Mae Bes Ae deere ae ie 
eingesetzt wird. 


/ 


Orillwinkelveriauf 


In der bei Rechenautomaten iiblichen algorithmischen Schreibweise lautet das oben Gesagte 
wie folgt ; 


Igl. far [(3d); 0’(0) = 91) = 0] = (9), | 
Ail + x) > 7, 


l (11) 
z#/(l+z)>x, ( 
—#)?+20pA> mit). J 


Bei Zugspannungen im Flansch wechselt das Vorzeichen von M,, A, x und p in (3d). Statt (3d) 
gilt 
% alg 2 2 
(1 + 38) songat | ee gears ee (36) 
In (3e) sind A und x positiv. Wir schreiben wieder 
fol dir [(.e)3. 4 (0j— 01) —0]-—.(10):, 
Pil=—x) >, 
x/l—x) >, 
(1—1)2— 20 8/, > mQ). 


(12) 


5. Numerische Behandlung der hergeleiteten Integralgleichungen. Die in den vorhergehenden 
Abschnitten aufgestellten Integralgleichungen werden fiir den Fall ® = 1,631176 numerisch aus- 
gewertet. Dieser D-Wert liegt auch den Untersuchungen von Meissner und Pettersson 
zugrunde!, so da bei zwei Lastfallen gewisse Vergleichsméglichkeiten bestehen. Die 

_Abmessungen des entsprechenden T-Profiles sind in Abb. 7 eingezeichnet. 

Zur numerischen Rechnung bringt man die Integralgleichungen durch Umwandlung 
der Integrale in Summen in die Form von homogenen Gleichungssystemen fiir # an den 
aquidistanten Stiitzstellen s;. Die Koeffizienten der zugehérigen Matrix enthalten aber 
den Eigenwert nicht in derselben einfachen Form wie die gewéhnlichen Matrizeneigen- 
wertprobleme. Zur Lésung dienen die beiden folgenden Verfahren. Fir gewisse ange- 


. : Abb. 7. Die 
nommene Groen J, werden die Werte der zu dem homogenen Gleichungssystem ge- —_Abmessungen 


hérenden Determinante AN mit dem Gaufschen Verfahren berechnet. Nach Wahl von dee ee reas. 


drei passenden AN-A,-Wertepaaren liefert eine quadratische Interpolation naherungs- 

weise den Nulldurchgang der Determinante und den zugehérigen /-Wert. Dieses Verfahren wird bei 
4-reihigen, symmetrischen und unsymmetrischen und bei 8-reihigen symmetrischen Matrizen ver- 
wendet. Im Falle unsymmetrischer 8-reihiger Matrizen wird das Verfahren der Matrizen-Eigen- 


1 siehe FuBnoten 2 bzw. 3 von S. 187. 
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wertbestimmung bei unsymmetrischer Koeffizientenmatrix in abgewandelter Form benutzt. Fiir 
Werte J,, von denen man annimmt, daf sie in der Nahe des wirklichen Eigenwertes liegen, ermittelt 
man die unsymmetrische Koeffizientenmatrix §t(A,). Damit ist ein spezielles Kigenwertproblem 
[R(A.) — E/A2] x = 0 gegeben, dessen kleinstes A, mit dem Matrizeniterationsverfahren! berechnet 
wird. Die aus drei gut liegenden /,-/,-Wertepaaren gebildete Parabel schneidet die Gerade 
dy =1q im gesuchten Wert J. Zu jedem A gehoren die zwei Parameter g und f. Fiir die Zahlen- 
rechnung werden gewahlt 


[ 0 (Last P oder p wirkt im Schubmittelpunkt), 
o = 4 — 0,625 (Last P oder p wirkt im Schwerpunkt), 
—2 (Last P oder p wirkt im freien Ende des Steges) 


und f? = 0,025; 0,050; 0,100; 0,250 und 0,500. % 

Die Rechenarbeit besorgt die von der Deutschen Forschungsgemeinschaft bereitgestellte 
elektronische Rechenanlage IBM 650. Sie fiihrt im Bell-Code die Berechnung, die Elimination und 
das Ausrechnen des Determinantenwertes einer 4-reihigen Matrix in 15 Sekunden und einer 8- 
reihigen symmetrischen in 45 Sekunden durch. Die Behandlung einer 8-reihigen unsymmetrischen 
Matrix verlangt durchschnittlich 4 bis 5 Minuten Rechenzeit. Hier wird lediglich die Integral- 
gleichung (11) mit diesem Verfahren gelést. Bei der nicht symmetrischen Integralgleichung (9) 
liegen die aus der 4-reihigen Matrix berechneten Naherungen so nahe an der Grenzhyperbel A = 1/ 
@ B, daB in den Fallen 9 + 0 die Berechnung der 8-reihigen Matrix unterlassen werden kann. Bei 
(9) und (11) ist eine Symmetrisierung wegen des durch die negativen 9- Werte bedingten Vorzeichen- 
wechsels von m(t) unméglich. Dieser ProzeB gelingt dagegen bei den mit n(t) behafteten Gleichun- 
gen (10) und (12). Es wird gesetzt g(s) = 0(s) /n(s) . Die Integralgleichungen (6), (7) und (8) werden 
nach Einfiihrung von y(s) = 0(s) M(s) symmetrisch. 

Die A, werden zuerst willkiirlich gewahlt. Bei linearem Momentenverlauf liegen die negativen A, 
die im Flansch Zugspannungen hervorrufen, stets unter der Grenzhyperbel 1 = 1/@f. Dasselbe 
gilt fiir die positiven A von (9) und (12). Weiterhin ist. meist 7 (6 = 0) bekannt. Damit hat man be- 
reits einen guten Anhalt. Mit den schnell rechnenden Programmen fiir die 4-reihige Matrix bestimmt 
man fiir verschiedene A, die A N-Werte und kommt bald zum Nulldurchgang. Die A,-Werte fiir die 
8-reihigen Matrizen werden dann in der Nahe davon angenommen. Fiir positive A, von (6), (7), 
(8), (10) und (11) bestehen keine Grenzkurven. Es mu8 solange probiert werden, bis AN das Vor- 
zeichen wechselt. Die Form der Eigenfunktion gibt dann Auskunft iiber die Nummer des Eigenwertes. 

Im folgenden Abschnitt 6 sind die berechneten Kippwerte nach Lastfallen in Tabellen zusammen- 
gefaBt, und ihr Verlauf ist in Diagrammen dargestellt. Die Abweichung der aus 4- und 8-reihigen 
Matrizen ermittelten 2 voneinander ist nur im Falle des Kragtragers mit Einzellast geringfiigig 
(siehe Abb. 8a und 8b und Tabelle la und 1b). Bei den anderen Lastfallen gehen die Kurven fiir 
bestimmte 0 teils wenig, teils aber ganz betrachtlich auseinander (siehe die entsprechenden Ta- 
bellen und Kurven des Abschnittes 6). Die Naherungen der 8-reihigen Matrizen liegen nicht immer 
tiber denen der 4-reihigen. Bis auf die Kurve 90 = — 0,625 von (9) [Abb. 11b] hat der Fehler der Ab- 
weichung der beiden Naherungen immer dasselbe Vorzeichen. 

_ Um fiir alle oben behandelten Lastfalle die in Abb. 2 angedeutete Naherungslinie zeichnen zu 
kénnen, miissen neben den bereits vorliegenden Kippwerten der anderen Lastfalle auch die des 
Zweistiitzentragers mit doppeltsymmetrischem I-Profil unter der Wirkung von zwei voneinander 
verschiedenen Endmomenten bekannt sein. Deshalb werden dafiir schlieBlich noch Naherungs- 
werte fiir 2 mit Hilfe von 4- und 8-reihigen Matrizen aus der folgenden Integralgleichung? 


1 
p(s) = 22 [ K(s, t) y(t) dt berechnet. Dabei gilt 

0 - 
__ Gin [a(1—s)] Ginat 


K(s, t) =| (Ves | Ms) M() fr 0 <i Zs, 


a Gina 
aa aA oe 
RG y= f (ge Gin [a . 2 eel my MG) fins = eg 
13 
1 MI M, i 


o= M(s) = [« + (l—a)s] 
und y(s) = #(s) M(s) . 


* Zurmihl, Matrizen, Berlin/Géttingen/Heidelberg 1958. 
® bag. deren Herleitung siehe FuBnote 5 von S. 187. 
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6. Diagramme und Tabellen. 


A-PL/YBD -A=P\VBD 


System: 


“4. o=0 
| Querschnitt ine - 9605 


ae Q=-2 


3 Q--2-——— (x2Lastangrifspunkt) 


Q=-0625 


° Wes 005 OF 05 05 ORs 006 af 025 Mato. aS 
a b 


Abb. 8a und b. Verlauf der aus (6) berechneten Kippwerte A. 


A aus 8-reihiger Matrix; ------------ 4 aus 4-reihiger Matrix. 


Tabelle la. Zusammenstellung der in Abb. 8a aufgetragenen Kippwerte 4. 


wane e = — 0,625 


8-reih. Matrix 


e=—2 


4-reih. Matrix 8-reih. Matrix 4-reih. Matrix 4-reih. Matrix 


8-reih. Matrix 


Tabelle 1b. Zusammenstellung der in Abb. 8b aufgetragenen Kippwerte i. 


e = — 0,625 
8-reih. Matrix 


e=—2 
8-reih. Matrix 


g@=0 


2 
B 4-reih. Matrix 8-reih. Matrix 4-reih. Matrix 4-reih. Matrix 


194 


L. Kraus: Integralgleichungen gerader Kipptrager Ingenieur-Archiv 


Abb. 9a und b. Verlauf der aus (7) berechneten Kippwerte A. 


—______— 4 aus 8-reihiger Matrix; 
ete ert tee A aus 4-reihiger Matrix. 


05 


20 


16 


=. 


G7 O25 


a 


o=1 


hex. — 


a[(x ® 6/2) + V1 + @% A/2)] 


4-reih. 8-reih. 
Matrix Matrix 


o=0,5 


A=Myl/ BD 


6005 


7504) 


5562 


th POCO 


Querschnitt: ik 


3He 


a 0005 005 G7 025 G5 


a— 0 


4-reih. 
Matrix 


8,280 


9,572 | 28,426 
11,486 | 35,508 
15,421 | 49,442 
20,135 | 65,130 


20) 


23,001 


8 -reih. 
Matrix 


4- und 8-reih, 
Matrix 


B hex, = 12/0 p 1B —1/0 1z—110 6 A= —1/9 p 12 —1/0 8 == Yo 6] A=— 1/0 6 
al(n@p/2)—YI+@Os) | Mawix | Matix | Macix | Matiz | Mariz | mente | Maries” 
0,025 — 2,117 — 2,684 — 3,877 
0,050 1831 — 2,275 9,749 
0,100 — 1,498 — 1,824 — 1,939 
0,250 — 1,081 — 1,226 — 1,226 
0,500 — 0,809 — 0,867 — 0,867 


HEP 
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_ p72 
ent’ ee” 


Yysten: es \ | 2p 
L—] ee] aA System: ee N 
A \t ee pS 
30}—+ A . 

uerschnitt: #88 = -0625 — -A=1]08 pene 

Q E = ie Ee ea ee : Querschnift: focus 
ae =- — See = 
(x SLastangriffspunkt) g Stee ae 0-0 
\ 


(x4 Lastangrifspunkt) 


20 


10 


ol 0 
0025 005 07 025 05 005 005 Of 005 05 


Abb. 10a und b. Verlauf der aus (8) berechneten Kippwerte A. 


4 aus 8-reihiger Matrix; ------------ A aus 4-reihiger Matrix;-++-+.-«.-+.- A aus Abh., v. Pettersson 
[s. FuBnote 3 von S, 187]. 


Tabelle 3a. Zusammenstellung der in Abb. 10a aufgetragenen Kippwerte i. 


= 05625 


@=0 
8-reih. Matrix 


4-reih. Matrix 8-reih.-Matrix 4-reih. Matrix 


Tabelle 3b. Zusammenstellung der in Abb. 10b aufgetragenen Kippwerte A. 


o = — 0,625 o=—2 


4-reih. Matrix 8-reih, Matrix 


=U 


8-reih. Matrix 4-reih. Matrix 8-reih. Matrix 


4-reih. Matrix 


Ingenieur-Archiv 


196 L. Kraus: Integralgleichungen gerader Kipptrager 


A=plf2 BD 


A-plj2yeo" 


55 =a] 

P E=a=_s=— 

System: thitite 
pate Z 30 _0=-2 
-. Querschnitt:  -9=-0625 
Querschnitt: BS ae “pe ~p=0 
\o=-2 Be (x4Lastangrifspunkt) ay 6 
4 (x4Castangriffspunkt) ——— Z 
Le o= 0 
7 


0 
o5025 005 a7 025 05 9005 005 G7 


Abb. 11a und b. Verlauf der aus (9) und (10) berechneten Kippwerte A. 


A aus 8-reihiger Matrix; ------------ A aus 4-reihiger Matrix. 


Tabelle 4a. Zusammenstellung der in Abb. 11a aufgetragenen Kippwerte A. 


oe = — 0,625 
4-reih. Matrix 


e=0 


8-reih. Matrix 4-reih. Matrix 8-reih. Matrix 


4-reih.Matrix 8-reih. Matrix 


Tabelle 4b. Zusammenstellung der in Abb. 11b aufgetragenen Kippwerte A. 


e=0 = — 0,625 
4-reih. Matrix 8-reih. Matrix 4-reih, Matrix 8-reih. Matrix 4-reih. Matrix 8-reih. Matrix 
8,464 
0,050 13,682 10,072 9,137 7,706 2,934 2,935 
0,100 17,411 12,567 9,329 8,471 2,179 2,212 
0,250 25,328 17,930 8,390 9,236 1,416 1,421 
0,900 34,649 23,977 6,870 9,025 1,012 1,013 
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tA=p1//B0 


6 
4 A=pl’ BD 
7=ép 
20 V System: eee AS 
\ hela 
A=19 ARE 
1710 Uatsd Querschnitt: ee 0605 
oe 
i (x4Lastongrif'spunkt) 
15 ; 


O5 


tee. 
Querschnitt: iss — 0625 
‘0= -2 
(x4 Lastangriffspunkt) 


Of 025 a5 Oo TR5 aE ar 025 a5 


b 


oes 005 


a 


Abb. 12a und b. Verlauf der aus (11) und (12) berechneten Kippwerte 4. 


A aus 8-reihiger Matrix; ------------ A aus 4-reihiger Matrix;-++----++-+0+ A aus Diss. Meissner 


[s. FuBnote 2 von S. 187]. 


Tabelle 5a. Zusammenstellung der in Abb. 12a aufgetragenen Kippwerte A. 


o = — 0,625 = 2 


8-reih. Matrix 


g=0 
4-reih. Matrix 8-reih. Matrix 


B 


4-reih. Matrix | 8-reih. Matrix 4-reih. Matrix 


Tabelle 5b. Zusammenstellung der in Abb. 126 aufgetragenen Kippwerte i. 


oe = — 0,625 @=—2 


4-reih. Matrix 8-reih. Matrix 


OEM 


4-reih. Matrix 8-reih. Matrix 


B? 


4-reih. Matrix 8-reih. Matrix 


e . - . r-Ar bhi 
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$A=M, 1/60 
ole Mors 
20 —— System:  <s0 ee 
—— nr 
Querschnitt: 
foe 
Jip a 
0005 005 = G25 OS 
Abb. 13. Verlauf der aus (13) berechneten Kippwerte d. 
A aus 8-reihiger Matrix; ------------ Aaus 4-reihiger Matrix. 
Tabelle 6. Zusammenstellung der in Abb. 13 aufgetragenen Kippwerte i. 
Pan | c—.0 a = —0,5 a=—l 
B 4-reih. 8-reih. 4-reih. 8-reih. 4-reih. 8-reih. 4-reih. 8-reih. 4-reih. 8-reih. 
Matrix Matrix Matrix Matrix Matrix Matrix Matrix Matrix Matrix Matrix 


7. Zusammenfassung. Fiir den geraden Trager mit einfach-symmetrischem, offenem und quasi- 
wélbfreiem Profil (der Wélbwiderstand ist Null) unter der Wirkung von linear und quadratisch 
verlaufendem Biegemoment werden die Kipp-Integralgleichungen fiir verschiedene Randbedin- 
gungen mitgeteilt. Fiir diese Integralgleichungen bestehen Programme, mit denen Kippwerte 
fir bestimmte Parameter (Festwerte des Querschnittes und des Lastangriffes) mit Hilfe des 
elektronischen Rechenautomaten IBM 650 berechnet werden. Die numerischen Ergebnisse sind 
in Form von Kurven und Tabellen zusammengestellt. Es gelingt nicht, die Integralgleichungen 
des einfach-symmetrischen I-Profiles mit bekannten Funktionen darzustellen. Fiir solche Profile 
wird ein zeichnerisches Naherungsverfahren vorgeschlagen. Dessen Genauigkeit hangt von den 
naherungsweise numerisch ermittelten drei Stiitzwerten und von der individuellen Annahme des 
mutmaflichen Kippwertverlaufes dazwischen ab. 

Mein Dank gilt den Herren Professoren Dr.-Ing. K. Marguerre und Dr. A. Walther fiir die freund- 


liche Férderung und der Deutschen Forschungsgemeinschaft fiir die finanzielle Unterstiitzung 
dieser Arbeit. 


Aus dem Institut fiir Praktische Mathematik der Technischen Hochschule Darmstadt (Prof. Dr. 4. Walther) 


(Eingegangen am 25. Juli 1959) 
Anschrift des Verfassers: L. Kraus, Wiesbaden, Idsteinerstrafe 11. 
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Uber den Einflu8 der Unterteilung von Schmierflachen 
auf die Tragfahigkeit von Schmierfilmen* 


Von H. Peeken 


1. Einleitung. Durch viele Untersuchungen theoretischer und auch experimenteller Art ist das 
hydrodynamische Verhalten des ebenen Schmierkeiles hinreichend geklart. Im Falle des unend- 
lich breiten ebenen Schmierkeiles geniigt der im Schmierfilm erzeugte Druck P bekannterweise 
der Beziehung von Reynolds (vergl. Abb. 1) ; 

dP H — H* 

qx = 67U <5 
Das heiBt, der Druckgradient enthalt die dritte Potenz der Spalthéhe H im Nenner. Bei groBen 
Spalthéhen wird also der spezifische Druck gering bleiben, wahrend er bei kleinen Spalthéhen 
groBe Werte annimmt. Aus dieser Aussage leitet sich fiir die Praxis sofort die F orderung nach 
engen Spalten ab, wenn die Erzielung einer méglichst groBen Tragkraft gewiinscht wird. Diese 
Forderung, die Schmierspalte klein zu halten, st68t dann auf Schwierigkeiten, wenn man sehr 
lange Lagerungen, wie sie z. B. in Werkzeugmaschinen auftreten, aus einem einzigen Schmierkeil 
herstellt. Der Grund fiir das Auftreten groBer Spalte bei groBen Lagerlangen muB darin gesucht 
werden, da der herstellbaren Schmierkeilneigung aus fertigungstechnischen Griinden eine Grenze 
nach unten gesetzt ist. Abb. 1 zeigt deutlich, daB8 der Druck am Keileinlauf zunachst sehr gering 
bleibt, um erst gegen Ende des Keiles zu technisch brauchbaren Héhen anzusteigen. Praktisch 
wird also bei einer Lagerausfiihrung nach Abb.'1 die gesamte Tragkraft gegen Ende des Keiles 
erzeugt, wahrend die Flache am Keileinlauf fiir die Druckerzeugung ausfallt. 


lk 
= Z ie 
= x YYLLL 3 
> 
=. ALLL) Ws YZ 
Ni J J N 
Xx U— x= 
vo 9 oy - 
st rq | 
p 
Abb. 1. Ebener Schmierkeil mit Druckverlauf Abb. 2. Unterteilter ebener Schmierkeil mit Druckverlauf 
(schematisch). (schematisch). 


Um hier Abhilfe zu schaffen, werden Lagerungen und Fiihrungen von gréBerer Lange aus meh- 
reren ,,Teilkeilen“’ aufgebaut. Die Anordnung einer ,,mehrflachigen“ Lagerung dieser Art ist 
gleichzeitig mit dem dazugehérigen Druckverlauf in Abb. 2 dargestellt. Schon seit langem wird 
in der Praxis das Prinzip der Mehrflachigkeit angewandt. Jedoch wird die Aufteilung der zur Ver- 
figung stehenden Gleitflachen meist vom Konstrukteur nach rein konstruktiven oder gefiihls- 
maBigen Gesichtspunkten vorgenommen ohne Kenntnis, wie sich im einzelnen Fall die Unterteilung 
auf die Tragfahigkeit auswirkt. Fiir die Praxis ist demnach die Frage von Bedeutung, wie ein 
ebener Schmierkeil von vorgegebener Lange und Breite zu unterteilen ist, um ein Optimum an 
Tragkraft zu erreichen. 

Dabei soll erwahnt werden, daB oftmals eine Unterteilung der Lagerflache nicht zur Erhéhung 
der Tragkraft, sondern zur besseren Olversorgung vorgenommen werden mu. Dieser Fall tritt bei 
sehr kurzhubigen Gleitschuhen auf. Die Lange des Hubes ist bei Gleitschuhen ein MaB fir die Teil- 
flachenlange. Wird die so bestimmte Teilflachenlange iiberschritten, so kann der Fall eintreten, 
da® Teile der Lagerflache nicht yom Schmiermittel benetzt werden. Dieser Gesichtspunkt soll je- 


* Auszug aus der von der Fakultat fiir Maschinenwesen der Technischen Hochschule Braunschweig 1959 
genehmigten Dissertation. 1. Berichter: Professor Dr. H. Schlichting, 2. Berichter: Professor Dr. G. Vogelpohl. 
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doch hier unberiicksichtigt bleiben. Den Untersuchungen liegt die Annahme zu Grunde, dal die 
Schmierspalte zu jeder Zeit mit Schmiermittel angefiillt sind. ; 

Das Schrifttum liefert zur Frage des Einflusses der Teilflachenzahl n auf die Tragfahigkeit nur 
sehr sparliche und sich zum Teil widersprechende Angaben. Wahrend von Giimbel! die Ansicht 
vertreten wird, daB in bestimmten Fallen eine Unterteilung des Schmierkeiles unendlicher Breite 
durchaus sinnvoll sein kann, halt Duffing? in jedem Falle die Unterteilung als Mittel zur Tragfahig- 
keitserhéhung fiir nicht geeignet. Bei Shaw und Macks? findet man die Aussage, dafs eine Unter- 
teilung der Laufflache immer zu einer Abnahme der Tragfahigkeit fiihrt. Fir den Fall des endlich 
breiten Schmierkeiles besteht unter den Autoren Talke*, Schiebel®, Drescher®, Fréssel’ und Niemann® 
darin Ubereinstimmung, da8 der maximale Flachendruck dann erreicht wird, wenn Lange und 
Breite des Schmierkeiles etwa gleich groB sind. 

Die Entwicklung der mehrflachigen Lagerung in der Praxis zeigt zwei Bauformen. Bei der 
weitverbreiteten ,,starren Bauform“ sind die einzelnen Keile durch eine spanabhebende Bearbei- 
tung der starren Lagerflache hergestellt. Die zweite Bauform ist nach Art des Michell-Lagers aus 
einzelnen Gleitsteinen aufgebaut, die auf einer balligen Flache gelagert, sich der Belastung ent- 
sprechend einstellen kénnen. Die vorliegende Untersuchung beschrankt sich auf die starre Lager- 
form. 

Nicht nur im Falle des ebenen Schmierkeiles, sondern auch beim Gleitlager findet man Baufor- 
men, die eine Vielflachigkeit aufweisen. Durch die Forderung des Werkzeugmaschinenbaues nach 
méglichst genauer Zentrierung der Wellen wurde besonders in neuerer Zeit die Entwicklung dieser 
Lager vorangetrieben. Daher sind gerade auf diesem Sektor des Maschinenhaues die verschieden- 
sten Bauformen mehrflachiger Gleitlager anzutreffen. Der Gesichtspunkt, wie im einzelnen Fall 
die Ausbildung der am Umfang gleichmaBig verteilten Schmierkeile erreicht wird, fihrt zu einer 
Einteilung der mehrflachigen Gleitlager in drei Gruppen: 

1. Die Schmierkeile werden durch elastische Verformung einer entsprechend gestalteten diinn- 
wandigen Lagerbuchse gebildet. Der typische Vertreter dieser Lagerbauart ist das Mackensen-Lager. 

2. Die einzelnen Schmierkeile werden durch eine entsprechende Verlagerung von Gleitsteinen 
erzeugt. Man hat hier das Prinzip des Michell-Lagers als Axiallager auf das radiale Gleitlager an- 
gewandt. Die Gleitflachen sind in diesem Falle nicht mehr Zonen einer geschlossenen Lagerbuchse 
wie noch beim Mackensen-Lager, 
sondern einzelne Segmente, die 
sich, auf balligen Leisten gelagert, 
je nach Belastung einstellen kén- 
nen. Diese Bauart trifft man in 
den verschiedensten Ausfiihrungs- 
formen meist im auslandischen 
Schrifttum an?. 

3. Bei der letzten Gruppe der 
mehrflachigen Gleitlager werden 
die einzelnen Gleitflachen durch 
entsprechende Formgebung einer 
starren Lagerbuchse hergestellt. 
In Abb. 3 ist ein Beispiel dieser 


ES Lagerbauart gezeigt. Es handelt 

GleitFlichenprofil sich dabei um ein fiinfflachiges 

Abb. 3. Fiinfflachiges konisches Gleitlager mit starrer Lagerbuchse konisches Gleitlager, das als Spin- 
(Fritz Studer A.G., Bern). dellager einer Werkzeugmaschine 


Verwendung findet. 


1 L, Giimbel u. E. Everling, Reibung und Schmierung im Maschinenbau, Berlin 1925. 
* G. Duffing, Schmiermittelreibung bei Gleitflachen von endlicher Breite, in: Auerbach-Hort ,.Handbuch 
der physikalischen und technischen Mechanik“ Bd. 5, S. 839, Leipzig 1929/30. é 
2 bes pone u. E. F. eae and Lubrication of Bearings, New York 1949. 
. Talke, Beitrage zur rodynamischen Schmiertheorie d i i 
POPE eye ek y y chmiertheorie des ebenen Gleitschuhes auf ebener Gleitbahn. 
5 A. Schiebel, Die Gleitlager, Berlin 1933. 
6 H. Drescher, Konstruktion 8 (1956), S. 94. 
‘ Ao yet ee Math. Mech. 21 (1941) S. 321. 
. Niemann, Maschinenelemente 1. Bd., Berlin/Géttingen/Heidelberg 1950 
° A. G. M. Michell, Lubrication, Principles and Practice, London 1951 se 
Analysis and Lubrication of Bearings, New York 1949. cheese ee Silas pees 
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Wie im Falle des ebenen Gleitschuhes erhebt sich auch hier die Frage nach dem Einflu8 der 
Unterteilung auf die Tragfahigkeit der mehrflachigen Gleitlager. Einige Hinweise, die diese Frage 
betreffen, finden sich im neueren auslandischen Schrifttum. So wird von Shaw und Macks besonders 
bei Lagern geringer Breite eine Aufteilung der Lagerflache in quadratische Teilflachen befiirwortet, 
um, wie die Verfasser angeben, auf diese Weise durch Herabsetzung des Seitenflusses eine Er- 
héhung der Tragfahigkeit zu erzielen. Eine weitere Untersuchung dieser Art wurde von Boyd und 
Raimondi} allerdings mit Lagern der Bauart 2 vorgenommen. Die Autoren kommen nach Diskus- 
sion der Zusammenhange zu dem SchluB, daB die Anwendung der untersuchten Segmentlager in 
gewissen Fallen vorteilhaft sein kann, daB aber eine bedeutende Uberlegenheit gegeniiber dem nor- 
malen Gleitlager, wie sie manchmal propagiert wird, keinesfalls besteht. 

Bei den mehrflachigen Gleitlagern ist demnach in erster Linie die Frage von Interesse, ob und 
in welchem Mae bei Gleitlagern geringer Breite durch die MaBnahme der Unterteilung eine Trag- 
kraftssteigerung zu erzielen ist. 


2. Der unterteilte unendlich breite ebene Schmierkeil. Abb. 4 zeigt einen unterteilten ebenen 
Schmierkeil fiir n = 4 Teilflachen und fiir eine Bewegungsrichtung. Es geniigt bei dieser Unter- 
suchung vollkommen, Schmierkeile mit nur einer Bewegungsrichtung zu behandeln, da sich La- 
gerungen fiir hin- und hergehende Bewegung sofort auf den dargestellten Fall zuriickfiihren lassen. 
Fiir die Behandlung des Problems miissen noch folgende Voraussetzungen gemacht werden: Die 
Schmiernuten stehen fortwahrend mit dem Schmiermittelzulauf in Verbindung, so daf jede Teil- 
flache geniigend mit Schmierstoff versorgt wird und der Druck P im Ol] mit Sicherheit an den 
Kanten der Schmiernuten auf den AuBendruck abfallt. Unterdriicke kénnen in der Schmiernut 
nicht auftreten. 

LaB®t man die Schmiernuten 
zwischen den einzelnen Teil- 
flachen unberiicksichtigt, so 
ergibt sich fiir die Gesamtlange 
des Schmierkeiles 


L=n-l. (1) 


Der Schmierspaltverlauf be- 
rechnet sich aus 


Hi (b = xX) <2) 


7 


S ES 
AuBerdem gilt nach Abb. 4 die > 
Beziehung 
Hy Abb. 4. Unterteilter ebener Schmierkeil fir n = 4. 
ee } 3 
pe eer (3) 


Der Druck im Schmiermittel des unendlich breiten Schmierkeiles laBt sich mit Hilfe der bekannten 


Gleichung von Reynolds 
dP H — H* (4) 


berechnen. Darin bedeutet 7 die Viskositat des Schmiermittels, U die Gleitgeschwindigkeit und H* 
eine Konstante, die die Spalthéhe an der Stelle darstellt, fiir die dP/dX = 0 wird (vergl. Abb. 1). 
Durch Einfiihrung des Spaltverlaufes nach (2) und durch Integration ergibt sich bei Beriicksichti- 
gung der Randbedingungen P = P, = 0 fiir X = 0 und X = /aus (4) die Beziehung fiir den Druck 
P(X) 

67 UX (== X) (6) 


HOM Eessp Or he 


Eine nochmalige Integration zwischen den Grenzen 0 und / liefert den Ausdruck fiir die Trag- 
fahigkeit T, einer Teilflache der Lange | 


ie 1 Bit eo} 6 
oO 2 (it pear eT: (6) 


1 J. Boyd u. A. A. Raimondi, Mechanical Engineering 75 (1953) S. 380. . 
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Fiihrt man in (6) die geringste Schmierschichtstarke H), die bei vélliger Reinheit des Schmier- 
mittels nur durch die Héhe der Rauhigkeiten gegeben ist, als konstant nach (3) ein, so ergibt sich 
mit 1 = L/n die gesamte Tragkraft T des Schmierkeiles zu 


Le 
Le One Le n Hy 7 
A ae In (1 + a) res ae (7) 
n Hy 
pill eke s : : : 
Es ist zweckmaBig, statt der Tragkraft T die Druckkennzahl Pare einzufiihren. Sie ergibt sich 
mit dem Flachendruck p = T/L aus (7) zu J 
Le 
Piel aos on ( Le n Hy 8 
or = =n ie a4 Saas é (8) 
“H2 n Hy 


Wie (8) zeigt, ist die Druckkennzahl als Ma8 fiir den mittleren Flachendruck p von der Teilflachen- 
zahl n und von der geometrischen GréBe A = L ¢/H, abhangig: 


A 


aa 6n in(1 fl n 
ey + —)—2 . (9) 
fas fe n 24 = 


Es liegt daher nahe, die Auswertung von (9) fiir verschiedene Teilflachenzahlen n in Abhangigkeit 
von A vorzunehmen. Die Ergebnisse dieser Auswertung fiir n = 1; 2; 3; 4; 5 und 6 und A = 1 bis 
100 sind in den Abb. 5 und 6 aufgetragen. Man erkennt aus dieser Auftragung, daB im Bereich 


Druckkennzoh! 
pig 
7UL 
016 : —— > al 
On 7} = 1 
n=1 
One 4 24 Al 
070 
008 a 
04 — 9 - 
oo} i 
a 
6 
ae 
| zp 
App = 1194 
; pie i per, | a= 
7 1194 45 (A 3 4 0. G=7-8 90 15 20 50 40° 50 60 70 8095100 


pe]. De 
Abb. 5. Druckkennzahl Ape? in Abhangigkeit-von 2 = a7 mit n als Parameter. 
° 


A = 1 bis A = 10,4 eine Unterteilung zu einer geringeren Druckkennzahl ae “| fiir die unend- 
4] 


liche Breite fiihrt. Erst fiir Werte von A, die gréBer sind als 10,4, bietet eine Unterteilung hinsicht- 
lich der Tragkraft Vorteile. Um zu erkennen, welche Tragkraftsteigerungen im Bereich 4 > 10,4 
durch Flachenunterteilungen méglich sind, sind in Abb. 6 noch einmal der Verlauf der Druck- 


spear plates 
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kennzahl fir n = 1 sowie diejenigen Druckkennzahlen fiir n = 2; 3; 4; 5 und 6, die die jeweils 
héchste Tragkraft liefern, aufgetragen. Aus dieser Darstellung ist zu ersehen, daB A den Base 
10,4 [A < 17,3 die Teilflachenzahl n = 2 , 
17,3 SA < 27,2 die Teilflachenzahl n = 3, 
27,2 SA < 35,0 die Teilflachenzahl n = 4, 
35,0 <A < 38,3 die Teilflachenzahl n = 5 
die gréBte Druckkennzahl und damit die groBte Tragkraft liefert. In Abb. 6 sind auBerdem die bei 
Wahl der optimalen Teilflachenzahl méglichen Tragkraftsteigerungen in % bezogen auf die Trag- 
kraft fiir n = 1 in Abhangigkeit vonA = L é/H, aufgetragen. Die Steigerungen betragen im Bereich 
n = 2 bis zu 15%, 
n = 3 bis zu 37%, 
n = 4 bis zu 53%, 
n = 5 bis zu 60%. 


Bei noch héheren Teilflachenzahlen nimmt die Tragkraftsteigerung weiterhin stetig zu, bis sie bei 
ganz groBen Teilflachenzahlen wieder abnimmt. 


Druckkennzahl 
Pit of 
Ut 
06 ~ 100 pee tat 
go} Magliahe Steigerung von ze int 
We go|- tur verschiedene 7 bezogen auf _| 
L. dle Druckkennzahl fiir n=1 
G72 
ai - 
G08 
Ant 
ae ; 
go4 
ae 
A gpt =1194 L 
h Aa=te 
| Pare 4 fo 
7 15 é 3 YOSSI EGO eR S910 15 ZOO 40 50 60 10 8 90100 
A= A=3 — A=tIl A=35 ; 


A=382 
Abb. 6. Verlauf der maximal méglichen Druckkennzahl. 


Die hier gefundenen hohen Prozentzahlen fiir die durch Flachenunterteilung mégliche Steige- 
rung der Tragkrafte diirfen jedoch nicht dariitber hinwegtauschen, daB es sich absolut gesehen in 
diesem Bereich um sehr kleine Tragkrafte handelt. Die bei weitem héchsten Tragkrafte ergeben 
sich ohne Unterteilung im Bereich A = 1 bis A = 10,4, wobei das Maximum bei A = 1,194 zu 


pis i. 
Pat| = 016 


gefunden wurde. Die dabei einzuhaltene optimale Schmierkeilneigung ergibt sich aus A = 1,194 zu 


Ay 
opt ES 1,194 Te ° 
Da der Bereich von A, in dem n = 1 die héchste Tragkraft liefert, eine Zehnerpotenz umfaBt, wird 
es sicherlich mit Hilfe der modernen Bearbeitungsverfahren immer miéglich sein, eine Keilstei- 
gung herzustellen, die die Bedingung ¢ = (1 ~ 10) H,/L erfillt. Unter dieser Voraussetzung kann 
festgestellt werden, daB eine Unterteilung einer unendlich breiten Gleitflache keine Vorteile bringt, 
; 14* 


Ingenieur-Archiv 
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sondern nur zu einer Verminderung der Tragfahigkeit fiihrt. Die eingehende Darlegung der Vor- 
ginge macht es jedoch méglich, auch die Grenzen anzugeben, in denen in bestimmten Fallen eine 
Unterteilung des Schmierkeiles zu einer Tragkraftsteigerung fiihren kann. 


3. Der endlich breite ebene Schmierkeil. Der Druck P im Schmiermittel eines endlich breiten 
ebenen Schmierkeiles geniigt der bekannten Differentialgleichung von Reynolds 


0 fH® OP\>( <0 (H* OP y @H _ 9 10 
ax ml al Be ox z (10) 


Die Ubereinstimmung dieser Grundgleichung der hydrodynamischen Schmiertheorie mit der Praxis 
ist durch zahlreiche Versuche bestatigt worden, so daB gleichfalls eine Ubereinstimmung der 
Ergebnisse dieser Untersuchung mit der Praxis erwartet werden darf. Fir die nachfolgenden 
Ableitungen erweist es sich als zweckmaBig, folgende dimensionslose GroSen einzufiihren: 
L ree 


x E 3 H => diy S- Lig __ 
= S — =h, =) ese e r if Pp. 


L am ~~ 


Der Schmierspaltverlauf ergibt sich zu 
H =¢(E—X) 
(12) 


h=e(e—x). 


oder dimensionslos 


asf 


AuBerdem folgt aus Abb. 7 


H, =eL, =e(E—L) 
(13) 


oder 
Abb. 7. Bezeichnungen am Schmierkeil endlicher Breite. hy =ép =£ (e — 1) 3 


Das halbe Breitenverhaltnis eines Schmierkeiles von endlicher Breite ergibt sich zu 


B 

[= oT (14) 
Fiihrt man die dimensionslosen GréBen nach (11) in die Reynoldssche Differentialgleichung (10) 
ein, so erhalt man eine dimensionslose Form dieser Gleichung: 


Of OP ad ép oh 
aa \™ ze(" rice 0 (15) 


mit den Randbedingungen: 
¢ == 6° und alles 


p=n=0{s=0 und alle€. 


all 


Zur Lésung dieser Gleichung soll hier ein Verfahren angewandt werden, das von Weber? fiir die Be- 
rechnung von Gleitlagern angegeben wurde. Es stellt ein Naherungsverfahren dar, hat aber hier 
den Vorteil, da sich die Lésungen geschlossen ausdriicken lassen. Wird die Reynoldssche Dif- 
ferentialgleichung (15) als Eulergleichung eines Variationsproblems aufgefaBt, so lautet das zu- 
gehérige Variationsproblem 


Behan ties, 5 . 
3 (oP 3 (oP op 
J [iG 38s 
= 2.0 
Die Lésung des Variationsproblems gelingt mit dem Naherungsansatz 


P =f (x) a) (17) 


g(t B) =0, 
FO) =S0% 
SO ieo 0 


1 Z. B. W. Niicker, VDI-F orschungsheft 352, Berlin 1932. — H. Drescher, Konstruktion 8 (1956), S. 94. 
® C, Weber, Z. angew. Math. Mech. 30 (1950), S. 112. 


dx dé = Min. (16) 


samt den Randbedingungen 


XXIX. Band 1960 


H. Peeken: Unterteilung von Schmierflachen auf die Tragfahigkeit von Schmierfilmen 205 


Schreibt man zur Abkiirzung df (x)/dx = f’ und dg(¢)/dt = g’, so folgt aus (16) bei Einfihrung des 
Naherungsansatzes (17): 


pad 
i Jie ie + fg —2hef) dx d& = Min. (18) 


Setzt man eine beliebige Funktion f(x), die die Randbedingungen f(°) = 0 erfiillt, in das Variations- 
problem ein, so kann man nach x integrieren und erfahrt, welche Bedingungen die Funktion g(¢) 
erfiillen mu, damit das Integral zu einem Minimum wird. Die bei dieser Operation auftretenden 
bestimmten Integrale werden wie folgt bezeichnet: 


1 l 1 
mafhfide, afi frde, a =fhftdr. (19) 
9 0 


Setzt man die bestimmten Integrale a), a, und ay, in (18) ein, so erhalt man 


+B 
if [a, g° + a, g’*— 2 ag] dC = Min. 


Aus diesem Variationsproblem folgt-die Eulergleichung 


deren Lésung sich unter Beriicksichtigung der Randbedingungen zu 


g(t) = C, (1 ean ) (20) 


mit a2 =a,/a, und C, = a,/a, ergibt. Es muB8 also, um eine optimale Lésung fiir den Ansatz 
p =f (x) g(¢) zu erhalten, fiir jede beliebige Funktion f(x) die Funktion g(¢) (20) gewahlt werden. 
Die gleiche Untersuchung wird jetzt fiir f(x) durchgefiihrt. 

Kine beliebige Funktion fiir g(C), die die Randbedingungen g(+ f) = 0 erfillt, in das Varia- 
tionsproblem (18) eingesetzt, l48t die Integration nach ¢ zu. Die dabei auftretenden bestimmten 
Integrale werden wie folgt bezeichnet: 


ve aa +6 
by =a g(c) dc , b, aah g(C) dC, by me g’(C) de. (21) 


Es mu8 folgendes beachtet werden: Die Naherungsfunktion p(x,¢) ist das Produkt der beiden 
Funktionen f(x) und g(¢). Bei der Zerlegung der Funktion p in diese beiden Faktoren kann die 
eine Funktion mit einer Konstanten multipliziert und die andere Funktion durch die gleiche Kon- 
stante dividiert werden, ohne daf sich das Produkt p(x,¢) dabei andert. Wird nun g(¢) mit (b)/b,) 
multipliziert und 6,/b, g(¢) fur g(C) in das Gleichungssystem (21) eingesetzt, so ergeben die ersten 
beiden bestimmten Integrale dieses Gleichungssystems denselben Wert. Man kann daher unbe- 
schadet der Allgemeinheit des Ansatzes annehmen, dab g(C) so gewahlt ist, daB 6) = b, wird. 
Fiihrt man nun 6, = b, und J, in (18) ein, so erhalt man 
1 


[ (ret RBs) de = Min 
1 
| 


Aus diesem Variationsproblem folgt die Eulergleichung fiir die Funktion f(x) zu 
Ee, fees dh 
EAA aaa Sete 
Differentiiert man aus und setzt die Beziehung fiir die Schmierspalthéhe ein [Gl. (12)], so wird 
3 rd 1 


(aA tA 


ae xe be ee — ay 


Setzt man b,/b, = 1/x und c? = e?/x und transformiert die Gleichung mit Hilfe von (e — x) Vx hs 
auf die neue Veranderliche €, so ergibt sich 


of Sin, eae 1 
fA ef fae (2) 


206 HH. Peeken: Unterteilung von Schmierflachen auf die Tragfahigkeit von Schmierfilmen —_ Ingenieur-Archiv 


Die Randbedingungen dieser Gleichung lauten mit e— 1 = 9 und e = +1 [vergl. (13)] 


fee ie (23) 
50 cs | ara 


Mit Hilfe einer zweiten Transformation erhalt man die Besselsche Differentialgleichung 


tas 1 be 1 — 1 
bs oe. —(l spe fae 
mit (24) 
y St 
f= af oder f= ae? 
deren homogene Lisung die Zylinderfunktion vom imaginaéren Argument und vom Index 1 ist. 
y =4 (16) 
mit 
Z,(i €) = hy Tv 5) + hy Ny (6) - 
Die inhomogene Lésung kann mit Hilfe des Ansatzes 
y= Dm, EC? 
0 
zu : 
1 1 
PEER ao ee Eee eS) 


ermittelt werden. Damit ist bekannt, wie die Funktion f(&) zu wahlen ist, um das Variationsproblem 


(18) zu erfiillen 


f®) =a; ZG) +I - (26) 


Mit Hilfe dieser Gleichung la8t sich der Druckverlauf des Schmierkeiles endlicher Breite geschlossen _ 
ausdriicken: 


c 
P= ayey AWE) +90] (1 Got p) 


oder 


Cof at 

P= a z Ue LEE) + ba NAGE) +991 (1 Gotz) (27) 
Die hier vorliegende Lésung fiir den Druck p im Schmiermittel wurde mit Hilfe eines Produkt- 
ansatzes der Form p = f(x) g(¢) unter gleichzeitiger Erfiillung des Minimalproblems nach Gl. (16) 
bestimmt. Man kann aus diesem Grunde die hier vorliegende Lésung als eine optimale Naherungs- 
lésung bezeichnen. Sie stellt die optimal mégliche (wegen Erfiillung des Minimalproblems) inner- 
halb der Bewegungsfreiheit des Ansatzes dar. Die Frage, inwieweit ein solcher Ansatz in der Lage 
ist, die tatsachlichen mechanischen Vorgange zu beschreiben, ist nicht einfach zu beantworten. 
Produktansatze dieser Form haben sich jedoch bei vielen ahnlichen Problemen zur Gewinnung 
von Naherungslésungen gut bewahrt!. 

In der zur Verfiigung stehenden Funktionstafel von Jahnke-Emde? findet man die Funktionen 


[—i (i &)] und ee HOY (i 5) tabelliert. Es ist daher zweckmaBig, die Lésung (26) auf die 
tabellarisch vorhandenen Funktionen umzuschreiben. Man erhalt dann fiir die Funktionen f (é) 


eee: {Coli HGH] + Gy Ls HM | + 99} 


7 H? 
oder mit ¢2 = 


FTE = Ee {Gi ES] + G [5 BMG 8] + al. (28) 


* Z.B. C. Weber, Z. angew. Math. Mech. 30 (1950), S. 112. — C. Weber, Profilriicknahme bei schragver- 
zahnten Radern, Heft 11, Bericht 141 der Schriftenreihe Antriebstechnik, Braunschweig: Vieweg 1953. 
2 E. Jahnke u. F. Emde, Funktionstafel mit Formeln und Kurven, Leipzig 1948. 
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Fiir den Druckverlauf ergibt sich damit 


a ree ; : % rr: Coj a 
Pye (AEF REA + G|—F HPCE + 9, 1G’, 
Bei dieser Umschreibung besteht zwischen C,, Cs, ky und k, der Zusammenhang 
Cj=ih—k, und 6, = 2th 


Ersetzt man den dimensionslosen spezifischen Druck p durch den wirklichen spezifischen Druck P 
nach (11), so ergibt sich mit 2 = H2/L? g? 


P H? C, 9? Pet ce Cojae 
Die Konstanten C, und C, bestimmen sich aus der Randbedingung (23) 
ae Jey oa 
Ve” 
== () 
f Y 


Zu 


(30) 


C, 
: peat ee rae iy Ae 
Bea) ea ral 
. -@P + 1 M4 Q) [+ Pp + 1 
ae eee Ee 
Fur die vorliegenden Betrachtungen interessiert weniger der Druckverlauf P als die GréBe der 


Tragkraft fiir einen Schmierkeil von vorgegebener Lange und vorgegebener Breite. Die Trag- 
kraft ergibt sich zu 


ro| & 


Te 
J f[ PaXd 
0B 

2 


I 


qT, 


oder mit P = p67 U/L, dX = L dx und dy = Ld¢ 


sal 1 +8 
te | [owe =| [ forse) deat. 


0 —B 0 —8 


67 


Mit Beriicksichtigung der Transformationsgleichung € = e — «| Ve wird 
g+1 
Ve +6 


T. i 
t= yu [| [ fie)ele)dtae. 
6nUL of J 


Ve 
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+6 
Fiihrt man nach ( ) fal 6) dé = by ein, so wird 
(3 


GEL = x by Jw 


Die Integration von f(&) nach (26) iiber & liefert den Ausdruck fiir die Tragkraft T, 


aa 
bo 
sli = heme ae—ae + a0) en 
Vu 
Dabei ist 
Ale) =1n6-+ 5 P 4 a pace rae oe) 


Schreibt man (31) auf die in der Funktionentafel Jahnke-Emde tabellarisch vorhandenen Funk- 
fone (7 LG S|, [HG OL ES HMI 8) und EB nn HOI is) um, so ergibt sich 


srr ~ % {Ca [Cues)) a — (i Le8) 
oan Gs | | Gey HO 5) dé + Gee H(i 5)| < ae 


Fiuhrt man in diese Gleichung noch den mittleren Flachendruck 


ee i Petes if : B 
PBL” € erie sowie boa 

BH 
ynUB 


ein, so gewinnt man damit eine Beziehung fir die Druckkennzahl 
ebenen Schmierkeiles 


BU |-IBs [ol frene—ceneey 


| des endlich breiten 


— 6] fig HOGS) ds +(—F HOGS) + ae} (33) | 
Vx 


Zur Ermittlung der Druckkennzahlen nach (33) werden noch die bestimmten Integrale 6), b, und 
b, benétigt. Sie ergeben sich zu 


4 = fof er 2) a = 2G, 2p Tq af), 


ap 
ts — let aco B 
: 4 i( Soap) E~ OS q (206 —38g08 +h, tap)? (34) 
et ate : Cia? 
a) Cojrap Gin? «C) dé Cee Gin 2 «0 B — 6). 
Die Bedingung 6, = b, liefert eine Bestimmungsgleichung fiir C,: 
C= Sn (35) 


2aB—3 ga + ——— CoPae 
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Fir x = b,/b, erhalt man aus den Gleichungen (34) 
Pe. /. 2 Cinta B 
cae a (sueeOoise 1 22 (36) 
an eee 


Die Ausdriicke fiir C,, x sowie by sind Funktionen des Produktes « 6. Nimmt « f gréBere Werte als 
4.6 an, so kénnen fiir die genannten Ausdriicke vereinfachte Beziehungen angegeben werden. 


Der dabei auftretende Fehler ist kleiner als 0,01%. Es gilt fiir o B > 4,6: Tq ape Il, Sof a Bx 


Gin «a B und Stag ® 0 und damit 
«aB—l 


Se ae 


1 2 =e] )2 
x= af —3), b= = a. (37) 


Das in diesen Gleichungen vorkommende a ergibt sich aus der folgenden Beziehung: 


Darin sind a, und a, bestimmte Integrale nach (19): 


party 

a, = ff? dx, 
0 
1 

a, = fh? f? dx 
0 


Die Transformationsgleichung e— «|x = € liefert zusammen mit (12) den Ausdruck fiir den 
Schmierspalt h in Abhangigkeit von & 


he ae (38) 


Setzt man f’ aus (26) zur Berechnung des bestimmten Integrales a, und f nach (26) zur Berechnung 
des bestimmten Integrales a, in die Gleichungen (19) ein, so fiihrt die Integration auf Zylinder- 
funktionen vom Index 2 und 3. Da diese Funktionen tabellarisch nicht greifbar sind, werden fiir 
die Auswertung die bestimmten Integrale a, und a, numerisch berechnet. Da selbst f’ eine Funktion 
von Z,(i €) ist, wird f’ durch numerische Differentiation aus der Funktion f(&) gefunden. Zur 
numerischen Differentiation lat sich die Stirlingsche Formel verwenden. Die Integration zur 
Berechnung der bestimmten Integrale a, und a, kann mittels der Simpsonschen Regel durchgefiihrt 
werden. Mit Hilfe der nunmehr vorliegenden Beziehungen lat sich die Tragkraft bzw. die Druck- 
kennzahl fiir den ebenen endlich breiten Schmierkeil berechnen. 

PG 
hiv 
g und B bekannt bzw. frei wahlbar sind, mufs die Auswertung mit der Bestimmung der x-Werte 
beginnen. Man erhalt sie durch iterative Rechnung, die gut konvergiert, so da man bei geniigend 
genauer Anfangsschatzung des x-Wertes meist schon mit einem Schritt auskommt. Da f fir 
die Berechnung der Druckkennzahl angenommen wird, ist x nach (36) eine Funktion von a, und a? 
ist, wie oben erwahnt, der Quotient der bestimmten Integrale a, und ag. 

Zur Bestimmung des x-Wertes fiir einen Schmierkeil einer bestimmten Breite f und eines be- 
stimmten Wertes L,/L = p geht man zweckmaBigerweise so vor, daB man zunachst einen Wert fiir 
~ annimmt und mit diesem geschatzten Wert die Konstanten C, und C, nach (30) sowie die Funk- 
tion f(&) nach (28) berechnet. Die zur Bestimmung der Konstanten C, und C; sowie der Druck- 


funktion f(&) benétigten Funktionen [— i (7 §)], —+ H(i €)| kénnen der Funktionentafel von 


4. Auswertung. Die Druckkennzahl nach (33) hangt von den GréBen g, 6 und x ab. Da 


Jahnke-Emde entnommen werden. Die zu diesem Zweck ebenfalls erforderliche Funktion y, wurde 
nach (25) berechnet und in Tab. 1 niedergelegt. Die Ableitung f’() kann nach der Stirlingschen 
Formel oder auch durch graphisches Differentiieren bestimmt werden. Mit den bekannten Werten 
fiir f(é), f’(€) und h(E) lassen sich dann die Integranden der bestimmten Integrale a, und a, berech- 


nen. Die Integration in den Grenzen 
g+l1 


Vx 


ee ee eee 


yx 


=O > = 
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Tabelle 1. Funktionen y,(&) und A(&) nach Gl, (25) und (32). 


0,00 + 1,00000 — oo 1,90} +2,52453| +1,3211 | 3,90] + 14,0767 | + 5,6314 

0,05 | + 100083 | —2,9954 | 2,00] + 2,7321 | 41,4558 | 4,00 | + 15,4022 | + 6,0030 
0.10 | +1,00333 | —2,3009 | 2,10} + 2,9607 | +1,5945 | 4,10 | + 16,8576 | + 6,4033 

0.20. + 1,01337 | —1,6028 | 2,20 | +3,2117 | +1,7378 | 4,20 | + 18,4529 | + 6,8251 

0.30 | +1,03018 | —1,1890 | 2,30) +3,4876 | +1,8866 | 4,30 | + 20,2035 | + 7,2806 

0,40 105390 | —0,8895 | 2.40 | +3,7904 | +2,0415 | 4,40 | + 22,3152 | + 7,7658 

0:50 1 41,0847) 4 La =20651 If 2504 Ada +2.2027 | 4,50 | + 24,2305 | -+ 8,2886 

0.60 | +1,12290 | —0,4501 | 2,60| +4,4876 | +2,3716 | 4,60 | + 26,5291 | + 8,8551 

0,70 | +1,16867 | —0,2737 | 2,70 | + 4,8875 | + 2,5482 | 4,70 | + 29,0783 | -++ 9,4413 

0,80 | +1,22260 | —0,1176 | 2,80 | +5,3263 | +2,7339 | 4,80} + 31,8623 | + 10,0888 
0,90 | +1,28494 | +0,0332 | 2,90 | +5,8071 +2.9293 | 4,90| + 34,9174 | + 10,7713 
1,00 | +1,35618 | +0,1724 | 3,00 | + 6,3350 | +.3,1346 | 5,00 | + 38,2707 | + 11,5063 
1,10 | +1,43696 | +0,3055 | 3,10} +6,9139--| +3,3515 | 5,20 | + 45,9949 | + 13,1668 
1,20 | +1,52797 | +0,4341 | 3,20) 47,5484 | +3,5813 | 5,40 | + 55,3022 | + 15,0585 
1,30 | +1,62985 | +0,5603 | 3,30) +8,2436 | + 3,8236 | 5,60 | + 66,5290 | + 17,2613 
1,40 | +1,74360 | +0,6875 | 3,40 -+9,0064 | +4,0829 | 5,80 | + 80,0652 | + 19,8448 
1,50 | +1,87024 | +0,8102 | 3,50) +9,8429 | +4,3540 | 6,00 | + 96,3866 | + 22,9004 
1,60 | +2,01005 | +0,9354 | 3,60, + 10,7595 | + 4,6434 

1,70 | +2,16500 | +1,0614 | 3,70 | +11,7648 | + 4,9504 

1,80 | + 2,33602 | +1,1899 | 3,80 | + 12,8669 | + 5,2810 


wird mit Hilfe der Simpsonschen Regel durchgefihrt. Sie liefert die bestimmten Integrale a, und 
a,. Aus der Division von a, und a, folgt a2. Das mit dem so bestimmten a aus (36) ermittelte x 
mu8 mit dem geschatzten Wert fiir x tibereinstimmen, gegebenenfalls ist die Rechnung zu wieder- 
holen. Ist eine Ubereinstimmung des geschatzten ~-Wertes mit dem gerechneten erreicht, laBt sich 
jetzt die Druckkennzahl nach (33) berechnen. 


ce ee ; p HG 
Beispiel: Es soll die Druckkennzahl ore 


fir B =0,5 und y = 0,8 berechnet werden. 


x wird geschatzt zu x = 0,0961, x = 0,31. Fir die Integrationsgrenzen ergibt sich o|Vx = 2,5800 
und (pg + 1)//x = 5,806. Mit diesen Werten werden die Konstanten C, und C, nach (30) be- 
stimmt. Es ergibt sich: C, = — 1,570977 und C, = — 2,399772. Die Ergebnisse der Berech- 
nungen von ¢, h L/H), f (&) H,/L?, f'(€) H,/L? sowie der Integranden (H,/L) f?(€) h? und (H,/L) f’2(&) h® 


der bestimmten Integrale a, und a, sind in Tab. 2 niedergelegt. Hieraus lassen sich die bestimmten 


Tabelle 2. Ergebnisse der Berechnungen von &, WT f= m 5 Foes Tt , sowie der tiasstaiden en Te h? und — Hy 7 f? h3, 
hL hL\3 He : | 
« ; 47 (=) fe) Fe) Fo pyeyns | Sepayey ns 
——— ee EE EE ee a a eee et 
0,0 5,8060 2,250 11,390 0,00000 — 0,01460 0,000000 0,002427 
0,1 5,4834 R125 9,596 0,004.69 — 0,01492 0,000211 0,002136 
0,2 5,1608 2,000 8,000 0.00963 — 0,01583 0,000742 0,002004 
0,3 48382 1,875 6,591 0,01456 — 0,01562 0,001397 0,001608 
0,4 4,5156 1,750 5,309 0,01969 — 0,01550 0,002077 0,001287 
0,5 4,1930 1,625 4,291 0,02436 — 0,01518 0,002546 0,000988 
0,6 3,8704 1,500 3,375 0,02821 — 0,00872 0.002685 0,000256 
0,7 3,9478 1,375 2,600 0,02967 + 0,00101 0,002288 0,000002 
0,8 3,2252 1,250 1,953 0,02660 + 0,02045 0,001381 0,000816 
0,9 2,9026 (hols) 1,424 0,01686 + 0,04075 0,000405 0,002364 
1,0 2,5800 1,000 1,000 0,00000 + 0,07130 0,000000. 0,005083 


Integrale a, und a, und daraus a berechnen: 


Aus (36) folgt dann x = 0,0995. 


el 


a 


oo cee & 6 0,526. 


- 0,3226 - 0 013732 , 


- 0,3226 - 0,015213 , 
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Wiirde man die Rechnung noch einmal wiederholen, so wiirde das neu berechnete x sich kaum 
von dem jetzt gefundenen unterscheiden. Auch der Wert fiir « wird in einem weiteren Iterations- 
schritt nur unwesentlich verandert. Mit dem x-Wert von x = 0,0995 kann jetzt die Druckkenn- 
zahl (33) berechnet werden. Es ist 


x=0,0995,  Vu=0,3154, g~=08, f=0,5. 


Zur Auswertung der Gleichung (33) muf noch die Funktion A(é) (32) bekannt sein. A(é) wurde 
fiir die Werte € = 0 — 6 berechnet und in Tab. 1 niedergelegt. 


Mit @/j/x = 2.5364 und (p + 1)//% = 5,707 ergibt sich 


% . @ 
il, (i ell = 26022, ees Ho) (i | = 0,07076 , 
= eS (i ASE | = 46,7512 , a Hw (i Pt ‘| — 0,001853 , 
2 i 
Pp 
—=|= 4,25368 , A({——|= 2,26418 , 
+(e) (i) 
Vp q as * = 73,47441 , A (? _ = 18,63152 
\x x 
Die Konstanten C, und C, werden neu zu C, = — 1,571512 und C, = — 2,321784 berechnet. Fiir 
die Integrale [ I,(i €) dé und fi + H(i ) dé ergibt sich aus der Funktionentafel von Jahnke-Emde: 
5,107 
f 1G &) dé = 54,5234, 
2,5364 
5,707 
[ i> HOGS) dé = 0,050669 . 
2,6364 2 


Mit « B-= 0,526 folgt aus (35) 

C, = 10,04714 
und nach (34) 

by = 0,834526 . 


Mit diesen Werten wurde die Druckkennzahl (33) zu 


pe. 
E = | — 0,0688 


gefunden. 


Druckkennzah/ 
py 


0 G2 Gt 06 08 40 42 


Abb. 8. Druckkennzahlen fiir den endlich breiten Gleitschuh in Abhangigkeit von 
mit dem Breitenverhiltnis # als Parameter. 


i i I ei fiir die Tragkraft des Schmier- 
Die auf diese Weise vorgenommene Auswertung der Gleichung 
keiles endlicher Breite erstreckt sich auf die y-Werte gp = 0,4; 0,6; 0,8; 1,0 und auf die Breiten- 
verhaltnisse B = 0,25; 0,4; 0,5; 0,8; 1,0. 
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Das Ergebnis dieser umfangreichen Berechnungen ist in den Abb. 8 und 9 dargestellt. Es ergibt 
sich aus diesen Auftragungen, daB die groBte Druckkennzahl und damit der gréBte Flachendruck Pp 
fiir die untersuchten Breitenverhaltnisse im Bereich y = 0,8 bis 0,6 auftritt. Das optimale Breiten- 
verhiltnis ergibt sich aus dieser Untersuchung zu B = 0,5. Der maximal erreichbare Flachendruck 
p und damit die gréBte Tragkraft des endlich breiten Schmierkeiles wird sich demnach bei einem 
Breitenverhaltnis 8 = 0,5, d.h. B/L = 1 und fiir p = 0,8 = 0,6 einstellen. Man wird daher die 
Unterteilung eines endlich breiten Schmierkeiles stets so vornehmen, daB quadratische Teilflachen 
entstehen. Da das Tragkraftmaximum in Richtung der beiden Veranderlichen PB und @ recht flach 


Druckkennzahl 


G02 + 4 


Breitenverhilinis 


0 Of 02 93 OF 5 06 G7 G8 99 40 


Abb. 9. Druckkennzahlen fiir den endlich breiten Gleitschuh in Abhangigkeit von 8 mit ¢ als Parameter. 


verlauft, ergeben etwas abweichende Werte fiir 6 und g noch recht hohe Tragkrafte. Vergleicht 
man die nach diesem Verfahren berechneten Druckkennzahlen mit den nach der exakten Lésung 
von Fréssel! berechneten Werten, so ergeben sich hier etwas héhere Flachendriicke. Der Unter- 
schied betragt etwa 24%. Das Ergebnis dieser Untersuchung, das maximale Flachendriicke beim 
endlich breiten Schmierkeil dann erreicht werden, wenn ein Seitenverhaltnis der Teilflache B/L = 1 
und ein KeilsteigungsmaB ~ = 0,8 — 0,6 eingehalten wird, stimmt mit dem im einleitenden Ab- 
schnitt zitierten Schrifttum iiberein. 


5. Das unterteilte Gleitlager. Bei der Betrachtung der unterteilten Gleitlager ist im Rahmen 
dieser Arbeit in erster Linie die Frage von Interesse, ob bei schmalen Lagern durch Unterteilung 
der Gleitflache eine Tragkraftsteigerung zu erreichen ist. Die Beantwortung dieser Frage gelingt 
durch einen Vergleich der Tragkrafte bzw. Belastungskennzahlen zwischen dem unterteilten und 
dem zylindrischen Gleitlager. Der Vergleich soll hier an Halbschalen, d. h. an 180°-Lagern durch- 
gefiihrt werden. Dabei wird fiir das mehrflachige Gleitlager die in der Einleitung angefiihrte Bau- 
form 3 gewahlt. Bei der Bauform 3 sind die einzelnen Gleitflachen durch entsprechende Form- 
gebung einer starren Lagerbuchse entstanden. 

Tragt man den Verlauf der Schmierspaltdicke bei einem mehrflachigen Gleitlager iiber der Teil- 
flachenlange auf, so stellt man fest, daB der Unterschied gegeniiber einem ebenen Schmierkeil mit 
linearem Spaltverlauf nicht groB ist. Man kann also die Teilflachen des Lagers in erster Naherung 
als ebene Gleitschuhe von endlicher Breite auffassen. Mit dieser vereinfachenden Annahme kénnen 
die Tragkrafte unterteilter Gleitlager relativ einfach ermittelt werden. 

In Abb. 10 ist ein unterteiltes Halblager fiir n = 6 Teilflachen dargestellt. Die maximal még- 
liche Tragkraft eines solchen Lagers wird dann erreicht, wenn man fir die als ebene Gleitschuhe 
aufgefaBten Teilflachen optimale Verhaltnisse wahlt. Die im vorhergehenden Abschnitt durch- 
gefiihrte Untersuchung hat ergeben, daB die maximale Tragkraft eines ebenen Gleitschuhes von 
endlicher Breite dann erreicht wird, wenn er ein Breitenverhaltnis B/L = 1 und ein Keilsteigungs- 
maB g = L,/L = 0,8 aufweist. 

Um das Maximum an Tragkraft beim mehrflachigen Gleitlager zu erhalten, denkt man sich ein 
solches Lager aus Gleitschuhen vom Breitenverhaltnis B/L = 1 und der Keilneigung g = 0,8 auf- 


1 W. Fréssel, Z. angew. Math. Mech. 21 (1941) S, 321. 
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gebaut. In diesem Fall ergibt sich fiir die Druckkennzahl einer Teilflache der Maximalwert P Hy 
= 0,069. Die Tragkraft einer Teilflache wird dann ae 


ULB? ; s 
T, = 0,009" mit p=. da B=L. 


Die Gesamttragkraft des mehrflichigen Lagers ergibt sich durch vektorielle Addition der Teil- 
krafte T,. Vernachlassigt man die von der Welle auf das Schmiermittel ibertragenen Schubkrifte, 
so missen die Teilkrafte T, durch den Wellenmittelpunkt gehen. Da simtliche Teilflachen den 
gleichen Schmierspaltverlauf haben, betragt der Winkel zwischen den Teilkraften a/n. Die Ge- 
samttragkraft T ergibt sich dann immer als Eck- 
maf} aus dem mit T, als Seite gebildeten (2 n)-Eck. 
Es gilt daher 


Setzt man fiir den mittleren Flachendruck p,, = 
T/B D und fir das Breitenverhaltnis des Lagers 
o = B/D, so ergibt sich der auf den Durchmesser D 
bezogene mittlere Flachendruck eines mehrflachi- 
gen Halblagers zu 


2 
0,069 « . (39) 


=F0—s) 


Pm Hj 
nUD 


Die fir eine Teilflachenlange einschlieBlich 
Schmiernut zur Verfiigung stehende Strecke be- 
tragt 2 D/2n. Nimmt man an, daf von dieser 
Strecke nur der i-fache Betrag als Teilflachenlange 
genutzt werden kann, dann ergibt sich die Lange 


einer Teilflache 


i 


lO] 


Mit B = L 1aBt sich hieraus das BreitenverhAltnis 


des Lagers berechnen 
Abb. 10. Bestimmung der Tragkraft fiir ein unterteiltes Halblager 
Beit (40) der Teilflachenzahl n = 6. 
0 => ;F >>: 
D 2n 


Setzt man zur Auswertung der Gleichung (39) i = 0,8, d. h. 20% der zur Verfiigung stehenden Ge- 
samtlange werden fiir die Schmiernutbreite verwandt, so lassen sich fiir verschiedene Teilflachen- 
zahlen die maximal méglichen Druckkennzahlen des unterteilten Halblagers nach (39) berechnen. 


Es ergibt sich 


n> Baseline A 5 6 


ey 8 
| 


o 0,420 | 0,314 0,252 0,210 0,180 0,157 


[para | 0,02435 | 0,0178 | 0,01418 | 0,01177 “ooo 0,00875 
0 : | 


Die Tragkraftkennzahlen des unterteilten Halblagers nehmen mit zunehmender Teilflachen- 
zahl ab. Fiir den Vergleich des unterteilten Halblagers mit der zylindrischen Halbschale kann die 
bekannte, Veréffentlichung von Sassenfeld und Walther! herangezogen werden. Hier liegen umfang- 
reiche Berechnungen der Sommerfeld-Zahlen So in Abhangigkeit vom Breitenverhaltnis o = B/D 


1 H, Sassenfeld u. A. Walther, Gleitlagerberechnungen, VDI-F orschungsheft 441, Diisseldorf 1954. 
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und von der relativen Exzentrizitat y vor. Die die Belastung kennzeichnende GréBe ist beim zylin- 
drischen Halblager die Sommerfeld-Zahl So. Es gilt 


noe Pa gs emt 8 Se oias 
SBD y a Y im R Y Ae Re Re 


Hierin bedeutet P die Tragflahigkeit des Gleitlagers in kg, y das relative Lagerspiel, R, den Radius 
der Lagerschale, R den Wellenradius, @ die Winkelgeschwindigkeit und ¢ die Exzentrizitat der 
Welle. Zum Vergleich beider Lagerarten ist es erforderlich, die Sommerfeldzahl in die hier ver- 


2 . . . 
wendete Tragkraftkennzahl Feed umzurechnen. Dabei ist zu beachten, daB ein Vergleich beider 
7] 


Lagerarten nur dann sinnvoll ist, wenn er auf der Basis gleicher minimaler Schmierschichtdicke Hp 
erfolgt. Mit 


oe pD a Hy ae P 
US Be ee ag es 
ergibt sich die Umrechnungsgleichung zu 
ie ae (41) 
nUD 2 


Mit Hilfe der so gefundenen Zusammenhange kann der Vergleich zwischen beiden Lagerarten vor- 
genommen werden. In Abb. 11 werden die Druckkennzahlen der unterteilten Halbschale denen 
der zylindrischen Halbschale gegeniibergestellt. Die letzteren ergeben sich durch Umrechnung der 
Sommerfeld-Zahlen aus der Veréffentlichung von Sassenfeld und Walther nach (41). Bei der Auf- 
tragung der Druckkennzahlen fiir das zylindrische Halblager in Abb. 11 erscheint die relative 
Exzentrizitat y als Parameter. y kennzeichnet in diesem Fall diejenige relative Exzentrizitat, bei 
der das zylindrische Lager die minimale Vergleichsschmierfilmdicke H, erreicht. Die Tragzahlkurven 
Kurven der grofen y-Werte stellen demnach die Lager mit groSem Lagerspiel dar, wahrend die 
der kleinen y-Werte Lager mit geringerem Lagerspiel bedeuten. 
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Abb. 11. Mittlere Flachendriicke fiir die unterteilte Halbschale und die zylindrische Halbschale mit y als Parameter. 


Die Auftragungin Abb. 11 zeigt, daB bei den y-Werten 0,95; 0,9; 0,3; 0,2, d. h. fir Lager der 
groSen und der kleinen Lagerspiele im ganzen untersuchten Bereich (0,42 <o < 0,15), das unter- 
teilte Lager hinsichtlich der Tragkraft im Vorteil ist. Die mégliche Tragkraftsteigerung, die das 
unterteilte Lager gegeniiber dem jeweiligen zylindrischen Lager zulaBt, ist in Abb. 12 dargestellt. 
Sie betragt in diesem Fall bei geringer Lagerbreite o = 0,15 und 7% = 9,2 ber 200%, wahrend 
sie bei y = 0,95 und o = 0,15 ca. 90°, ausmacht. Das Lager y = 0,7 nahert sich am meisten 
der Tragkraftkurve des unterteilten Lagers. Es stellt unter den zylindrischen Lagern etwa das 
Optimum dar. Fir ein Breitenverhaltnis o < 0,3 wird aber auch dieses Lager vom unterteilten 
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Lager iibertroffen. Vergleicht man nun das optimale zylindrische Halblager y = 0,7 mit dem opti- 
malen mehrflachigen Halblager, so ergibt sich, daB die Uberlegenheit des unterteilten Lagers etwa 
bei dem Breitenverhaltnis o = 0,3 beginnt und mit kleineren Breitenverhaltnissen Weiter anstelgt 
Die méglichen Tragkraftsteigerungen durch Einsatz von unterteilten Lagern betragen bei dem 


oO = ] a 0 3 25 4 9 
° 0, 0, 0 0, l 5 


Steigerung = 2% MIG 2896 7400 


Aus diesen Prozentzahlen geht hervor, daB der Einsatz von mehrflachigen Lagern in der Praxis 
wohl erst fiir Breitenverhaltnisse, die kleiner sind als 0,25 zu erwiagen ist. 
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gegeniiber dem zylindrischen Halblager mit x als Parameter. 


Abb. 12. Mégliche Steigerung des Flachendruckes | bei Anwendung des unterteilten Halblagers 


Von Shaw und Macks! wird die Ansicht vertreten, da es vom Standpunkt steigender Trag- 
fahigkeit bzw. abnehmenden Seitenflusses wiinschenswert sei, méglichst breite normale, d. h. nicht 
unterteilte Lager einzusetzen. Der Hinsatz breiter Lager stéBt jedoch auf Schwierigkeiten, da die 
auftretende Wellendurchbiegung das Lagerspiel zum Verschwinden bringt. Aus diesem Grunde 
empfiehlt es sich nach Shaw und Macks, schmale Lager zu verwenden und zur Herabsetzung des 
Seitenflusses, d. h. zur Erhéhung der Tragfahigkeit die Gleitflache in eine Anzahl von quadratischen 
Teilflachen aufzuteilen. Diese Ansicht, daB die im Lager auftretende Wellendurchbiegung den Ein- 
satz schmaler Lager erfordert, ist heute nicht mehr haltbar. Buske? hat in seinen Versuchen ein- 
deutig nachgewiesen, daB sich der Lagerwerkstoff den Durchbiegungen der Welle weitgehend 
anpaBt. So wurden z. B. im Priiflager (B/D = 1), obgleich der Biegepfeil bereits den vierfachen 
Wert der geringsten Schmierschichtstarke betrug, noch hydrodynamische Betriebsbedingungen 
nachgewiesen. Bei auftretenden Durchbiegungen dringt die Welle in den Bereich der Lagerschale 
ein, d. h. der Lagerwerkstoff gibt nach. Wie auch von Vogelpohl? formuliert, lassen sich die Wellen- 
durchbiegungen beherrschen. Dies gilt im allgemeinen auch bei gréBeren Lagerbreiten, wenn die 
von Buske erarbeiteten konstruktiven MaBnahmen Beriicksichtigung finden. 

ZusammengefaBt ergibt sich, daB nicht die Wellendurchbiegung im Lager als verantwortliches 
Kriterium fiir den Einsatz schmaler Lager anzusehen ist, sondern ihre Verwendung erfolgt nur bei 
besonderen konstruktiven Gegebenheiten, wie z. B. Raummangel. In diesem Fall ergibt sich eine 


1 M. GC. Shaw u. E. F. Macks, Analysis and Lubrication of Bearings, New York 1949. 

2 4. Buske, Stahl u. Eisen, 71 (1951), S. 1420. 

3 G. Vogelpohl, Uber die Tragfahigkeit von Gleitlagern und ihre Berechnung, Forschungsbericht Nr. 263 des 
Ministeriums f. Wirtschaft u. Verkehr des Landes Nordhr.-Westfalen, Koln u. Opladen 1956. 
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Uberlegenheit der mehrflachigen Lager hinsichtlich der Tragkraft nur fiir Breitenyerhaltnisse 
B/D < 0,25. 

Von Boyd und Raimondi! wurden mehrflachige Gleitlager nach Bauart 2 mit 4, 6, 8 und 12 Gleit- 
kérpern untersucht. Das dort gewonnene Ergebnis, daB eine Uberlegenheit der mehrflachigen 
Lager bei einem Breitenverhaltnis von B/D < 0,3 beginnt, stimmt etwa mit dem hier gefundenen 
Wert iberein. 


6. Die Berechnung von Gleitlagern geringer Breite. Zur Aufstellung von Unterlagen fir die 
Berechnung schmaler Gleitlager sind in Abb. 13 die Werte So (1 — x) tiber y mit dem Breitenver- 
haltnis o als Parameter aufgetragen. Dabei wurde die Erkenntnis von Kochanowsky’, da8B der Grenz- 
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Abb. 13. Auftragung der Kennzahl So (1 — x) tiber x mit dem Breitenverhaltnis o als Parameter. 


wert So (1 — x) fiir y—> 1 und alle Breitenverhaltnisse dem Wert 1,224 zustrebt, beriicksichtigt. 
Als Grundlage fiir diese Auftragung dienten die von Sassenfeld und Walther berechneten Sommer- 
eld-Zahlen. Abb. 14 enthalt aus den gleichen Unterlagen die hier verwendete Druckkennzahl 
Pm Ho 
nUD 

Die Ergebnisse der in einer neueren amerikanischen Arbeit von Kreisle* experimentell unter- 
suchten Lager der Breitenverhaltnisse o = 0,502; 0,334; 0,202; 0,1008; 0,0412 und 0,0187 sind in den 
Abbildungen 13 und 14 den rechnerischen Werten gegeniibergestellt. Es besteht, wie die Auftra- 
gungen zeigen, eine recht gute Ubereinstimmung. Abb. 13 und 14 bilden die Grundlage fiir die 
Berechnung schmaler Gleitlager, deren Gang abschlieBend kurz skizziert werden soll. 

Ein Gleitlager arbeitet dann betriebssicher, wenn die Belastung kleiner ist als die Belastungs- 
grenze P,,.n,, bei der das Lager den Zustand der Mischreibung erreicht. Sind das relative Lager- 
spiel y und die Gesamtrauhigkeitshéhe von Lagerschale und Welle h,,, bekannt, so lat sich daraus 
die relative Exzentrizitat 7,,.,,. bei der das Lager in die Mischreibung geht, berechnen. - Es gilt 


2 Fey 5 
Xgrens — Noe “yD : (42) 


iiber y mit o als Parameter. 


Mit Yrenz kann bei dem jeweiligen Breitenverhaltnis o in Abb. 14 die zugehérige Druckkennzahl der 
Pm Hj 


Tragfahigkeitsgrenze il 3 = konst. abgelesen werden. Hieraus lat sich dann die Grenztrag- 
fahigkeit zu a 
konst D® Bx 
P renz = 43 
f 60 heul ( ) 


1 J. Boyd u. A. A, Raimondi, Mechanical Engineering 75 (1953), S. 380. 
* W. Kochanowsky, Schmiertechnik 2 (1955), H. 4. 
8 L. F. Kreisle, Transactions of ASME (1956) S. 955. 
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ee Darin ist 1) die im Lager auftretende mittlere Viskositat des Schmiermittels und n die 
= ince in ae ae es der Temperatur abhangt, ist die Kenntnis der sich einstellenden 
gertemperatur er ich. Sie la i ar i S i i 
a : ener ae et erlic ie laBt sich aus der Warmebilanz des Lagers ermitteln. Es ergibt 
: On 2 UPR fe 
9 —9,=— =" PUB. (44) 
In dieser Gleichung hat « die Bedeutung der Warmeiibergangszahl und A die der warmeabgebenden 
AuBenflache des Lagers. Da 7 eine Funktion der Temperatur @ ist, muB zur Lésung der Gleichung 
(44) das Viskositats-Temperaturverhalten des verwendeten Schmiermittels bekannt sein. 
Beispiel. Als Beispiel aus der Praxis fiir die Anwendung schmaler Gleitlager sei hier die Kurbel- 
wellenlagerung einer amerikanischen Pipeline-Pumpe erwahnt: Es handelt sich um eine 5-zylindrige 


Kolbenpumpe Quintuplex VQE 2, 3,5” x 5”, 


der amerikanischen Worthington-Gesell- Pry 
schaft. Das als Beispiel gewahlte Lager hat gues} 
einen Durchmesser von D = 9!/,’” + 241 mm, 9.020 ce or oh) 
eine Breite B + 1°/,,/’ + 48,2 mm und da- re . 
mit ein Breitenverhaltnis von o = 0,2. Die 0015 
maximale Lagerbelastung ergibt sich bei 
der Nenndrehzahl n = 360 U/min unter 0070 gu gg. ha 
Beriicksichtigung der Massenkrafte zu P ast - 
= 4120 kg. Bei der Férderung einer Fliissig- On - Me s = 
keit; die eine maximale Temperatur von 9.006 L 
150 °C haben darf, kommt es zu einem 0005 4 tal 
WarmefluB zu den Lagern. Es kann aber i ‘6 | 
angenommen werden, da’ die Umlauf- aa 
schmierung in der Lage ist, die Lager- 0003 “08 
temperatur auf 60° zu halten. Mit Riick- ? aoe 
sicht auf die Erwérmung wird ein relativ ie 
groBes Lagerspiel yp = 2-107? angenommen. By 
Die Rauhigkeitshéhe wird zu h,,, = 10 u outs as! | 
geschatzt. Mit diesen Werten ergibt sich 
mach (42) Ypgrenz = 0,96. Das heiBt, bei 001 a ae Sl eed 
x = 0,96 geht das Lager in die Mischrei- 0009 | 
bung. Bei einem Breitenverhadltnis von aa |__Versuchswerte vonircisle 
o = 0,2 betraigt die Druckkennzahl nach 40006 | “Oa te 
Abb. 14 0005 cae o7008 
Ay Ma Ze =| = 0,006. go0o4 Beals govte, 
WO pranks.) UD ee i a relative 
Die Grenztragkraft ergibt sich dann nach SU “o7 4203. OF OS 06 G7 08 et 
g G2 038 OD 6 OF 08 08. 4 
(43) zu ro ge 
ae = 2120 Hn. Abb. 14. Auftragung der Druckkennzahl = iiber x 
Mit einem Ol 90° E/50° (Se das bei 60 10; mit dem Breitenverhdltnis o als Parameter. 


eine Zahigkeit von 7 = 80 c P hat, ergibt 
sich die Drehzahl, bei der das Lager in die Mischreibung geht, zu n = 243 U/min. Das Lager 
arbeitet demnach betriebssicher im Zustand der FlieBreibung. 

Berechnet man in dhnlicher Weise die Tragfahigkeiten schmaler Gleitlager und vergleicht sie 
mit denen von Walzlagern gleicher Breite, so stellt man fest, dafs von einer allgemeinen Uber- 
legenheit des ,,Kurzwalzlagers“ gegeniiber dem Gleitlager nicht gesprochen werden kann. Schmale 
Gleitlager sind durchaus in der Lage, bei gleicher Lebensdauer gleiche oder sogar héhere Krafte 
als das Walzlager aufzunehmen. Dabei muB allerdings besonders bei héheren Drehzahlen die 
eréBere Reibungswarme des Gleitlagers durch geeignete KithlmaBnahmen abgefiihrt werden. 


7. Zusammenfassung. Die vorliegende Untersuchung hat gezeigt, daB im Falle des unendlich 
breiten Schmierkeiles Flachenunterteilungen zum Zwecke der Tragkrafterhéhung in der Praxis 
nur dann Sinn haben, wenn die Keilsteigung ¢ sehr groBe Werte annimmt. Im Bereich der Keil- 
steigung 1 < L ¢/H, < 10 wiirde eine Unterteilung zu einer zum Teil wesentlichen Tragkraftver- 
minderung fiihren. Da sich in der Praxis Keilsteigungen, die in den Bereich 1 < L ¢/H, < 10 fallen, 

15 
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fast immer herstellen lassen, hat die Untersuchung des Bereiches L ¢/H, > 10 nur theoretischen 
Wert. Sie zeigt jedoch, daB hier eine Unterteilung sinnvoll sein und zu héheren Tragkraften fiihren 
kann. Es lassen sich, wie in Abb. 6 dargestellt, die Bereiche angeben, in denen die jeweilige Teil- 
flachenzahl n die héchste Tragkraft liefert. 

Aus der Untersuchung des endlich breiten Gleitschuhes folgt die Erkenntnis, daB die optimale 
Tragkraft nur dann erreicht ist, sobald ein Breitenverhaltnis B/L = 1 und eine Keilsteigung 
9 = 0,8 eingehalten wird. Die gréBten Tragkrafte bei Schmierkeilen von endlicher Breite werden 
immer bei einer Unterteilung in quadratische Teilflachen erreicht. Jedoch liefern auch die benach- 
barten Werte von B/L = 1 und g = 0,8 hohe Tragkrafte. 

Der Vergleich des mehrflachigen Halblagers mit dem zylindrischen Halblager bestatigt ein- 
deutig, da® fiir kleine Lagerbreiten das unterteilte Lager dem zylindrischen Lager hinsichtlich der 
Tragkraft iiberlegen ist. In der Praxis diirfte sich jedoch der Einsatz mehrflachiger Lager erst fiir 
Breitenverhaltnisse o < 0,25 lohnen. Dabei betragt die mégliche Tragkraftsteigerung 17%. Sie 
steigt fiir ein Breitenverhaltnis von o = 0,15 auf 74% an. 

Die Auftragungen in Abb. 13 und 14, die in guter Ubereinstimmung mit Versuchswerten von 
Kreisle stehen, gestatten die Berechnung sehr schmaler Gleitlager. Eine Gegeniiberstellung von 
Gleit- und Walzlagern zeigt, daB die Gleitlager geeignet sind, an Stellen wirklichen Raummangels 
Walzlager hinsichtlich der Tragkraft vollwertig zu ersetzen. 


(Eingegangen am 7. August 1959.) 
Anschrift des Verfassers: Dr.-Ing. Heinz Peeken, Kiel, Hasseldieksdammer Weg 27. 


Aus dem Max-Plank-Institut fiir Stromungsforschung, Abt. Reibungsforschung 
(Leitung Prof. Dr.-Ing. G. Vogelpohl). 
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Analysis of Plane Stress in Cartesian Coordinates 
and with Varying Thickness 


By H. D. Conway 


1. Scope of the Investigation. The biharmonic stress-function equation for plane stress in 
cartesian coordinates x, y is generalized to take account of varying thickness perpendicular to the 
plane of the stress, which thickness is a function of the coordinate y. For stress functions in the 
important forms Fy) sin y x and F(y) cos y x, this equation is found to be particularly tractable if 
the thickness is proportional to e’, since it leads to an ordinary, homogeneous, fourth-order differen- 
tial equation in F(y), the characteristic equation of which is, moreover, easily factored. If ky is 
everywhere small, a linear variation of thickness can then be approximated to. A numerical 
example is given for a square-plate problem, the stresses being compared with those from the 
corresponding constant-thickness solution due to Girkmann!. 


2. Analysis. If it is assumed that the thickness, perpendicular to the plane of the stress, is 
varying but everywhere small compared with the other dimensions of the plate, the usual method 
of analysis employing stress-functions may be used. We proceed in the usual manner but with this 
generalization. 


normal stress—p sin %x/1 


ae 
hoy+ 9 (Roy)dy 

try tee (toy) dy 

Itey +a (Rtey) det 
hOp ho,+32(ho,)dx 
SIG a fr) thickness h=h (2) yf 
Wry 
y 
Nigy 
LZ hoy : 


Fig. 1. Forces per unit length on sides of element. Fig. 2. Loaded Square Plate. 


Fig. 1 shows a rectangular element of thickness h = f(y), the equations for horizontal and 
vertical equilibrium being 


3 a 8 ee re 
2 hod +Z hte; =0, — Flhay) + 5 (hts) =O. (1) 
It is seen that these equations are satisfied by a stress-function P(x, y) such that 
aD aD PD 
ace = =— ‘ 2 
Paes pas? bx by @) 


The stress-strain equations are, with the usual notation 
ou 


6v : 
eye hee (3) 


éu ov 
(ay + ae) — tr 


1 K. Girkmann, Flachentragwerke, 4th edition, page 53, Wien 1956. 
ag Tis 
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and, by substituting Equations (2) in (3) and eliminating the displacements, the compatibility 


equation is 
cf / 2 22 2 92 D2. e? 1 2D 
CP MO ene. one ! é ACT eae es ee ae oy) ae (4) 
6x? \h_ Ox? h éy? dy*\ h dy? h 6x? x dy\ h ax dy 


A thickness h is now selected in the form 
hi hie (5) 


where h, and k are constants, which, after substitution, gives 


a@ aD oD aD 
4 racy Pe *, 2 LS | I ne 
A’@ — 2k ays +k ay yk aa 2k ax? dy 
where (6) 
aD a+®D ot®D 
1h 2 
A‘? oxt te Ox? Gy? dy? 


The stipulation that h is everywhere small naturally limits the magnitude of k y. It will be noted 
that, if k is zero and consequently that the thickness h is constant, Equation (6) becomes the well- 
known biharmonic equation for the stress-function. If k y is everywhere sufficiently small for higher 
powers to be ignored, the thickness variation is approximately given by the linear equation 


h = hg (1 + ky). 
For loadings which are discontinuous, stress-functions in the form of Fourier series or integrals 
are commonly used, and accordingly we consider stress-functions in the form 
m 7 


P= Fly)sinyx, F(y) cosy x with y = —— (7) 


where m is an integer and F(y) a function of y only. Substituting (7) in (6) then results in the follow- 
ing differential equation for F(y): 


FV 9k FM + (4222) FU +2 hy? FI + (yt +) ky) F=0 (8) 


where the superscripts denote differentiation with respect to y. For solutions of the form F = A e¥, 
the characteristic equation is 


M—2k#B + (M—2y)24+2kyA4+ (pt tvky?) =0. (9) 
Very fortunately, this quartic can be factored without much difficulty. We write, tentatively 
(hat fy) +a—Bf)(A+0 +f) +a—f,) =0 (10) 
where «a = —k/2. Multiplication then gives 
3k? Be 2 
M22 + ("5 — A — Ba)? + (K+ HG) + (F—6i)(7 A) = 0. ) 


Comparing the respective coefficients of 7? and A in Equations (9) and (11) results in the same 
equation, viz 


k2 ; 
+ mB poy, (12) 
whereas a comparison of the terms independent of J yields, with Equation (12) 
xt k4 k2 y? 
ee DL erg ah yore ea (13) 
From (12) and (13) it may be shown that 


Kk? age ts 
B=Gtrthyhi, pak + ky pri (14) 
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where i = /— 1. Equation (14) with a = —k/2 gives the roots of Equation (9). Thus F(y) is 
known and the stress-function may finally be written in the form 


@D = ekyl? [e—4y (C, cos By + Cy sin By) + e47 (C, cos By +(C,sin By)]sinyx, cosyx (15) 
where C,, C,, C, and C, are constants of integration, and 


(ek k2  y2\2 22 1/2) 1/2 = k2 2 ke p22 2 2 y]1/2) 1/2 
eS ree r (s+'s) + 4 | | : LS | ocr (e 2] 4 ae) 


If y is zero, Equation (9) has multiple roots, giving 


k2 12 
see ”) bf 


PD — esr a 


kt 
4 


\ 12 
cee ”) ene 


siny x, COSy Xx. (17) 


3. Example. As an example, the numerical solution of the problem shown in Fig. 2 of a square 
plate loaded on its top face by a normal compressive force distributed so as to give a stress Oy 
there of — p sin z x/l, and supported by shearing forces on its vertical faces. The thickness varies 
from h, at the bottom to 2 h, at the top, whence the law of thickness variation is h = hy (2)¥!!. The 
constant-thickness version of this problem has been analyzed by Girkmann. 


The appropriate stress-function is given by Equation (15) with y = 7/I, the boundary conditions 
necessary for the determination of C,, C,, C,; and C, being 


y=0: op =t,=9, 


yas (18) 


y=l: of =—psin i > Tey = 0. 


The boundary conditions of zero normal stress on x = 0 and x = | are automatically satisfied by 
the factor sin z x /] in the stress-function. The stresses from Equations (2) and (15) are then 


—ky/2 ; 
Gees oa A ee B P| [CO euiiay oC cin By) 
0 


—2B1l 


2] A IC, sin B y — C, cos By) — atr(((F + A eee P| [C,cos By 


+ C,sin By] —2 BIlT) + AI][Gysin By — Cy cos By)}} sin, 


(19) 
o me Bay (C, cos By + C,sin By) + e4¥ (Cz cos By + Cy sin B y)] sin “*, 
y ho? 
t et fey [B1(— C,sin By + C,cos By) +(BI—A 1) (C, cos By + C,sin By)] 
xy hE? 


+ e4y [B1(— C,sin By + C, cos By) + (B1+ Al) (C, cos By + C,sin By)]} cos. 


Assuming » = 0,25 and kl = In 2 = 0,69315 for the required thickness variation, A and B 
are obtained from Equations (16) and the constants of integration obtained by solving simultane- 
ously the four equations obtained from the boundary conditions in Equation (18). Hence 


Al = 3.16533, Bl = 0.17199, 
C, =— Cy = — 0.0291808 ph, 2”, 
C, = — 0.9433000 phy 2, 


C, = —0.1307951 phy I?. 
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Fig. 3. Graphs of on. 1/2" 
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Fig. 4, Graphs of o. 7 
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Numerical values of Ox, — 1/29 Oy, 12 and Txyy — 9,1 ate plotted in Figs. 3, 4 and 5, respectively, 


where they are compared with the corresponding values for a square plate of constant thickness 
obtained by Girkmann. The values of the stresses for the latter case are, of course, independent of 
the value of y. 


(Eingegangen am 13. August 1959.) 


Anschrift des Verfassers: Prof. H. D. Conway, Thurston Hall, Cornell University, Ithaca, New York, USA 
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Die Singularititen der frei drehbar gelagerten orthotropen 
Halbplattenstreifen 


Von G. Nitsiotas 


1. Einleitung. Die Singularitaten der unendlich ausgedehnten Voll- und Halbplatten wurden 
eingehend von Stein! behandelt. In der vorliegenden Arbeit sollen unter Anwendung des von 
Pucher? entwickelten Verfahrens fiir die isotropen Plattenstreifen geschlossene Formeln fiir die 
Singularitaten der Kriimmungen und der Verwindung von frei drehbar orthotropen Voll- und Halb- 
streifen abgeleitet werden. Den Untersuchungen wird die Hubersche Differentialgleichung 

o4w 64w , Ow 


Bs axt Eo oy? i dy* 


0 (1) 


zugrunde gelegt. Hierin sind K,, K, die Biegesteifigkeiten der Platte in der x- und y-Richtung. 
Bezeichnet man mit (4,, 4y die bezogenen Poissonschen Zahlen in denselben Richtungen, so ist die 
Drillsteifigkeit des vorliegenden Systems 


ap 2C =(1—/u, uy) VK. K, ; 
somit ist 
2H = (Kyu, + Kyu.) +4C€. 


Das Belastungsglied p(x, y) wurde in (1) gleich Null gesetzt, da wir uns auf die Untersuchung von 
Einzelangriffen beschranken. 


Die allgemeinen Gleichungen fiir die Biegemomente und das Drillmoment lauten 


O4w 6?w Ow O?w ew 
ae | syn one 
m, = — K, (Ga fe hy rd) : my = Ky (2 [by 73) ; My =—2C fee 


Andererseits gelten fiir die bezogenen Querkrafte die Beziehungen 


La 


?w , 2C\ aw ae K 6? w | ,2cG\" 6% 
el tele Ga Xl t+ hem 


/ 


q, =—K 


x 


Mit Riicksicht auf die Aufstellung der Randbedingungen der Halbstreifen mit freiem Querrand 
schreiben wir die fiir die statisch gleichwertigen Ersatzquerkrafte geltenden Formeln auf: 


eek. K 6? w _4C 8 w poe) K 6? w es 02w 
Ie ~ x | 38 + (Hs | Ke) maa Ee » lays ae K) ay 


Zur Integration der homogenen Differentialgleichung (1) erweist sich als zweckmafig die Be- 
nutzung des Bernoullischen Ansatzes? 


28) 5 
w=Ae ° sin “~ x, (2) 
welcher die Randbedingungen 
0?w 
ey und aye (m, = 0) 
entlang der frei drehbaren Rander 
a= 0 und x=a 


des unendlichen Plattenstreifens von Haus aus erfiullt. 


1 P, Stein, Dissertation TU Berlin, Stahlbau-Verlag 1959. 
2 4. Pucher, Ing.-Arch. 14 (1943), S. 259. 
8 Grundsatzliches iiber orthotrope Platten bei K. Girkmann, Flachentragwerke, S. 300. Wien 1956. 
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Nach Einsetzen des Ausdruckes (2) in (1) ergibt sich zur Bestimmung der freien, noch wahl- 
baren Konstante A die charakteristische Gleichung 


K,—2H#2+K,iM=0. 


Fihrt man nun in dieser den Kennwert 


— tei fie 
ak 
und das Biegesteifigkeitsverhaltnis 
B = |/Ke 
/ Ky 


ein, so bekommt man folgende vier Wurzeln: 


ha =p Ve + Vet = +6 (V/A + Vs). (3) 


Im Hinblick auf diese Werte, haben wir die nachstehenden drei Falle zu unterscheiden. 


Erster Fall: Uberdrillsteife Platte 0 +x > 1. Alle vier Wurzeln sind reell. Davon 
scheiden die negativen A; , aus, da_w fiir y = 00 nicht unendlich werden kann. Somit verbleiben 
die positiven Wurzeln /,, A, welche durch die Formeln 


Ac ae =p) 
gegeben sind. 


Zweiter Fall: Isodrillsteife Platte x = 1. Fiir diesen Wert liefert (3) zwei Doppelwurzeln 
+ B. Im allgemeinen ist dieser Fall gleichwertig mit dem Fall der isotropen Platte. In der Tat, 
durch die Koordinatentransformation 


R= und y = py 


verwandelt sich, unter Beibehaltung der Randbedingungen, die Differentialgleichung (1) in die 
biharmonische Gleichung 


Dritter Fall: Unterdrillsteife Platte 1+>x>0. In diesem Falle sind jeweils zwei 
Wurzeln konjugiert komplex, davon kommen nur die mit positivem reellen Teil in Frage: 


A=optiy, Ag =~9—tY, 
ia Sy 
e=b/ 3". v= pV 


2. Die Singularitaten des unendlichen Plattenstreifens. Betrachten wir nun den frei aufliegenden 
Plattenstreifen mit der Einzellast P = 1 im Aufpunkt A(u, v). Die Belastung p(x) entlang der 
Geraden y = v stellen wir mittels einer Fourierschen Reihe als ungerade Funktion von x mit der 
Periode 2 a dar. Die Entwicklung lautet 


(2 a= z sin y sin x 
PO eg >) ny ae 


Fiir die Biegeflache wahlen wir den Ansatz 


ba A F “Za (y—») A ¥ “7a (y—») % nz 3 nw 
w on ye rene Fe. (4) 
welcher die Differentialgleichung (1) und die Randbedingungen entlang der Geraden x = 0,x =a 
bereits erfiillt. Das negative Vorzeichen im Exponent gilt fiir den Bereich y >v. Der Index IV 
bezeichnet den Fall des frei aufliegenden Plattenstreifens. 


Der Ansatz (4) stellt die gesuchte Singularitat dar, wenn entlang der Bereichgrenze y = v die 
Randbedingungen 


wobei 


ist. 


Ow 1 
(35), Meme ee ie 


XXIX. Band 1960 G. Nitsiotas: Die Singularitaten der frei drehbar gelagerten orthotropen Halbplattenstreifen 225 


erfiillt sind. Beriicksichtigt man die Symmetrie um die Gerade y = v , so bekommt man 


Paes = 0) Ay,v 
dx" Oy a : | 


Mithin transformiert sich die zweite Randbedingung zu 


Sw 1 “A 
Sere a 
Yao phe 
| bu 
Hiermit ergeben sich zur Bestimmung der zwei Unbekannten A,, Ap ye od 
? u 
die linearen Gleichungen 
A, — A, =0, 
n 7 \3 1 
—K,("F) (—Aj A, Ag 42) =. ; 
pucieeas 
Beriicksichtigt man nun die Gleichung Lo : 
c Q@—-——+ 
2 Abb, 1. Plattenstreifen 
i — Ay = 2 B? \e—1 Be 2 Vel ? iti) Finuallast seduce (u, v). 
und setzt man K = VK KS K, fiir die mittlere Biegesteifigkeit der orthotropen Platte, so bekommt man 
=/1 1 SSS 
Agama == — : Ke=).K, K 
1 2 eK: 2 a | nw 3 ( / 4) 
a 
Diese Werte in (4) eingesetzt liefern die Singularitat der Durchbiegung des Plattenstreifens 
nam nz 
- 1 1 ie AO) Ve so) nz ni 
TV nok A ee 4 in —— in 
S a oH mo ob Ae sin —— usin —— x. (5) 
a 


Durch zweimalige Differentiation nach den Koordinaten u bzw. v des Aufpunktes A ergeben sich 
die den Kriimmungen und der ashlng zugeordneten modifizierten Singularitaten. Es wird 


Sis ss tee pr ey AG) a ng 
> ——e sin — usin —%*, 
du* Seren A, a a 
aglV = AG») FA aly —») 
= = > Apes ; tA,e = sin “* wsin“~ x, (6) 
mt Cae 


iia ete ¥ “= Aly —») SF “3 Aly — v) nw . NTA 
5 == cos —— UW sin ——_ %. 
ou ov esi a a 


Diese Ausdriicke kénnen auch in geschlossener Form dargestellt werden, wenn man die Identitaten 


nm 
1 Tatty 9) ‘ ; a 
I(a,y—v,x—u)=4 ye sin ““ usin “™ x 


Coj = A(y —v)— cos — (& + u) 
=n: ») 
Coj = A(y — v) — cos — (x — u) 


(7) 


; FAT Aly — 9) 
G(,y—v,x—u) =2 De ¥ cos "™ wsin “= x 
sin — (x + u) sin — (x — u) 
= arctg = + arctg 7 a 
oe a Sree ckeg Pau : — cos — (x —u) 
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in Betracht zieht. Diese Formeln gelten sowohl fiir reelle, als auch fiir komplexe Werte von A, wenn 


$F ay —)| <1] 


| 
| 
|e 
ist. Da der reelle Teil der Wurzeln 4 immer positiv ist, wird obige Konvergenzbedingung erfiillt. 
Nach Einsetzen der Werte (7) in (6) erhalt man die geschlossenen Formeln 


ese 1 Te sok: 

du* 8a KYe—1 Ay Ay}? 

e2sl¥ i 

Si ae eee (A td Ys 8 
ove aK Veni | ify atl 2 2) (8) 


opus 42 


dudv 8 aK V2—1 


wobei I,(G,) und I,(G,) die Werte der Ausdriicke I(G) fiir A = A, und / = A, sind. 

In dieser Form kénnen die Singularitaten nur fiir die Berechnungen im Bereich co > x > 1 
beniitzt werden. Im Falle von unterdrillsteifen Platten, (komplexe Wurzeln A) erweist sich als 
zweckmaBig die Umbildung der Formeln (8). 

Fuhrt man die Abkirzungen 


A = oj = p(y —») cos — p(y — 0) — cos — (x +u), 
B =Goj = oy —») cos — p(y — ») — cos = (x—u), 


C = Gin = y(y—») sin = p(y — 0) 


ein, so bekommt man 


g2gl¥ ee! y2 1 - | A? + C2 1 (A— B)C 
bu2 82K Bp yl +x is 4.G2"3 jie are tg AB AC? |? 
Lies a 1 A? + C2 1 6 PAZ BG 
aa CD) ee BANG Ets See a 
dv 82K \ yi+tx BHO et Sa © AB+ C?]° 


Die Singularitat der Verwindung 1aBt sich auch schreiben in der Form 


Ren Fae st cos [p(y —») — (x + w)] —6oj = oly —v) 
du év 8a KYl—# 


n 
cos [p(y —v) + + u)] — Goi = p(y—») 


4 
en [vy —v) — (&—u)] —Coj = p(y —2) 


cos [v(y — v) + (x—u)] — Gof = oy —2) 


Um die Singularitat fiir x = 1 baw. 4, =/, =f au erhalten, muB der Grenziibergang lim ge- 


“ol 


bildet werden. Allerdings auch fiir Werte von x in der Nahe der Einheit (« = 1+ Ax) erweisen 
sich die Formeln (8) ungeeignet zur Berechnung, da das Resultat sich als Differenz zweier fast 
gleicher Zahlen ergibt. Aus diesem Grunde erscheint es zweckmaBig, die Ausdriicke (8) in Reihen 
nach den Potenzen von AA = 4, —/, zu entwickeln. Wir schreiben zur Abkirzung 


A, =B +A, Ay =B —Ah, 


Ah = hy Iq = Ady + Ady = 26 Y=, 


Ah — Aly =}, +A — 28 =28()/* 4) —1) 4 
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und entwickeln [,//, in Reihen. Es wird 


f _ FB+ 44) _ (B— AA) (2 +44) jot ee " 

Ay B+ Aa, pt Ip I Ad, Ip+ > AA Ip +..| 
see f ea: Any Has Ig | Boon 
or [Te + 4a (1h | 3 ( Tp 4 palics Ii). 


wobei Tp, Tz, Ip die Werte des Ausdruckes fiir 1 =f und seine Ableitungen nach f sind. Ganz 
ahnlich la8t sich der Bruch I,//, entwickeln. Daraus und unter Bericksichtigung von (9) folgt: 


e2glv 1 y Ig eo De epee u 
ef — 8a K (15 2) 3 Vi al 2) | ( V5 ial B Is). (103 
Por = 1 (A= 6) bekommt man die Singularitat der isodrillsteifen Platte zu 
Geri io.2 begat 
Se Eee Se ee 


Diese Formel stimmt fiir 6 = 1 mit der von Pucher! fiir die isotrope Platte abgeleitete Formel iiber- 
ein. Durch Einsetzen des Wertes (10b) in (10a) ergibt sich 


sagt 2 ese” . BIg Beas 
ou? 26/5 2) éu2 | a V-n)- (11a) 
In gleicher Weise erhalt man fiir die zwei anderen Singularitaten folgende entwickelten Aus- 
driicke : 
2G1V highs uarae! IV 2 ig u eye 
os myo wei Ey, any 
Gv x+1 = év? 8aK l+x 
eesi¥ aS Beg = e255" 9 paige”, / 2 
meee Viet 1 ba bo oe #(1— et? ; ae) 
wobei 
255. 1 vent OS Rive: 
av? = igh tas du dv tae (i?) 


ist. Dies sind Singularitaten der isodrillsteifen Platte. 

Die Ausdriicke (8) stellen modifizierte Singularitat dar, welche neben den klassischen Singu- 
laritaten S der unendlich ausgedehnten Platte eine regulare Funktion erhalten, deren GréBe im 
Aufpunkt bestimmt werden soll. Im allgemeinen gilt 


Reg D[S!V] = D[S!”] — D[S] . (13) 


D{_ ] bedeutet eine Differentialoperation nach den Koordinaten des Aufpunktes A(u, v). 

Wir gehen von der klassischen Singularitat von Stein? aus, da bereits fiir diese die Berechnung 
des Restvolumens des in den Héhenlinienplanen nicht mehr darstellbaren Teiles des unendlich 
langen Schlauches im Aufpunkt durchgefiihrt wurde. 

Setzt man in der Formel (13) einerseits die Werte (8) der Reihe nach, andererseits die von Stein 
ermittelten klassischen Singularitaten ein und entwickelt man die trigonometrischen und Hyperbel- 
funktionen in Reihen nach den Potenzen von y — v und x — wu, so lassen sich bestimmte Grenz- 
werte leicht nachweisen. Demnach wird im Aufpunkt 

I 4 
a 


e2.SlV 1 if Z z ( 2 
Ree (*)=— 5244 as Top eae 


eet Ine +11), (14) 


p2gh¥ 1 2 2 au 2 alle 
Reg ( = te Vena me (1 — cos 22*) + 6 V3 (3 in al 


Da die Singularitat der Verwindung im gesamten Gebiet endlich bleibt, kann sie durch die Hohen- 
linien vollstandig dargestellt und ausgewertet werden. 


1 Siehe FuBnote 2 von Seite 223. 
2 Siehe FuBnote 1 von Seite 223. 
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Aus der Singularitat des Plattenstreifens kénnen auch die Singularitaten der unendlichen Platte 
abgeleitet werden. Es geniigt in (8) die trigonometrischen und Hyperbelfunktionen in Reihen 
nach den Potenzen von 1/a zu entwickeln und den Grenziibergang fiir a—> oo zu bilden. Hiermit 
ergeben sich bis auf eine belanglose Konstante die von Stein aufgestellten klassischen Singularitaten. 


3. Halbstreifen mit frei drehbarem Querrand. Der Spannungszustand des Halbstreifens ist 
derselbe wie fiir den Bereich y > 0 des vollen Streifens, der in Punkte (u, + v) die Last + 1 und im 
Punkte (wu, — v) die Last — 1 tragt. Wegen der Antimetrie verschwinden entlang der x-Achse die 
Durchbiegung w und die Kriimmung 0?w/dy?. Nach dem Uberlagerungsgesetz bekommt man mit 
Hilfe von (8) folgende Singularitaten: 


2oV 

Oe = Fe eI BO + 7 ho — 9) — By + ol, 
Sei ay [I I de [I 

at ~ snk Yen | 1 [hy —*) — Ey + »)] +42 [hy —») — Ly +») ]}. (15) 
ozs 2 


Fa Bs = 5K PFET 9 — GO $91 + 1G 2) — Gly + | 
Der Index V bezeichnet den Fall des frei drehbaren Querrandes eines Halbstreifens. Fiir den Fall 
komplexer Wurzeln A lassen sich obige Ausdriicke entsprechend umbilden. Andererseits gelten im 
Falle der quasi-isodrillsteifen Platte («~ 1) ahnliche Gleichungen wie (14). Diese ergeben sich, 
wenn in (11) der Index IV durch den V ersetzt wird und im Restglied Ig(y — v) — Ig(y +») 
statt Iz eingesetzt wird. Die Singularitat 0*S¥/du? stimmt fiir 6 = 1 mit der Singularitat der 
isotropen Platte tberein. 

Wie beim vollen Streifen, kénnen auch hier die regularen Teile abgeschieden werden, welche 


; | eee a ae 1 2 
b f eat £5 : z 5 ; 
is auf das Glied — 5 F / 71 Aund —,—> 8 \; TT A mit den Ausdriicken (14) tiber 
einstimmen. Es ist 
Pera 2 
oj “a — cos = 2 
Tea — (16) 
Cofj —1 
a 
Wye 
| 
A) 
+7 te 
| 
ID 
u me 
| | | 
sean 
| aes 
h-7 
| | 
| | 
: bee 
Abb. 2. Halbstreifen mit frei drehbarem Querrand. Abb. 3. Halbstreifen mit eingespanntem Querrand. 


: 4. Halbstreifen mit eingespanntem Querrand. Die Singularitat SV! der Durchbiegung setzt 
sich aus der Singularitaét S!V des vollen Streifens und einem regularen Integral w der Differential- 
gleichung (1) zusammen: 


SvI= SV 4, (17) 


he S'V bereits die Randbedingungen entlang der Geraden x = 0 und x = a erfillt, kénnen wir den 
nsatz 
ae i 1 A SAG ea 
Pre ean a 2 a - nae - n 
serie pl Tee aie oe “ae (18) 


a 
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wahlen. Langs des eingespannten Randes y = 0 verschwinden die Durchbiegungen und die Nei- 
gungen: 


ay Ow Sly 
YN acs 
w Be ober les Flue 


Setzt man nun in diesen die Werte (5) und (18) ein, so bekommt man zur Bestimmung der Unbe- 
kannten A,, A, die Gleichungen 


| A l —*2Ae ea det 
=a ate! 2 =e is y == ie 
A A A hy 
ee he ee eee 
A,—A,—e “ +e ° =—() 
Hiermit ergibt sich 
2 Bay | aD 27 ——A,v 
ef ees St ee es BAL 
pea ey She ieeaag 
3 an. ie ts hi sad, Pa 
a ds 1A, 


Der Regularteil w wird demnach 


Ss —*Za(y +») — "aly +) 
ae >> 1 jouer 1 Ae ike a 
2 aK y2—1 nm\3 Ay — Ag Ay Ay 


nim n 

9 ag Cid Ae) — AF (aay +40) nr na 

+ =—_ |e ee 7 sin — usin — x. 
Ay — Az a a 


Die héheren Singularitaten folgen durch Einsetzen dieses Ausdruckes in (17) und Differentiation 
nach den Koordinaten des Aufpunktes. Unter Beachtung von (7) kénnen diese in geschlossener 
Form angeschrieben werden, indem man die Summen durch die identischen geschlossenen Aus- 
driicke ersetzt. Es wird somit 


gis tos 1 je see 
du? 82a K)x2—1 Ay As 


Aaaies v I 2, 
phthy wt) 4 sn oe Ary + a0) + May + Ae) 


OSs: 1 
av? 8a Kk ye—1 


PERG AO +) +My +) + 7p My +h) + Bay +4, 


AS Ay L(y — 0) +g L(y —»)] | 
(19) 


Ay 
e281 ) 
gudv 8nKye—l1 | ee, 
ee na » [Gay +) +6,(y +v)] a A. 7, Ve G(Ayy Ag v) +A, 6 (day +A, | . 


Liegt der Aufpunkt auf dem eingespannten Rand (v = 0) so bleibt der Singularanteil des Ein- 
fluBfeldes des Stiitzmoments im gesamten Bereich einschlieBlich des Aufpunktes endlich und kann 
vollstandig dargestellt werden. Es wird somit die entartete, modifizierte Singularitat des Stiitz- 


momentes 


2gVt Rees. a 
sates Bee ee. (1, i) = shee (i 1). 
ou? \v=0 8a KYe—1 20K y2(% —1) 

Fiir kleine Werte x — 1 gilt naherungsweise die Formel 

e2s¥! ee Ste es ( x+1 

ov? Bi! av? p20 «820K 2 ‘i 
wobei 

ay" 1 


oS ae 
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die von Snidaric! ermittelte Singularitat des Stiitzmomentes des isotropen Halbplattenstreifens 
(Ba) ist: 
Zu der klassischen Singularitat 02S1/dv? des Stiitzmomentes gelangt man, wenn man den Koor- 
dinatenursprung in der Mitte des eingespannten Randes verlegt (d. h. x durch x — a/2 und u durch 
u — a/2 ersetzt) und die trigonometrischen und Hyperbelfunktionen nach Potenzen von l/a ent- 


wickelt. Es wird 


lf 2 2 1 biallearz 2 
Sirus 1 Se ee ee peels RA | Sree 
Cof —- Ary cos — (x + a) 1 (a4) +1 9 ery) a 
lim J, = lim In Wee MSE = lim In ae 5 I tx a 
a 00 esac” (ai) | cos a u) 1 9 (79) see 1+ aes (x | 


= lim In 45 —In [2 y? + («—u)]. 


Demnach folgt 
ae ES: {In [22 y® + (« — u)®] —In [23 y? + (2 — uy} 
ov? 29n K Y2(x—1) : 
= lage Re ee, 
a K y2(x—1) B? (x V2 — 1) y? (x He 


5. Halbstreifen mit freiem Querrand. Zur Vereinfachung der Berechnung wird der EinfluB der 
Querdrehung auf die Randschnittkrafte vernachlassigt. In diesem Falle lauten die Randbedingungen 


entlang der Geraden y = 0 


e2sVll) poe 2 A is 0 20 
( oy” j= eye + 2%8 Ox? dy /y =0 4 ( ) 

wobei 
Svit — SIV 449 (21) 


ist. Fur den regularen Teil wird der Ansatz (18) gewahlt. 
Setzt man nun die Werte (5) und (18) in (20) ein, so bekommt man zur Bestimmung der Un- 


bekannten A,, A, die Gleichungen 


nz nm 
—A,v — —A,0v 


AYU 4 ra et a 


BG eet es Avil + Ga ee 


Fir die zweite Randbedingung wurde die Identitat 2% 6? = 7? + 72 in 
Betracht gezogen. Durch Auflésung des linearen Systems erfolgt 


nz nw 
AEA gS DM A re 


Sih Ans ogee. A, = 6 + — @ 
. eas: RH 
Abb. 4. Halbstreifen ne aie 
mit freiem Querrand. te ese poe ese, 
ee pele a ay ge a 
2g aa Teeny ats 
emer R35 


Der Regularteil w wird demnach 


na nia 
sd 1 1 B+ As | I pees alr 2”) PS re) 
w= : 12 > Y 
sae ea pay ELE cag Pak 


E 


2A, A. eee aes — "Fey +412) 
oa ag 


Oe pee es 
sin — usin — x. 

a a 
Durch Einsetzen dieses Ausdruckes in (21) und durch Differentiation nach u und v sowie unter Be- 


riicksichtigung der Identitaten (7) bekommt man folgende geschlossene Ausdriicke fiir die Sin- 


1 R. Snidaric, Dissertation T. H. Wien 1944. 


4 


sees ee eee ar ‘ 
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gularitaten der Kriimmungen und der Verwindung: 


oe ees 1 I(y—v)  L(y—») 
Gi anKyes ill & i 
AFR [LO+% , Lhy+»] . 244 
taal te | ty a Ha) + ay +A, 
e2gVl 1 


a ~ ga Kyat 4 bY —9) + ho —9) 


23 — 23 21, A. 
+53 TH A, L(y +2) +A, L(y +) ja 59 [AS L(A, y+Ag v) +2? I(Ag +A0)]}, 
Se 2 
gudv ~ gn K Pel 


23 + 28 Dhl 
Figg +) + Gay +0) — pF, Vay +he0) tC lay +40) 


(22) 


[G0 —9) +G0—2] 


Wenn der Aufpunkt auf dem freien Rand liegt (v = 0), entartet die Singularitat der Kriimmung 
der x-Richtung. Es wird dann 


VII 
a2s 2 Pea + A) 1) 1 pe eee: i 
tut leno Sa KOLA TARTAR 2H) oS eerie 2) 


Fir kleine Werte von x — | gilt naherungsweise die Formel 


azgvll 3 asyt 1 “+1 B Sus 
Ons 9 =055 ee On ly = 0 ! aK (1 + 2x) (/ 2 i)(4 facets : 3) =) 
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die von Dworzak} ermittelte Singularitat (6 = 1), des isotropen Halbplattenstreifens ist. 

Da die Singularitat im Aufpunkt nicht endlich bleibt, ist der Regularteil in der Nahe dessen zu 
bestimmen. Wir gehen wieder von der klassischen Singularitat von Stein fiir die Kriimmung am 
freien Rande einer unendlich ausgedehnten Halbplatte aus. Diese lautet 
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~ Demnach wird im Aufpunkt 
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Zu dem verbliebenen Glied des rechten Teiles der Gleichung (25) gelangt man ebenfalls, wenn man 
in (26) den Grenziibergang fiir a > co vornimmt. 


6. Schlu&bemerkung. Bei der Aufstellung von geschlossenen Ausdriicken fiir die Singularitaten 
der Kriimmungen und der Verwindung frei dehnbarer orthotropen Halbplattenstreifen kénnen 
grundsatzlich die fiir die isotropen Platten geltenden Verfahren angewendet werden. Ebenfalls 
laBt die sich von Pucher? angegebene Lésung des isotropen eingespannten Plattenstreifens auf den 
orthotropen erweitern. Dieser Fall, soll in einer spateren Arbeit behandelt werden. AuBerdem 


1 W. Dworzac, Osterr. Ing.-Arch. 1 (1946), S. 66. 
2 A, Pucher, Federhofer-Girkmann-Festschrift, Wien 1950. 
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bereitet keine Schwierigkeit die Behandlung der zweiseitig gelagerten rechteckigen Platte. In 
diesem Falle wird zu der Singularitat des Halbstreifens das regulare Integral ; 
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hinzugefiigt werden. Die Konstanten A, B sind aus den Randbedingungen entlang der Querrander 
y = 0 und y = b zu bestimmen. 

Da die Singularitaten vom Kennwert x abhangig sind, wird im allgemeinen erforderlich sein, 
in jedem Einzelfalle die Berechnungen von Grund auf durchzufiihren. Selbstverstandlich ist auch 
die Méglichkeit vorhanden, fiir ausgezeichnete Werte von x eine Reihe von EinfluSfeldern in Vor- 
rat zu haben; beriicksichtigt man aber die veranderliche Lage des Aufpunktes, so erscheint dies 
dann unpraktisch. 

Jedoch ist es im Falle reeller Wurzeln A (iiberdrillsteife Platte) méglich, fiir alle Aufpunkte 
einer Paralelle zur y-Achse und alle Arten von Querrandern mit einem einzigen EinfluBfeld 
auszukommen. Allerdings muS man dann eine umstandlichere Auswertung des EinfluBfeldes in 

Kauf nehmen. Das Verfahren stiitzt sich auf die Tatsache, 


da®B die Einflu8funktionen der Kriimmungen lineare Aus- 
(m) + fathe t : driicke einer Einzigen Funktion sind. Betrachten wir 
eae z.B. die Kriimmung 0?SY!/du? (19). Wir zeichnen das 
AB EinfluBfeld der Funktion 
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Pine welche der Geraden (6,0) zugeordnet und unabhangig 
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Ganz ahnliche Beziehungen gelten fiir die Funktionen J,. Liegt nun die Last P im Belastungs- 
punkt (x, y) , so ist der Wert der Kriimmung in der x-Richtung laut (19) gleich zu 
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Die Auswertung des Einflu®feldes erfolgt, indem man auBer dem Punkt I auch in den Punk- 


ten II bis VI mit der Last P belastet und die EinfluBwerte mit den entsprechenden Koeffizienten . 


multipliziert (Abb. 5). Fiir Kennwerte in der Nahe der Einheit (x ~ 1) versagt das Verfahren. In 
diesem Falle kénnen aber die nach den Potenzen von x entwickelten Formen angewendet werden. 
(Eingegangen am 19. August 1959.) 
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